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A User's Guide to ARSPIQ: The Univariate 


(| Time Series Aanlysis Program for Autoregressive : 
a Spectral Information Quantile Identification 
" 1. Introduction 


ae - The ARSPIQ (AutoRegressive SPectral Information Quantile 
identification) program is a modified version of ARSPID, a 
univariate time series program in the TIMESBOARD Computing 
Library, (Newton, 1983). ARSPIQ provides various diagnostics in 
the time and frequency domains to help one determine whether a 
time series is long, short, or no memory. 

ARSPIQ is written in FORTRAN and consists of a main program 
and 58 subprograms, many of which are contained in the TIMESBOARD 
Computing Library. Many of the subprograms used in the quantile 
analysis are modified versions of those used by the ONESAM 
program, (Parzen and Anderson, 1980). The current version of 
ARSPIQ was developed on an Amdahl 470V/6-II and an Amdahl 470V/8 
operated by the Data Processing Center of Texas A&M University. 


The version of FORTRAN used should be compatible with most 


FORTRAN IV or FORTRAN 77 compilers. 

The basic goal of ARSPIQ is to provide diagnostics to aid 
in modeling a univariate time series. ARSPIQ is not intended 
to be a modeling or forecasting program, however. The models 
produced by ARSPIQ are intended only as guides to formulate more 
complete, rigorous, or valid models. Gur approach is to run 


ARSPIQ to obtain useful model building diagnostics and then to 


Pe ela ee a rc ee Baraat 


employ ARSPID in the more formal model building stage, i.e., 


estimating parameters and checking residuals. With this goal 

in mind, ARSPIQ has been made as fully automatic as possible 
with "logical" defaults provided when input options are specified 
to be zero. 

For a more complete discussion of the statistical functions 
mentioned in the next section, see Parzen (1981) and Parzen 
(1983). Also, for a complete discussion of the TIMESBOARD 
computing package and useful time series results, see Newton 


(1983). 


2. Stages in an Analysis 

2.1 Quantile Analysis 

The original time series Y(t), the periodogram E(w), the 
correlogram Pp lv), and the partial autocorrelations a (m) may 
be treated as data batches and exposed to a quantile analysis 


to gain further insight about these quantities. For a data batch 


{X(t), t=l,...,T}, one defines the informative quantile function 
IQ(u) by 
X(t) - X 
50 2t-1 
TOU) ote ee ’ = fs eg ’ 
2X75 - X95 2T, 


where X5o is the median and X55 and Xyo5 are the upper and lower 
quartiles respectively. Plots of IQ(u) are useful in describing a 


data baten and are primarily used in this work to diagnose probability 


characteristics. In ARSPIQ, the informative quantile function 
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for the original time series is used in the remaining analyses 
since it serves to ‘normalize and detrend"” the series 
analogously to subtracting the mean and dividing by the standard 
deviation. 


For a specified null value for the density-quantile function 


fQ(u) = £(Q(u)) where f is a probability density and Q is the 
corresponding quantile function, and for a raw estimate of the 
quantile density q defined to be the derivative of Q, one 


computes the weighted spacings 
d(u) = £,Q,(u) q(u)/o, 
where 
o = ft £.Q (u) au) du 
fe) ooo) 


The weighted spacings d(u) hence comprise a probability density 


and By is treated as an estimate of scale (compare to standard 


deviation). ARPSIQ takes 


-1 
£,Q,(u) = (@7>(u)) 


for the normal case where ¢ is the standard normal pdf and ie 


is the standard normal quantile function; and 


£,Q,(u) = l-u 


for the exponential case. When a data batch comes from the 
specified null distribution, d(u) oscillates about a uniform 


value of 1 and the cumulative weighted spacings 
D(u) = f¥ d(v) dv, Oxucl 


oscillates about the identity line. ARSPIQ provides a plot of 
D(u) as a goodness-of-fit check. For independent normal data, 
the sample autocorrelations and partial autocorrelations will 
also be normally distributed while the raw periodogram will 
have an exponential distribution, hence the normal and 
exponential cases are automatically checked where appropriate 
py ARSPIQ. 

A quantile analysis also produces some useful descriptive 
statistics, but discussion of these will be withheld until the 


appropriate time. 


2.2 Frequency Domain Analysis 


The raw periodogram £7 (w) is computed as a Fourier Transform 


a of the time series and then normalized to integrate to one. 
a Specifically, 

ee 7 A 

$ f£,(@) = (AL + BEO7/T +, kel,...,Q 
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where o is the sample variance of the informative quantile of 


the data and A. and BL are the Fourier coefficients evaluated 


at frequency k/Q. (A, is the amplitude of the cosine term and 
a By is the amplitude of the sine term at frequency k/Q when a 
sinusoid is fit to the data set. See Newton (1983).) As 
mentioned previously, white noise data produces a periodogram 
that has an exponential distribution, and hence a quantile 
analysis of E(w) uses £,Q, (w) = l-u. 

As discussed in Parzen (1983), the key descriptive statistics 

for the raw periodogram are the median, variance, and spectral 
log range (or dynamic range) defined to be 


max min 
SPLR = w_ log £7 (w) - w log £7 (w) 


Large values of periodogram variance or SPLR indicate moderate 
to long memory, while small values of periodogram median have 
the same interpretation. For normal white noise, the median 
is log 2 = .693, the variance is one, and SPLR = 0. 

A nonparametric periodogram is also computed and is 
labeled "Local Quantile Periodogram" on the output. For a 
neighborhood length NQMP, a batch of periodogram ordinates of 
length NQMP produces a lower quartile, median, and upper quartile 
value. These values are labeled G25, G50, and G75 and are 
associated with the median frequency in the neighborhood. G50 
is interpreted as a nonparametric smoothed estimate of the 
spectral density while G25 and G75 serve as "confidence limits." 
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The log of these values is plotted. 
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A Parzen window estimator using truncation point 


M = {largest even integer < (T/2)} also provides a smoothed 


periodogram. This estimate of a spectral density is used in 
creating the cepstral correlations described later. Parametric 
spectral estimators are also computed and will be described 
later. These estimators will have SPLR computed as well as a 


quantity known as the delta memory function defined by 


k : 
} log £ (ED = Gog eee) 


1 
eK, 0 © 


J 
where m = 0 and m = Q/XSEAS with Q = NFREQS and XSEAS specified 
by the user. By default, NFREQS is chosen large enough for the 
FFT routine and XSEAS = 12. These two delta sequences are 
plotted on a scale from -3 to 3 and have apparent memory 


classification properties as discussed in Parzen (1983). 


2.3 Time Domain Analysis 

Besides the interpretation provided by Box and Jenkins 
(1970), the sample autocorrelation function (acf.) and sample 
partial autocorrelation function (p.acf.) provide added insight 
when viewed in the quantile domain. 

The autocorrelation function (correlogram) Pp (v) is 
computed as the Fourier transform of the raw periodogran. 


Besides the plot of Pp Cv) and its quantile function, one also 


examines the value of 
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Small values of this quantity indicate no memory or short memory, 

while values larger than say 0.1 typically indicate long memory. 
The partial autocorrelation function is given by the value 

of a, (m) estimated via the Yule-Walker equations for the general 


AR(m) model denoted by 
Y(t) + a 61) Y(t-1) +...+ am) Y(t-m) = e(t) 


where e(t) is white noise and a, (im) estimates a (m) . 

A second partial autocorrelation function is computed 
using the Burg algorithm which appears to provide better 
estimates than the Yule-Walker equations when roots of the 
characteristic polynomials are near the unit circle (near 
nonstationarity). 

Only the Yule-Walker estimates are exposed to a quantile 
analysis with normal £52: The Burg estimate3 are primarily 
intended to produce an alternate approximating AR model for 
ARMA select modeling purposes and spectral analysis. Of interst 
is the number of values of. the IQ function for the p.acf. that 


are greater than one in absolute value. Outliers in both the 


acf. and p.acf. data batches indicate moderate to long memory. 


2.4 Diagnostic Modeling 


Five parametric models are provided to give insight into 
the nature of the time series examined. Srectral densities are 
computed for four of the models and residual variances are 
stated for all five. In all cases the parametric models are 
actually treated as approximating models and hence are really 
nonparametric in nature. This suggests that ARSPIQ should not 
be used to provide THE model for a time series although one of 
the five models examined may be considered appropriate by some 
criterion. 

The first two models considered are approximating auto- 
regressive models with orders given by the first and second 
relative minima of the CAT function. Values of AIC and log 
a residual variance are also provided in this AR modeling stage. 

For the best order model, a spectral analysis is carried out 


with computation of SPLR and 8 and display of the various 


spectral quantities. Only the coefficients and residual variance 
(RVAR) are produced for the second best order model. 


A second modeling stage actually produces five mcdels, 


‘ but only three are examined in detail. This second stage is 
primarily intended to suggest subset ARMA models to help better 

e understand the nature of the series. In many cases the ARMA 

i models are merely AR models mimicking the previous CAT derived 
ones, and hence little new information may be provided at this 

e stage. 


For subset ARMA modeling, the default strategy is based 
upon the Burg sample p.acf., but a second strategy based upon 
cepstral correlations may be requested by the user. The p.acf. 
obtained via the Burg algorithm is fed to a recursion algorithm 
which produces AR(p) coefficients where p is chosen to be the 
larger of the best and second best order as determined by CAT. 
Inverse autocorrelations are also computed and displayed. A 
spectral analysis is then performed on this AR model. The 
model is then inverted to an infinite MA and truncated to order 
M (NCOVM) corresponding to the number of computed autocorrelations. 
The residual variance and coefficients of this model are input 
to a select regression routine to produce an appropriate 
covariance matrix. Along with cross-correlation values, a value 
called PVH is also displayed. This is analogous to the 
prediction variance horizon defined in Parzen (1981) but is 


actually 1-PVH(h). For the infinite MA model 
Y(t) = e(t) + By (t-1) + Boe(t-2) +... 
we define 
PVH(h) = of (1 + 6, +...+ 8) 


where a? is the normalized mean square prediction error 


(innovation variance) estimated by RVAR for the Burg AR model. 
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The horizon HOR is the smallest value of h for which 
PVH(h) > .95. As a rule of thumb, HOR = 0 ‘~vlies no memory, | 
HOR/ (AR order) > 4 implies long memory, end values in between 
have specific interpretations as suggested by Parzen (1981). 


The select regression procedure automatically picks a diaenostic 


ARMA model displaying various criterion values, RVAR, and 
spectral density. 

A second select strategy is based upon obtaining truncated 
infinite MA coefficients directly from the cepstral correlations. 
The cepstral correlations are essentially the Fourier transform 
values of log spectral density. We define cepstral correlations 


by 


1 


5 log fu) e2tivv dw 


v(v) = f 
where f(w) is estimated by the Parzen window estimate of the 


spectral density mentioned earlier. A recursion algorithm 


provides the MA coefficients and may be expressed as 
n 
(n+1) 8(m7l) = Lo (k+1) y(k+l) 8(n-k). 
k= 


Truncating to order M, we proceed as before to the ARMA select 


stage. In addition, an estimate of c? based upon the formula 


log o2 = {5 log £(w) du 


is provided using the fact that the right hand side is actually 


¥(0) and hence 


oZ = exp [y(0)] 


This value is labeled "SIGMA INFINITY SQUARED" on the output. 


Note that the values of RVAR may also be treated as 


estimates of ce For the cepstral correlation truncated MA 


model, the residual variance is 


~ 2 2 
RVAR 1/(1l + By Fe ak Bap - 
The AR residual variances are obtained from the Yule-Walker 
equations, and the ARMA residual variances are obtained from 
the suitable element in the covariance matrix after a series 


of SWEEP operations has been performed. 


Input Options 

One of the primary goals of ARSPIQ is to make a complete 
memory analysis of a time series as automatic as possible. 
Hence, when a zero is specified several options will default to 
"acceptable" values. The input to ARSPIQ is described below. 


All values should be right justified. 
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3.1 First card of a run (in I5 Format) 


Cols Quantity /urpose 


1=5 NANS Number of data sets to be analyzed 


3.2 Second card of a run [in (615, 2F5.9) Format] 


Cols Quantity Purpose 
1-5 NTAPE File number where data set to be 


analyzed resides (e.g., if NTAPE = 11, 
specify FT11F991 DD card describing 
data file) Default: NTAPE = 5, data 
follows card 2 of input options. 
6-10 IRWND l if NTAPE to be rewound before being 

read, 0 otherwise. Default: IRWND = 0. 

11-15 NCOVM Number of autocorrelations to compute. 
NCOVM < 250. Default: NCOVM = 
(sample size)/2 

16-20 NFREQS Number of equally spaced frequencies at 
which to compute spectral density. : 
1152 > NFREYS > (sample size) + NCOVM. 
Default: Determined by subroutine FNFREQ 

21-25 IOPTMX 2 if select ARMA modeling desired, 3 is 
ARMA modeling using cepstral correlations : 
desired, 4 for both, 1 otherwise. 


Default: IOPTMX = 2. 


Cols Quantity 
26-30 NQMP 
31-35 XSEAS 
36-40 —-REXP 
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Purpose 


Length of neighborhood in computation 
of local quantile periodogram. If no 
local quantile analysis is desired, 
specify NQMP = -l. Default: 

NQMP = (sample size)/XSEAS 

Seasonal period, e.g., 12 for monthly 
data. Default: XSEAS = 12. 

Power or log transformation. 

Y(t) replaced by Y(t)**REXP if 

QO REXP>1. If REXP = @. or 1. on input, 
no transformation is performed. If 


REXP 1., a log transformation is produced. 


3.3 ARSPIQ data sets 


ARSPIQ employs ARSPID data set conventions. All time series 


data sets must be stored in standard format as described below. 


Card 1: A data set title (anywhere in the 80 columns) 


Card 2: Columns 1-5: 
Columns 10-29: 


Sample size right justified 


A FORTRAN format statement describing 
how the data has been entered on the 


remaining cards, e.g., (8F19.@). 


Card 3,4,...: The data entered in the format specified on Card 2. 


3.4 JCL for executing ARSPIQ 
// Job Card 
//* JES3 Control Cards 
//PROCLIB DD DSN=USR.R579.TW.PROCLIB, DISP-SHR 
// EXEC ARSPIQ 
//SYSIN DD * 
--input option cards -- 


--//FTnnF991 cards describing input files. 


4. Sample Output 
Following is an abbreviated output from ARSPIQ analyzing 
the Wolfer Sunspot numbers. The JCL that generated this run is 


given by: 


//ARSFIQ JOB (R5791006C2S60710> TW)» “WOOLIFIELD’ »MSGLEVEL=(190) 
//*FORMAT FRe DDNAME=» DEST=XEROX » JDE=JFMT i sFORMS=11100 
//¥LEVEL 1 ; 

//PROCLIB [Lt DSN=USR.RS79, TW. FROCLIBs ISP=SHR 

// EXEC ARSFIQ 

kay Dn x 


13 0 0 0 4 O 11. 1.0 
//FT13F001 Dl ISN=USR.RS79. TW. TSDATACWOLFER) » DISF=SHRrLABEL=( 995 IN) 
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7 f) c PROGRAM ARSPIQG (ARSON) ARSPIO 
é 2 c ARBPIO 
ie 3 Cearvece eeseaes eesene eeneeues eeseece aasPlg 
s c ARSPIO 
ry c ARSPID SPECTRAL BNALYSIS VERSION USING CUANTILES and AREPIO 
: e 4 INFORMATION FUNCTIONALS. (cooe 1, wury 1983) AREPIO 
bh ? c AASPTO 
» ry € TH1$ PROGRAM CARRIES OUT THE MEMORY ANALYSIS OF A&A TIME SERIES ARSPIO 
‘ ® € AS DESCRIBED IN PARTZEN {DEPARTMENT OF STATISTICS TECH REPORT ARSPIO 
c WO 0-22, Texas M UNIVERSITY. JULY 1983} & QUANTILE ARSPIO 
€ BNALYSIS 18 CARRIED CUT POR THE FOLLOWING DATA BATCHES BRBPIO 
t 1 THE aw DATA aAsPie 
c 2 THE SAMPLE SPECTRAL OENSITY (PERTOFOCRAM) AnsPle 
c P tHe CoO tocaam AaRSPIO 
i4 a THE TIAL AUTOCORRELATIONS aRsPig 
¢ ARSPIO 
c THE CURRENT OPTIONS ane: ARSPIO 
c NTAPE - TA NUMBER WHERE CATS RESTOES IN “ARSPIO FORMAT” aaseio 
c a Tape 18 TO BE REWOUND, © OTHERWISE AASPIO 
c “TOVM - NUMBER OF COVARIANCES TO BE COMPUTED ARSPLO 
C NFREOS - NUMBER OF FREQUENCIES POR PERIDDOCAAM [>*NONCOVM} ARSPIO 
t }OP TMK - 2 TF SELECT ARAMA MODELING DESIRED, 3 IF ARMA ARSPIO 
c MODELING USING CEPSTRAL CORRELATIONS DESIRED, BRSPIO 
c @ 17 BOTH METHODS DESIRED, 1 CTMERWISE ARSPTO 
(3 HOMP - LENGTH OF METCHBORHOOD IN LOCAL CUANTILE ANALYSIS ARSPIO 
c OF PERIODOCAAM fe-1 ITF NO LOCAL QUANT. ANALYSIS) ARSPIO 
c XSEQS - SEASONAL PERIOO, 1.8.. 12 FOR MONTHLY DATA ARSPIO 
c REXP - POWER TRANSFORMATION. IF REXP GT. 1, LOG TRANS. ARSPIO 
c ARSPIO 
€ OEFAULT INPUT OPTIONS [PLACE A ZERO IW COLUMN FIVE OF DPTION ARSPIO 
c CORD): NTAPECS TRWNOLO wCOVMEN/2 NFREOS FOUND BY FNFREO ARSPIO 
c MSEOSH12 lor rmx? NOMPEN/XSERS REXPS1.0 ARSPIO 
¢ SPECIFVIWG & ZERO FOR ANY OPTION PRODUCES CEFAULT VALUE. ARSPIO 
c ARSPIO 
c ARMA MODELING PERFORMED USING INFINITE MA REPRESENTATION OF ARSPIO 
c MIGHM ORDER AR MODEL [ORDERSMAX(BEST OROER, 2NO BEST ORDER)). ARSPIO 
c BURC‘S ALGORITHM USEO TO CENEMATE AR COEFFICIENTS. ARSPIO 
€ WLTH FOPTMKCS OR 4, CEPSTRAL CORRELATIONS USED TO FORM TRUNCATED ARSPIO 
¢ INFINITE ORDER MA TC PRODUCE COVARIANCE MATRIX FOR ARMA SELECT ARSP1O 
ic AWALYSIS ARSPIO 
c ARSPIO 
c PROGRAMMER: TERRY Uv. WOODFIELD ARSPIO 
c ARSPIO 
Cree eseavreereseseeeseseneeeuseeneseneeteseseauseseeaseuugesetsenese ARSPIC 
c ARSPIO 
COMMON /DAT/ YT( 1000) ,.CTL1200) ,$T11200} ,.FT( 1200) .wKM(61,63), ARSPIO 

a 1EAR{ SOO} A280), CORRIZ$O), INDT(6O} .COFTIS0}, ARSPTO 

a AIC{2$0) .ALAv(250), ARSPIO 

a S1G(2).MINVCI2). NAMESI 250), ARSPIO 

a LABV( 20}, 101110), ALPNTP( 280) ARSPIO 
COMMON /UNTT/ ITUNTT.NSCRCH ARSPIO 
OimenStow 1422.2) ARSPIO 
CIMENSION RCOEFFI6) X1(650,6) ,.wKeM1(6,.6),COVMAT(6, 6) ARSPIO 
DIMENSTON LiS).La8D( 20) ARSPIO 
Soiegeieieers eS esd ehetedereiars Saree See ras ARSPIO 

CATA LLZ/4HBEST, ,SM2ND , AHBEST/ : aRsPic 

DATA NIN/S/ ARS?1O 

DATA LABD/' ORIG’ "INBL'.* DATS ‘Sk +, 1be OF: ARSPIO 
PuMITs6@ ARSPIO 
wscacHs2 AAsPlo 
TwOP1s® eaTaNnt? oO) ARSPIO 
WALTEL LURIT, 300) ARSPIO 
WRITEL ITUNES T, 3S) ARSPIO 

S  FORMATI(//,T10, 814 1H8},/,T10,°8°.790,78#°,/, 710, aRSP10 
ote BRSP1G - AVTIORECRESSIVE SPEC TAAL INFORMATION SCUANTILE ‘. ARSPIO 

+ TOENTIFICATION’ T9O,°8' ./,710,°8°, TPO, 18°, /, THO, ‘*# EMANUEL °. ARSPIO 
*'PARZTEN, TERRY J WOOOFIELO, AND H. JOSEPH NEWTON’. T90,'0°,/, BRSPIO 

*T10 ORPARTMENT OF STATISTICS, TEXAS ASM UNIVERSITY, *, ARSPIO 

*T90 ‘,/,T10,"8 COLLEGE STATION, TEXAS 77843°,790,'#°./, aRsPro 
*T10, JULY 1983, 790, °R'  /,710,%8' 190, 8°, /, TIO, BIL IMs) ARSPTO 

c ARSPIO 
c INPUT WOMBER OF SERIES TO BE ANALYZED : ARSPIC 
c ARSPIO 
READ(NIN, 10) NSEYS ARSPTO 

10 «~FORMAT(1615)} arsPero 

c ARSPTO 
c ARSPIO 
c sTagmt LOOP - aRSPIO 
id ARSPIO 
c aansPio 
O06 1300 NSET:1.NSETS ansPig 

c ARSPLO 
c REHD OVERALL GPTIONS : aRsPlO 
c ARSPIO 
REACT HIN, 20) NTAPE, 1RWNO.NCOVM,NFREOS, IOP TMX NOMP,XSEDS .REKP ARSPIO 

20 FORmaTISi&8 32FS oO} anseigo 

c ARSPIO 
c Reap OATA SEY - ARSPIO 
c ARSPI1O 
IPT HTOPE BO O} NTAPEDS aRsPlg 

TEL TAWND £6.17} REWIND NIAPE ARSPIO 

CALL DATAIN(INTAPE, VT WYT,L.LaBY) ARSPTO 

c aansPig 
c S#T OEFAULT VALUES aRSPLO 
c ARSPIO 
IPL NCOVM EQ O) NCOVMe(WYT/4)02 ARSPIO 
TE{NCOVM GY. 280) NCOYMt#250 ARsPig 

TFPI MPREOS EQ ©} CALL FNFREOIMYT .NCOVM, NFREOS) BRSPIO 
TFURSEAS.£0.0 0} KSEaS112.0 aRrsPto 
LF(TOPten €9 6} LOM tTHKs2 aasPia 

16 (MOMP £O 0} NOMPSTPIRIFLOATINYT)/XSEAS) ARSP1O 

TET REXP £€O.0.0) REKPS1.0 ansPle 

le (reuPm CTY 1.0) GO To 100 aasPig 

LF RERP £6.4.6} CO TO 420 aRsPtg 

OO 9O fs1 MYT arsPio 

VTL ERY TEL pe sMe xP arsPio 

eo CONTINUE @RnsPig 
Go TO 120 aRSPIG 

100 CONTINVUE arnsPid 
Of 1tO Tet wet ARSPIO 
yr(rpeacoccvrTii)) ansPie 

1t0 Continue AAsPIe 
120 COWTINUE aARsPle 
WATTECLUMIT $601 LASy MTAPE, LRWNO NYT HCOVM HFRECS RSEAS NOMP, anePrio 

onene aaAsPrto 

180 PFORMAT(///,110,2086,/,120,°IMPUT TAPE © °.33,760,°'IRWND 8 °, ARSPIC 
OLD TO, LENGTH GF SEMTER © ° 86, /, TIS, “HEOYM © *° ES, THO, aRBPID 

* MPFREOS © *° 14,760, RSERS 8 ° PK. 2,/,.720, "ROMP «¢ *°,13, 780, aasPrg 
o°MERP © * eS) aAsPiro 

€ @RSPIO 
7 PERFORM GNESAM ANALYSIS OF TIME SERTES is aRnsPie 
c aRePig 


CALL QUENTINYT YT LaBY LABD. 1,0, MYT°4 Y2S,YSO. ¥78) aRSPIC 


e 


ed al Doll te ee a i cee dee 


ee ee ee gM OR Le ep 


€ 
¢ DETRENC BY BSUSTRACTING MECSAN AND "WORAMALIZE™ BY DIVIOING 
€ By TWICE THE 10 RANGE. 
¢ 
weueTer. /(2.¢1Y98-¥26)) 
OO 180 101, MYT 
100 VIET PCIVTETD-YeO}evmMULT 
c 
¢ CeScATee Y-TiLDA ; 
c 


246 CONTINUE 
IPUiNCOVM NE, O)CALL DESYTINYT NFREOS ESEAS NCOVM, NOMP, 
SLL2,10P 1M8) 
IP UNSET WE NSETS) WAIYECIUNIT, 300) 
Joo FORMATI IAT) 


¢ 
c *rNiSm LOOP 
€ 
t) 


3200) «CONTINUE 
WRITE(IUNTT, 6) 
STOP 


OUTIME AICLEIAY, M,N. ALAV,AIC,LAIC) 


(4 
¢ 
¢ SUBROUTINE TO CALCULATE AKAIKES AIC CRITERS ON AICIK}® 
rq ALOGIRVIKRDIO2K/N AND ALSO ALAV(KDSALOGIAVIK)},KOd,.. LM. 
¢ FURTRER, LADCOK THAT MINIMIZES AIC 18 OETERMINED. 

ra 

¢ 

€ 


Cm meee cones ene HeR ae ee ner eats eeeedeaesase sees eesasessassssenes 


OIMENSTOW R¥(M),AICiM)  ALRVIM) 


OnNeNn 

00 10 Ket,™ 

ALAV(KPPALOGCIAVIK)) 

BICUK CARVE KOE 2 ePLOaT(K))/0ON) 
1g CONTINUE 

CALL MIMIAIC,M,CO.LASC) 

t#(CO.GT.0.0) taiczo 


AETURN 
end 
SUBROUTIAE AMA(AL PH, NP LNG, OETA) 


Core veeccnseneesasesanenscenteeeaesaeceanesssesseeenseaee 


t 
¢ PunPase + TO FING THE COEFFICIENTS BETA(1)...., 
¢ GETAINO} OF & MOVING AVERACE REPRESENTATION 
(4 OF ORDER WO CIVEN THE AUTOREGRESSIVE 
c co TCTENTS ALPH(1),..., ALPHUNP) 
c 
c ARGUMENTS ALPM - INPUT VECTOR OF LENGTH NP CONTAINING THE 
cf COEFFICIENTS OF THE AUTOREGRESSIVE PROCESS 
c OF OROER WP 
(4 ORDER OF THE AUTOREGRESSIVE PROCESS 
c QUTPUT VECTOR OF LENGTH NO CONTAINING THE 
c COLFFICIENTS OF THE MOVING AVERAGE PROCESS 
c OF ORDER WO 
(4 Lr) + ORMOER OF THE MOVING AVERAGE PROCESS 
c 
4 RECO. TIMESBSOARD SUBROUTINES - NONE 
4 
Coesvsen verse sasssesus eenes seeessezasasas 
c 
OUMENSION ALPHINP) ,BETAINO) 
¢ 
Seva rpe-acen(s) 
TF(mo,.€0.1) RETURN 
©0 1 112,m0 
cero. 
170) .GT.0P) GO TO 2 
Ce-aLPmiy) 
2 MemMINOUNP, 1-1) 
0o 3 Jttlm 
B] CrC-ALPml uy) sOETAlI- J) 
1 weTarryec 
€ 
RETURN 
eno 


SUBROUTINE ARSP(ALCPH, SIG, WP NFREOS ,CT,SPOC) 


Cee meccece eens ss eBoeasseceeezrsunoseeesseesessesenneeees 


c 
c SUBROUTINE TO CALCULATE THE AUTOREGRESSIVE SPECTRAL 
c CENSITY (AT TRE MFREOS FOUMLLY SPACED FREQUENCIES 
¢ SE TWEEN © AND TWOP1) FOR THE AUTOREGRESSIVE PARAMETERS 
c 2 ALPH( WP), AND SIG (RESIDUAL VARIANCE). 
© 
c inpur 
c mp MLPMt dd... ALPH(NP) SIG, NFREOS 
c 
c QuTPuT : 
c SPECTI),.,.. SPECI NFREOS) 
c 
€ MUNTLL TARY 
c cr 
c 
c SUBROUTINES CALLED : FFT 
c 
Coscacee eeace eenneuvccneresen eeeene 
c 
DIMENSION ALPM( 1), CT(NFREOS) .SPECINFREOS) 
c 
FAC TSIC/I@.cAaTAN(t.0)) 
te (me. GT.0} GO TO 10 
00 & 191, NFREOS 
& sPecti)sFac 
CO 10 98 
10 CONTINUE 
c 
po ft Trt, wrREes 
crtipeo. 
1 SPECII}e0. 
ctiijyes. 
DO 2 fe1ne 
2 CTC1etpearPotie 
c 
CALL PRTICT, SPEC. WFREOS NOREOS ,WEREOS, 1} 
c 
DO 3 164, mFREOS 
JD SPECTIbAFAC/ICTCA)oCTELssPect(iiesPectiiy 
c 
$9 CONTINUE 
RETURN 
eno 
SUBROUTINE CEPSTC(FT WEREOS NCOVM,CT,S7.8156) 
Ceaces sooee eee @ecces eee eaees 
c 
¢ SUBROUTINE TO CALCULATE THE CEPSTRAL CORRELATIONS GIVEN 
c THE SPECTRAL DENSITY FUNCTION FTI1),. |... 7T( PROS). 
c 


ansPlce 
aRsPic 
aasPla 
aasPio 


ansPia 
ARSPIO 
aasPig 
ARSPIO 
aasPlo 
ARSPIO 
ARSPIO 
AaRSPICO 
ARSPIO 
ARSPIO 
ARSPIO 
ARSP1Q 
ARSPIO 


a3 


ced rRtovencias 
OETWEEN © Auo Two] 


ImPUT NFREOR - WUMBER GF EOUALLY 


SuUBFPROGMAMS CALLED: FORTIER 


6 

c 

rd er: CTRAL OENSITY FUMCTION 

e uCOVM - NUMER CEPETHAL CORMMALATIONS DESIRES 

c 

€ OuTPur: CT - CREPSTMAL CORRELATIONS : 
e S1G - SIGMA IHPIMITY SOUARED FROM INTEGRAL LOC SPECTRUM 74 

c 

€ auNliL lary sf arsPio 

€ ansPlo | 
c 

¢ 

£ 


POeeeee ee RS UEP CeCe eee Or reer re re rercerseerrrrrrrrr rrr rr errr ety) 
OIMENS 10m FT(NPREOS) CTINFREOS) STI NFAROS) 
COMPLEX £11200) 

FAT. /SFLOATINGREOS) 
DO 10 161 wFREOS 
cr¢ibeacocivrir)} 
tO STIDPePLOaATII per Ac 
NOVO teNCOvMes 
CALL FORTERICT ST. NFREOS A. NCP) 
06 20 Tet wCovmM 
20 CTD PsReartai ety) 
SiCeemP(meari acid) 
RETUAN 
END 
SUBROUTINE CHIP[IOPT.CHISO} 
OIMENSION CHISOl62),CH1162) 
18(}OPY £O.1) GO TO to 
1et1oP1 @€06 21 CO TO 20 
TFL 1GCFT £0 3) CO TO 30 
1O CONTINUE 
DATA CHI/6 64,9.227,11.32,13.28, 
1168 O9,16 Bt, 


.27.92,29. 
40.31.01. 
21,53 
71,64 
ga. 76. 
9t.87. 


ansPio 


20 


-31,23 
-20,34. 
78 


.398,62.03,63 
174,03,72.18.73. 
co TO 99 

30 CONTINUE 
Data CHI/1.64 


16$.73,66 81,67. 8 -97,70.00,71.92/ 
99 CONTINUE 

oc 40 121,62 
40 CHISO(T)FCHIEI) 

RETURN 

eno 


QUTINE CLPLTOIX.N. INIT, WAME, XMAX, KXMIN) 


Crsees eaesese aeee suave ate eseseaes eanene 
c 
c SUBROUTINE TO PRINT GND PRINTER PLOT THE N-VECTOR XK 
c THIS ROUTINE IS & MODIFIED VERSION OF CLOLTI DESIGNED TO 
c PLOT THE OELTA MEMORY FUNCTION 
c 
c INPUT 
c nix ARsPlo 
c art PRINTED JNOEN OF FIRST PRINTED © ARSPIO 
c WAME - @ CHARACTER LITERAL CONSTANT GIVING ARSPTO 
¢ LABEL FoR x ARSPTO 
c MMAN.XMIM > MAX ANO MIN VALUES FOR SCALING PLOT ARSPIO 
id RSPIO 
t SUBROUTINES CALLED - NONE arsPiog 
¢ ARSPIO 
CeVeasecerseeeseseseurtasetnerseeoeeghaenseesseseteessseseenest ARSPIC 
c aRSPIo 
OIMENSION imp. AL( 101) ARSPIO 
Date NOUT/E/ ARSPIO 
OBTA BLANK, DOV, 2, SL/IM (IM. THe LIMES aasPlo 
c 
1orprzgr16o q 
iv-iw GY t) Go TO 10 
WRITE MOUT tt) wWAmME BET) 
WY PORMAT{1On AG, '{ti} © °F ts. a) 
co ro 98 
10 CONTINUE 
c 
¢ ImMETia@rize ar 
c 
poms GT 
Owelmm- 13/2 
00 20 vei, mm 
20 « aL(ylsoor 
WRITE NOUT. 26) MAME. (AL( UT, utt, MM) 
ZS FORMATI/, 140 161, 6R, 84/F0K, 1S IM), 2K rat} 
DO JO srt mm 
JO ALI) BL am 
c 
c FIND SCALING vatues 
is 
ReeRMAN- NMI M 
KzeROr 3? 
VERA LT 18-20} LOOP TAT 
c 
t PLO? 
c 
PPEe ews 
oo 40 Jr}, 
ir( torts @0 +) GO TO 36 
Creel el sh - mms) sRN 
cre2 efCr- 8) 
co TO 37 
a6 cre uJ 
37 meQOmelCrer pore 
MUU RTERO) CSE arsPig 
Muemde? aRsPlo : 
weiTeiwour.3 JS MESS, GACCTD Ded eee ARSPTO 
38 FORMAT(10xK,18.F10.4,2K,10ta1} aRsPlg 
JJtsset aRsPioe 
ML (CR) SOL am ARSPIO 
40) «Continue ARSPIO 
€ arnsPto 
os COMTINUE : aAsPIoO 
RETURN ARSPIO 


eno arsPric 


BUOMOUTIME CAPLTEI RM, THIT, WOME MM, SORT) 


Cece veseae ser er sarees eee sss eeaeeeeSSSASEROSOREESESSESHESOEES 


c 
c SUGROUTINE TO PRINT AND PRINTER PLOT YHE M-VECTOR XK. 
€ 
¢ Input 
c hon 
€ IMIT 1 PRINTED 1uDex GF FIRST PRINTED 
ce WAME +: @ CMAMACTER LITERAL CONSTANT GIVING 
i4 Rael. PER 4 
c mm: oNUMe OF COLUMNS IM PLOT ILE. 101) 
c torr ¢ 1,2 [PStwT on BAR PLOT) 
i4 
¢ SUGROUTINES CALLED - Man MIN 
¢ 
Ceseeeesnseerencaee eeanes eeneese eaeene esee 
id 
OIMENSTOW MIM), ALL 101) 
a DATA NOUT/E/ 
bs Oats Bi ANR. DOT. F/tmH . Im ,IHOs 
© 
1oPTESO 
ite ct +) GO To 10 
WAITE(WOUT. 91) mame, NLT 
Tt FORMATION 6,540} 8°, 
co TO @8 
10 «CONTINUE 
¢ 
c PmrTIMLrre at 
4 
Owetmm- 11/2 
DO 270 yet mM 
20 «aLisieoor 
7 WRITE(NOUT, 28) mame CaLivs ued, mm) 
‘ 38 POommaAT( / 10K. IML, 6". A6/10n TELIM-}, 2K. 105a1) 
‘J 00 30 set, mem 
30 ALi J) BL awe 
if 
= c f1ND Man BNO MIM 
€ 
CALL MAX(R  RMAR. SND) 
CALL MINCE MME TNO! 
Abe RMAN- HMI 
etiam 27. ds. B20) LOOP THE 
c 
c PLOT 
ct 
Pret Bas 
he 00 40 vei.n 
IPC 3OPTR.€0.1) GO TO 36 
ee CreCML SD MMIND/RX 
y cilez2.0(6%-.8) 
co 168 37 
= ae cio 
327 meOne(Ctier, Jot.6 
acme? 
te(}0e7 £O.t) GO to 38 = 
oe oo 39 151.K 
~ 39 alti? 
- 328 CONTINUE 
c WAITECNOUT, 38) JI.KEUP act), ret, mm) 
= 3e FORMAT(ION,IS.F10 4.2K, 10181) 
af duesuet 
ALU KD OBL ANK 
tri 1OPT.€0.3) GO To 40 
8O 41 FIL K 
a @cfl) +e. anK 
60) «CONTINUE 
t 
99 CONTINUE 
' RETURN 
. END 
SUBROUTINE CLPLTZIK, VN, INIT. NAMEX, NAMEY MM, 1OPT) 
Ceoeeeesaenes ee seeesPoneaatsseeuasesetesseaeteaatesessseeeese 
c 
c SUBROUTINE TO PRINT AND PAINTER PLOY THE W-VECTORS K,Y¥. 
c 
c INPUT 
c N,M,Y 
¢ IMTT > PRINTEO INOEX OF FIRST PRINTED XK AND Y 
c NAMEN ,NAMEY | ACHARACTER LITERAL CONSTANTS 
c Givincg c.4BEL FOR & AND Y¥ 
¢ lal WUMBER OF COLUMMS IH PLOT (LO.91) 
c 1OP7 : 1,2 {POINT OR BAR PLOTS) 
c 
c SUOMOUTINES CALLED . MAX, MEN 
c 
ce aeotesn seoes anaes eraesen eenecen eesce 
c 
ioe OCIMENSIONW MINT. Vim) aL en) 
77 OBTA NOUT. BLANK, OOT, 23, 22/6,96 1M. Ime, MKS 
c 
ir tm GY.1) GO TO 10 
WRITEINOUT, 11) MAMER, NET) WAMEY, Y(t) 
1) FORMAT(ION, 6,701) 8 5 P98 8.6K, ae Et) 8 ee 8) 
+ co TO 99 
i 10 «CONTINUE 
- c 
“ c Iwitrarcize ab 
c 
= FOPTEO 
. sOrTY90 


Lect mms2pet 
Ont (mm- 1) /2 
po 20 sel, mm 
20 «alL(sls50T 
WRITE(NOUT. 26) NAMEN WAMEY (ALIS), Jet, mm) 
28 FOMMATI/, 16m, INT, 6x, 04, 68 BOSTON, 2610M-), 20. 81A1) 
DO JO set, one 
30) AL Ct OL aK 


. tc 
a c FIMO MAX AND MIN 
c 
: CALL MANNIE N, NMAR THO) 
_ CALL MENCN,6.NMIN THO) 
CALL MARCY, YMA INO) 
7 CMLL MINTY. RL YMIM, TPO) 
RReKMAR- EMI 
RYDYMAR-YMEM 
IFEAY LT $.6-20) [oPtHs) 
IPtmy CT §.8-20) 1ERTVes 
. c 
. £ PLor 
c 


ssermit 
Do 6s vrt lm 
iF (1oPte. £O 1) GO TO 36 
CreC (eC ud-wmimeysmay-t 
co ro 37 

38 Cie- 8 

37 «tF(ToPTy.@0 1} GO TO 38 
C2civisi-vmrMels/Ry 


AASFIC 
aaselio 


aARSPIO 
aaseig 
aRSPIC 
ARSPTO 
aAsPlg 
ARSPI1O 
@RSP10 
ARSPIO 
aRSPIQ 
ARSPIO 
ansero 
aasPig 
aRSPTO 
ARSPIO 
anselo 
aRSPI1O 
ARSPIO 
arses 
ARSPIC 
ARSPIC 
aRsPi1o 
ARSPIO 
ARSPIO 


75 


oe 


el 


s 


co Te 3¢ 

38 cae 6 

320 «Kteomel(eton. por.s 
M2sOMe (CII. por.s 
acimipest 
ALUK2 De 22 
IviloPy.@0.1) GO To a1 
Oo 42 let. ut 

2  aeetpors 
Co 43 relt,K2 

ay ae(rper2 

41) CONTINUE 
PRIFBINOUT. 40) Ju. xiv). vis). 

420) FORMAT(10K,18,2F10 4,20,0141) 
Juegset 
BLIRD OBL ANK 
AL (K2) BLANK 
s#(10PT. £0.41) GO TO 48 
00 44 1es,K1 

a6 ALI) BLANK 
OO 46 Tedt Kn 

oe ALI I) BLANK 

46 CONTINUE 


Ce ee 


c 
¢ 
98 CONTINUE 
RETURN 
eno 


SUBROUTINE CPMAI(FT NFREOS .NCOVM,CT ST OETA, R¥MA, SIG) 


Ceewcenvecanasesssesseenececesssescesessareretesensszevasesesens 


SUBROUTINE TO COMPUTE CEPSTRAL CORRELATIONS FROM SPECTRAL 
DENSITY FUNCTION FT VIA CEPSTC AND THEN TO COMPUTE WCOVM 
BETAS FOR GN INFINITE MOVING AVERAGE MODEL TRUNCATED TO 
OROER wCOVM. 


IMPUT Ft - SPECTRAL OENSITY FUNCTION EVALUATED AY NFREOS 
EQUALLY SPACED FREQUENCIES BETWEEN © AND TwO P1. 
NFRECS - WUMBER OF FREQUENCIES 
NCOVM - NUMBER OF CEPSTRAL CORRELATIONS TO COMPUTE AND - 
GCRDER OF TRUNCATED MOVING AVERAGE MOCOEL 


OUTPUT: CY - CEPSTRAL CORRELATIONS 
BETA - VECTOR GF MOVING AVERAGE COEFFICIENTS 
RVMA - RESTOUAL VARIANCE FOR M&A MODEL 
$31G - SIGMA INFINITY SOUARED FROM INTEGRAL LOG SPECTRUM 


AUXILLIARY: ST 
SUBPROGCAAMS CALLED: CEPSTC,FFY,CLPLTS 
Pewee ee eRe ROR TER SH SHES TEREST RES ROTHER EEE HEE H DER SE EE EeTes EEE eSeEBESE 
common /UNIT/ IUNIT, NSCRCH 


CIMENSION FTL 1), CT(1),S709), METAL) 
CALL CEPSTCIFT NFREOS .WCOVM,CT,.ST,SI1G) 


AnAnaAnnnAnAANnAANAANAAANA 


¢ CT CONTAINS CEPSTRAL CORRELATIONS 
WRITE{IUNIT, 10) SIG 
10 FORMAT(/, 10M, °SIGMA INFINITY SOUAREO {VIA SMOOTHED PER.) « °‘, 
Te UOne Rui lERi Ne DuMed CeneePE: aie Vi 
: COMPUTE THE BETAS 
c 


BETAiI)eCT{( 1) 
DO 823 172, NCOVM 
BETAI1I)+#¢. 
Temreg-1 
00 820 Kr1, IMI 
420 PETAL IL eBETALTIOCTIK)sBETALIT-K)eFLOAT(K) 
828 BETALI} © BETACT)/FLOAT(L} © CTr{T) 
aviary 
90 830 1't,wCOvM 
830 RVMACRVYMACBETAL I) *BETALI) 
RvMAr1 /RYMA 
WRITELTUNIT, 632) RYMA 
832 FORMAT(//,10%,‘'MA MOOCEL VIA CEPSTRAL CORR., RVARt', 
oFTO S$,//,3ORK, FIRST 10 COEFFICIENTS OF INFINITE MA: ‘,//) 
CMLL CLPLTI( BETA, 10,1, @HBETA, 64,1) 
RETURN 
end 
SUBROUTINE CUMSP(SPEC .WFREOS. 1OPT .wK,NI)} 


Coecennavornranesesesessesesscecssescerceenesnssesesssese 


SUBROUTINE TO PRINTER PLOT CUM SPEC DIST FCTN 


1MPUT 
WFREOS ,SPECII}...., SPECINFREOS) 
1? WFREOS.LT.O, CUM SPECTRA 15 CALCULATED 
eur NOT PLOTTED. 
TOPT : 1 MEANS FREOS ARE © TO 281 
© MEANS FREOS ARE © TO PI 


OuTPuUT 
wkit), . WRENN) WHERE Nie(NFREQOS/2)¢1 IF 
JOPT+1 OR WFREOS IF 1OPT#O : CUM DIST 


Peeve nenessrereeresasessunecearnessereeseeeueeseuaesesses 


ananaannaannnanna 


DIMENSION SPEC( 1) weit} 


an 


WPROCNFREOS 
WFREOS: [ABS (NFREOS} 
PFREO?0 6 
Wis (MFREOS/2Z}01 
DELF ee eATAN(' OO) /FLOAT(INFREOS) 
TF(t0oPT £0.19) GO TO 10 
NICNFREOS 
ORL FsB eATAN( 1.0) /FLOAT(2°(NFREOS-1)) 
10 CONTINUE 
cespecit) 
wat tdespecty) 
00 20 12,1 
wKEEDswa(l-1besPpEectt) 
cecespecti) 
20 «CONTINUE 
BO JO Eet.mt 
wutipewniiise 
30 «CONTINUE 


IF{erag 1.t.0)] GO TO 40 
CALL SPPLTI WK, WT, PPREO, DELP. O} 
ao MFREOSSNFAD 


RETURN 
eno 
SUBROUTINE CYARPI CORR, 4 NP wk AL PM 


° 
aRsPiCO 
aanselo 


76 


ARSPIO 
arnsPie 
aRSPIO 
ARSPIO 
ARSPIO 
aRrsere 
arsric 

aASPSO 
aRSPTO 
ARSPIO 
ARSPIO 
arsPig 
ARSPIO 
ARSPIO 
ARsPtO 
ARSPIO 
ARSP1OQ 
ARSP1O 
ARSPIOC 
ARSPIO 
ARSFIO 
ARSPIO 
ARSPIO 
ARSPIO 
ARSPIO 
ARSPIO 
ARSP1OQ 
ARSPIO 
ARSPIO 
ARSPIO 
ARSP10 
ARSPIO 
anseio 
ARSPIO 
ARSPIO 
ARSPIO 
aRSP10 
ARSPIO 
ARSPIO 
ARSPIO 
QaRSPIO 
ARSPIO 
ARSPIO 
ARSPIO 
ARSPIO 
ARSPIO 
ansPio 
ARSPIO 
ASPIO 
ARSPIO 
ARSPIO 
ARSPIO 
ARSPTO 
ARSP1O 
ARSPIO 
ARSPIO 
ARSPIO 
ARSP1O 
ARSPIO 
ARSPIQ 
ARSPIO 
ARSP1O 
ARSPIO 
anseta 
ansPia 
ARSPIO 
ARSPIO 
ARSPIQ 
ARSP1O 
ARSPIQ 
ARSPIO 
ARSPIO 
ARSPIO 
ARSPIO 
ARSP1Q 
ARSPIO 
ARSPIO 
ARSP1Q 
ARSPIO 
ARSPIQ 


ARSPIOQ 
ARSPIO 
aRSPto 
arsPig 
aRSPIO 
arsPig 

asPro 
sP10 
aRnsPrag 
ARSPTO 
aRSP1Q 
aRSPIG 
aRSPTO 
ARSPIO 
AASPTO 
ARSPIOQ 
aasPle 
ARSPIO 


NnAnnAAAAn An AA AANA 


a 


aAnannnanannanannaa 


a 


AHANAANHAAANNNRAAANANANARARAKH AAR ANA NANANAAAAN 


0 


27 
20 


Cause 


SUBROUTINE TO CALCULATE THE COEFFICIENTS ALP HII)... 
ALPMI NP) OF Bw AUTORECAESSi ve PROCESS OF OADER NO 


Givew COMRIF},...,CORRIM (me CE we) 
inpur 

a> mm COMAITI.... ,CORRIM) 
ovreur 

MLPNI Td, .. ake ae) 


SUBROUTINES CALLED - WONE 


Been seers eonsese eS Penedsazacarssaseessseesssesessans® 
DIMENSION CORRIM) .WKI( NP) AL PMI NP) 


weittbe-coam{s) 
Ifine GY +) GO TO 10 
BLPmL dd ewntin) 
RETURN 


DO 20 192,NP 

C+eCORRI I) 

czes, 

Imved-t 
oO 28 vet, Imi 
CHeCrownr( sss CORR I -y) 
CPeCaewKil( ype CORal sy} 

ca-qerse2y 

AcPH{ ise 
OO 26 Jet, tm 
BLPMISbewK Its) *CowKitt-J) 
BO 27 se1,3 
WKIT SD ALP WIS) 

CONTINUE 


RETURN 
END 
SUGROUTINE CY¥CRIA,AO.M} 
OIMENSTON R{M) 
oo t t1.m 
RIUPeRITISRO 
RETURM 
end 
SUBROUTINE CVPARTIRA, AO.M WKS wK2, PART) 
EPEC ORTECOOCEPOREOCrererererrr er errr reer errr te etter 


SUBROUTINE TO CALCULATE THE PAATIAL AUTOCORRELATIONS 


PART(V)..... PARTIM) GIVEN THE AUTOCOVARIANCES RO, 
Berd, 2. ROD, 
twPUT 
WMO MITTS... RIM) 
ouTPuT 
PART( V1... .  PARTEM) 
SUBROUTINES CALLED ° NONE 


. sevasceausencete 


CIMENSION Ri mM). WK Mm) .wKI( mM) PART (mM) 


wert rys-ait) 
PARTE) © -Alt) 
Fim 806.3} CO FTO OD 
CO 10 ft2,m 
Crismer) 
Cazey 
featet- 4 
OO 18 vet, amt 
CrsCrowmr(ypemit-s) 
CIeCPommr( ypeRis) 
PMat( yds fessor) 
wear) Paari is 
OO 16 vet, im 


16 WHIZ SI eWRTETISPART I Tewari t-y) 


00 27 vei.t 


27 wail slewezis} 
"0 CONTINUE 


OO 100 Tet. 


100) MC TPMT PERO 


RETURN 

exo 

SUBMOUTINE CV¥SPWIR.RO,1OPT NTRUNC,NFREO. WORK! WORK2. SPEC) 
Chvemw MOM TIL. ..., RUIMTRUNC), THE FIRST NTRUNC OD 
AUTOCOVARIANCES OF A TIME SERIES KI.}, SUBROUTINE CY¥SPW 
CavpeoucarTes §& ClID,- ., SPECI NFREQ), & SMOOTHED ESTIMATE 
Or Tme SPECT Lk DENSITY OF KI.) AT THE WFREO POINTS 
Oo TWOP L/MFREO, ..., TWOPTSINFREO-1)/NFREQ, TME WEIGHTING 


FUNCTION 3S SPECIFSED BY THE USER BY 1OP7T (SEE BELOW) 


METHOD © THE CHOSEN WEIGHTING FUNCTION WORKII1),... 
WORKTINTRUNC) 1S CALCULATED FOR A SPECIFIED TRUNCATION 
POIMT WIRUNT. THEN S VECTOR OF LENGTH NFREO IS 
CETEAMINED AS (RO, WORKH( 1 FeMmti),... WORKI(NTRUNE )e 
RINTAUNC).O, -.0) 

Tee SUBROUTINE FFT IS USED TO FINO THE COSINE 
TRANSFORM OF THIS VECTOR THIS TRANSFORM IS DIVICED BY 
TWOPL CIVING THE CESTAED sPectaum 


TmeuTt 
MIRUMC, MFREO; RO,RITE,... , RIMTRUNC) 


1oeT 


aS GIVEN BY T w ANDERSON © THE STAT ANAL OF TIME 
SEMIES, THE FOLLOWING OPTIONS OF WEIGMTING FUNCTIONS 
AE AVAILABLE 


TOPrey eaatierr 

1ORPTr? MODIFIZO GARTLETT 
torts) NANNING 

1OPYTes © MAMMIAG 

TOPTeR - PARTEN{1) 


LOPTes aatenis) (NTAUNC MUST BE EVEV) 
OutTeur - 

sPecty). SPECI NFREO) 

AURILL TARY WORK! WORK? 


SUBROUTINES CALLED wry 


anerico 
aasPic 


aRSPto 


ARSPIC 
aASPIC 
ARSPIO 
aASPIO 
aRsPto 
aRSPIO 
arsPig 
aASPIO 
BREPIO 
ARSPIO 
aREPIO 


AREPIO 
aRSPIiC 
AaRsPio 
ARSPIO 
ARSPIO 
aRSPIO 
ARSPIO 
ARSPIO 
ARSPI1O 
QRSPIO 
ARSPIO 
aRsPilo 
aRsPio 
aRSPIO 
ARSPIO 
aAsPig 
ARSPIC 
@asPiog 
ARSPIO 
ARSPIO 
aasPic 
AASPIC 
arsPrico 

RSPIO 
ARSPIO 
asPlo 


aRBP ICO 


AaASPIO 
ARSPIO 
AaRsPIo 
ARSPIO 
aRsPio 
ARSPIO 


asPig 
ARSPIQ 
ARSPIO 
ARSPIO 
aARSP10 
aasPico 
ARSPI1O 
ARSPICO 
SSPIO 
~ SP10 
ARSPIO 


ARSPIO 


CIMENSION RINTRUNC) WORKT(NTRUNE)  WOMKIINFREO) SPECT INF REOT 


Imitsactze 


ann 


Pied eaTanit oO} 
TwOP1e2 oP) 
oo t bet. weRreo 
BPEC(I}e0 

1 wORK2Z[ 1} 00. 
seecit)erno 
Crs mTAUNC 
CTemFREO 


CALCULATE SPECIFIED WEIGHTs 


ann 


IFC 1OPT #0. 
1F(1OFT £0 
teLtoryY a9. 
tri tor 
1F{10°T @€o. 
IPE LORY €O. 
IPL IOP £9. 
we proP 
co To $$ 


co To 2 
co to @ 
eo T0 6 
co To 6 
Go 19 10 
co TO 12 
Go TO 20 
co TO 22 


Seveaaun- 


BAATLETT 


anan 


2 00 2 vet NTRUNC 
Oreg 
ozeor1/ct 
O3+01/C2 
DB woarKisfule2 sift -O2)/11. -03)) 
co To 14 


MODIFIED BARTLETT 


ann 


4 pO $ tt NTRUNC 
BS woanttsdez eft -CFLoarisds01)) 
co ro 14 


HANN ING 


ann 


6 00 FF Jet, NTAUNE 
Fo wWORKIC ID) eCOS((PrleFrLaar(vd)d/Cs) 
co TO w4 


HAMMING 


nan 


a 00 9 Jtt,. NTRUNC 
BP WORKIG I} 81 O8% SIT9COSLIPIAFLOSTIVI)I/E1) 
co TO 14 


PARIEN( I? 


ann 


tO =DO Ft vt, WTRUNC 
PY WORKICUEt2 811 -EFLOaTIVI/Ct)ee2) 
Go TO 14 


PaRrEen( ul 


nnn 


$20 MUENTRUNE/2 
OO 13 ver.mt 
Jteseoms 
OrsFcoOarty)/c) 
O2eO1e92 
Oleoissy 
woRKi(vlez of 
12 WORKIN (G1)2a 8 
Go To 14 
20 OO 21 Jt1,NTRUNC 
21 WORKI( UD 2.O/ 1 1 OO1FLOATIVI/C1) O26) 
co ro 14 
22 OO 22 ser. wteRuNc 
23 woRKi(s)22 O/(1 OF fFLOATI VY) /C1} 988) 


1.°6 80266 203) 
Cr. -CRROATIUNI/C1) Deed 


c Cai frT 


14 90 18 yet. wTRUNC 

1B SPECI Vet) ewORKI( IORI) 
CALL FFT(SPEC .WORK2 NFREQ NFREO MFREQ, 1) 
00 t6 srt, NFAEO 

16 SPECI UITSPECI UY) /TwOPs 

99 CONTINUE 
RETURN 


OUTINE DSTAININTAPE,X,N,1, LAB) 


Ceeres eesas ager seene sesenansee 
c SUBROUTINE 10 READ A DATA FILE FROM TAPE NTAPE AS FOLLOWS 
c 
c camotr : caert(t)...., Las(zo) (2000) 
c CaRO2 : SAMPLE SIZE NW, FORMAT LIt),....LIS).(15,6¥, 50a) 
c CMROD.CARDG.... : DATA RE1),..., EN) TN L FORMAT 
¢ 
Ceseres senpsece eeeges eeanae sertesaves 
c 

OMrmMENSTOW K(1),. 018) Laslzo} 
c 

REROINTAPE,1) LAB 

y PFORMAT(20~a8} 
REMOUCNTAPE, 23 WAL 
2 FORMAT(IS, 6x, 8846) 

Read wrare Ly CEE) ret) 
c 

RETURN 

eno 

SUOPOUTINE DESTAT(N, NW. NAME, IMEAO.L, INIT, 10UT, 028,080,078) 
Cesces seeee tennn eecuseresseans Seeorccee euees 
c SUBROUTINE TO PRINT CROERED ARRAY BY QUAKTILES ANO COMPUTE 
¢ CESCRIPTIVE StaTISTICS 
c iepur: 
c X- MARAY OF OCROER STATISTICS 
c N- DOIMENSTON OF ARRAY Xx 
c NAME WaME OF DATA SEY. MUST SE ARRAY OF DIMENSION 20 16 
c CALLING PROGRAM 
c TUNITT. MUM@ER OF UNIT GUTPUT IS CESIREO GH. 
4 1m17 © FOR FIRST CMLL, 1 THEREAPTER 
c BOUT. 1 TF OUANTILES TO BE LISTED, 0 OTHERWISE 
t wk WK? WORK VECTORS OF LENGTH WNIsH8203 
c OuTPUT PRINTED OUTPUT IS ON 1UMET, 
c MO SUBROUTINES CALLED 
Creeeveessersevessecnsasteessseseetenetesendersessrenascetestasesesetensensenes 


COMMON SUNLT/ 1UNITT RSTREM 
DIMENSTON Kim). MAMEL 20), SUM( 4), SUMSO( 4) 
CIMENSION Lia) 
DIMENSION IMEAD( ZO} 
mumit + TUNIT 
c COMPUTE L, THE ARRAY OF QUARTILE SIZES 
seCIMIT £6 OF GOTO & 
tTeCcl@t .€9 1 Guta 25 
8 ator mse 
kioeoan 


SP ee 


ARSPIO 


AREPIO 


78 


eat. 
pez. 
ee) 


962. 
63. 


ges. 


c , 


30 


38 
c c 


40 


$0 


c t 
ct i 


vo 


yet 
1002 
1003 


10048 
1008 
VOC 
1007 
1008 
1009 
3010 
yor 
1012 
rors 


yore 


1020 


o 


AMAR ANRANAANAAANAANKN AAA NAAN ANA AAANAANAN 


2 4 ue 

Poeun 

teeeean 

1ORAw 6 MODIN,&) « 1 
COTO 120,14.92,13), 18k aM 
CONTINUE 

r@ieauredt 

Go tO 20 

CONTINUE 

ktoeure ft 

ee 

co Ta 20 

CONTINUE 

22 8 hee 1 

doe ee 
t@ eaupe d 

CONTINUE 

Livbous 

REPITeLLs) © b2 

Liddertay ¢ 4&3 

chavoutay * La 
RINT DATA BRRAY - OWE COLUMN FOR ESCH QUARTER. 
WRITE(MUMIT. 1001) NAME 

WRITEC NUNIT, 1020) 1HEAD 

IPC 1OuUT.£0.0) COTO 38 

WALTEL NUNIT, 1002) 

WAITE( NUNIT, 1003} 

OO JO Te s.ae 

@MITELNUNTT, 1008) TimePd mCat sd) © FILME e Zoe. MERC a © 1) 
WRITEL NUNIT, 1008) 


TF(L1 .6T. BL) WRITEINUNIT, 1006) xELE1}) 
VRLL2 .GT LL? WRITEEMUNIT, 1007) wie (2)) 
TRF(LY GY Lt} WRITET NUNIT, 1008) K(LI3)} 
LEiLG@ GT LED WRITE NUNTT 1009) CA. K(LI4)) 
TPLINST €O, 1] RETURN 
OMPLTE ANC PRINT DESCRIPTIVE STATISTICS 
Ke) d 

so 

$30 50 

po $0 1 + 1.4 

$12 0 0 

Sol » 0.0 

eR 8 Ltt) 

00 40 Jos KKK 

Btoe Sr ¢ KES) 

SO] + SOI © RlJ}ex(y) 

CONTINUE 

Koaveow utd) 

Ss © S$ ¢ $I 

$50 + $80 * SO? 

Sumit} «© $2 

SUMSOLI) & SO! 

CONTINUE 

JF(TOUT.€0.0) GO To 8S 

WRITE( NUNIT, TOTO) (sum {yy,323,8) 
WRITECNUNIT, TO1ST (SUMSO (1),271,4) 


K@AR «© S$/FLOATIN} 

VAR «© ([S$SQ0 -S*RBAR}/FLOAT(N-1) 

$D + SORT( VAR) 
W THE ORIGINAL VERSION ©2858 .0530,078 ARE COMPUTED MERE. 
N SRSP1QO, SUB. OCVUENT COMPUTES THESE VALUES. 


O10 © OF8 - O26 
3BGRIG » (KBAR 
so10 © so / 
S$O1OLG © ar OGi 


- 950) / (2. © O10) 
(2 
so 
SSONISSO/FPLOATIN 
12 
13 
14 


a) 


WRITE( NUNIT, 10 
WRITE( NUNIT, tO 
WAITEL NUNIT, tO 
cOnTinue 
RETURN 
FoAaAmMAT(///T20, 2006) 

FORMAT(/T40, ORDER STAHTISTICS IM QUARTERS ‘'/T40,28{ tHe)! 
FoamatT(1T20,°' SEQUENCE’ /T20,° WITHIN’ /T20,° QUARTILE Wes 

oF ERST OUBRTER SECOND OUARTER THIRD QUARTER FOURTH QUARTER’ 
* /TIS OC ime de, ZEIKHL LIC Med 2K, sah oMe dd) 

FORMAT(T20,18 . 411K. F1S 41) 

FORMATI IX} 

FORMOT[ Ime, T29.FIS 4) 

FORMAT( Ihe TAS. F IS 4) 

TORMAT( IMe TE1LF IS 4) 


1 
J 
i} 
to W,028.080,0768,010 

} SSOM, NEAR, VAR, S50, XBARIO.SOIO, SDIOLG 


FORMAT( Ne T2O,18 T77,F 18.4) 

FORMATIL INO T2O0,°SUM: Be, 15K. 718. a)) 

FORMAT(1HO T20.°SUM OF '/T20,° SOUARES’ IH. STIX, FITS. 4}) 
FOMMATI///T4S, “DESCRIPTIVE STATISTICS’ /T4S,2311me)) 
FORMAT(//,120, SAMPLE’. T28,° LOwEeR' ,152,° upper = °, 
« Tee, INT _ OUARTL'/,T20,° SIZE’ 

° .728.° QUARTILE’. 140,° MEDIAN’, 162.° QUARTILE’, 
’ T64,* RANGE’. //, 720,18, 729. 81611 4, 9K)) 
FORMAT(/// TIS. "SUMSO/KH' , 130, MEAN’ 740,‘ VARIANCE’, 782, 
. “ STO DEV’ 166, "MEAN 10°,T76.‘STO DEV 10°,TBB, 
. “LOG STO 10°//,786,71GH1.4, 9K) 

FORMAT(T20, 2084) 

eno 


SUBROUTINE DESYTINYT NFREOS _XS8AS,NCOVYM, NOMP, 
@.r2, lor rex) 


teen 


VOROUTIMNE TO DESCRIBE Y-TILOA VIA PERTOCOGCRAM, AUTOCORRELATIONS. 


ONESAM ANALYSIS OF PERTODOGCRAM, CORARLOGRAM, OND PARTIAL 
UTOCOMAELATIONS 18 PERFORMED TESTING FOR WHITE MOISE. 
mPUT: MYT - NUMBER OF OBSERVATIONS IN Y-TILOA 

MFREOS - NUMBER OF FOQUALLY SPACED FREQUENCIES AT 
WHICH TO CHTAIN SPECTRA 
SEAS - SEASONAL PERIOD OF DATA 


NCOYM - “UMBER OF COVARIANCES TO COMPUTE 
ROMP + SITE OF WEDGCMBOCRMOGD FOR LOCAL QUANTILE 
ANALYSIS OF PERI OOCOCRAM 
a7 MAMES OF BEST, SECENO BEST AR MODELS 
1OPTmMn - SELECT MODELING OPTION 
1 - MO BAMB SPLECT MODELING 
20 BAMA SELECT USING AR PROM BURG ROUTINE. 
AA INVERTED TO LARGE OROER MA TO PRODUCE 
COVARIANCE MATRIX FOR SELECT 
3 - BAMA SELECT USTMS TFRUNCATED INFINITE MA 
DERIVED USING CEPSTRAL CORRELATIONS. 
@ - @0TR 2 aw 3 PERFORMED 


ureur VamiOUS MOCELS AMO DIACwOSTICS. 
MODELS. BR GY BEST CAT OROER 

OA 8: SECOND BEST CAT CROER 

Ah @y BURG ALCGCRITHM USING LARCHER OF BEST AND 
SECOND Best OROERS 

MA BY SHVERTING BUAC AA 

ARMA Br USIMC COVARIANCES PROM ABOVE MA IN 
SELECT MPECRESSION ROUTING 

Ma By CEPSTRAL CORRELATIONS 

ARMA Gy USING COVSRT ANCES OF CEPSTRAL MA IN 
SELECT RECRESS ICH ROUTINE 


ansPso 


aRnsPieo 
aasPrig 
BaSPIO 
BASPIO 
ARSPIO 
ARSPICO 
AaRSPIO 
ARSPIO 
AASPI1O 
ARSPIO 
aRSPIG 
ARSPIO 
ARSPIO 
ORSPIO 
BRSPIO 
aRsPid 
ARSPIQ 
aRsPis 
arsPrid 
ARSPIO 
ARSPIO 
ARSPIO 
ARSPIO 
ARSPIO 
aaseia 
ARSPIO 
ARSP10 
ARSPIO 
aRsPlo 
BRSPIO 
ARSPIO 
ARSP1O 
ARSPIO 
ARSPIO 
ARSPIO 
ARSPIO 
ARSP1O 
BRSPIO 
ARSPI1O 
AASPIC 
BRSPIO 
ARSPIO 
ARSPIO 
ARSPIO 
ARSPIO 


@RSPIO 
aRSPIO 


ARSPIC 


BRSPIO 
arsPig 


79 


610 


625 


140 
wa 


BSUOPROCRAMS CALLED: CTCVF CYCM, OUEMT DESTAT 40 
minman, QUITR FOFNED OTOFO, WEP 
PRPLOTC . FTEAP Min MA OREPS 
PCOAR CYPART CLPLTS. TRvia 
RELMN CVYBRP ARSP OMEMAY, MAXS 
PARTAR MACY CLOLT2, ARMA, MAME 
musrse 


eee esocesccecs eae 
COMMOW /OAT/ VYTI1000),C1T11200) ST( 1200) ,7T1120 
DERALEOCO) AIL 280) COMRI 280}. INOT(8O} COF 
BICI280) acMvi2zso). 
SICC7I) Minvel 2) waMes( zo}, 
LAB 17O} 1C1110), ALPMTP( 280) 
COMMON SUNIT, 1UNIT NSCACH 
O'MENSTON BETSL( 260) 1aBPA(20). . ABARZ( 20} 
OimemsiOw 662(2.2) 1 88P 120). LaOli 70) LAasaRi( 20 
OaTa Ps Row PERI: ,°OD0C:.” 1 Si 
oata cy CORR," 8.0C', "mam -. 
oata PAs PERT’ “LAL”, 
CATA LABARI/ BEST. 
ry Tensi. 
DATA LABARZ/ SPEC »' OBN 'SITV',* 
o'R MO’, ‘OEL ‘,1te 
OATA LABW/'SMOO' THEO’. ° PER’. *1000', ‘GRAM’, 
. "Mow NDOWS tie” “7 
Oa2TR MOIM/61/ 
TwOP S68 OF ATAN|1 ©) 


COMPUTE PERI COOCRAM COMRELATIONS 


CALL OTCVFIVY NYT NFREOS .NCOVM, CYST FT, M,RO) 
OLv¥e@O/ TWOP] 

Met WPMEOS/2Ie1 

DELP eTwORI/FLOCATINFEREOS) 

COLL CVCRIFT OLY NPREOS) 


PERFORM OWNESAM ANALYSIS ON Raw PERI COOGCRAM 
CALL QUENTINFREOS FT CaBy,k PLZ, V NFREOSSEa PR? 
LOCAL QUMNTILE ANALYSIS OF PERIOCOOCRAM 


171 NOMP.LE.O) GOTO 608 

NOLS ENFMEOSSNOMP)/2°10 

Olvort /FLOATINFREOS) 

eT(y)eD1vO 

DO 600 1:2.NOL 
CTETPCETIL-t}eo1ve 

CALL OLOCALIFT. CT NOL, NOMP,LABP} 
CONTINUE 


PRINT PERSOCOCGRAM AND CUMULATIVE PERIDOOCRAM 


CALC PREPSP(FT NFREOS,1,1.0,.CT, SPMIN, SPMAK) 
SPLRCSPMON-SPMIN 

WRITE(TUNIT. 610} SPMIN, SPMAN, RO, SPLR 
FORMOT[1M1, 130" 26NFOR PERIOCOOGRAM, SPLMIN © | F 
TOMSPLMAK © ,F12 &. 2K. 7HRIO) © .Ft2 &8,/,729,°SP 
FFREO'O O 

CALL SPRPLTIFT WI, FRREO. OELF, 13 

NFROCNFREOS 

WRITECIUNTT, 620) 

FORMAT(//. TOK 2OMCUMULGTIVE PERIGCDOCRAM : 

CALL CUMSPIFT NERO,1 CTT) 


PRINT CORRELATIONS 
WRITE (SUNT 628) 
FORMAT( 1m), T10 "TIME DOMAIN ANWOLYSIS- CORRELOE 
COLL PCORRIRA. RO.NCDYM CORR) 
PERFORM OQUANTILE ANALYSIS OF CORRELOCRAM 
CALL OUENT(NCOVM. CORR, LABY LABC.3,0,NCOVMS4 C2 
CREATE INFORMSTION LAC FUNCTION 
00 7oe8 1:1, NCOVM 
FT{T)+-O S*ALOG{) -CORR(T}*9CORMITI) 
WRITELIUNIT, 708) 
FORMAT(/. 10K. “INFOAMATION L4G FUNCTION: -O.S5*L 
CALL CLPLTICFT, NCOWM, 1, @MINFO, 41,9) 
walTe(tuNtt, 710) 
FORMAT( IM) /. 10K, 27HAR DESCRIPTION OF Y TILDA 
PARTIAL AUTOCOMRELATIONS VIG THE VULE WALKER 
CMLL CV¥PART(R.RO,NCOYM,CT,ST.FT) 


FT CONTAINS PARTIAL AUTOCORRELATLIONS 
PRINT PARTIAL AUTOCORRELATIONS 


CALL CLPLTOC FY. NCOVM,1,4HPACF,1.0,-1.0) 
PERFORM OQUANTILE ANALYSIS OF PARTIAL AUTOCORRELA 
CALL CUENTINCOVM FT LABY, LABPA,1,0,NCOVMOS,C2S 


COMPUTE PRECTCTION VARIANCES AWD AIC 


COLL PAaRTRAV( FT NCOVM, 4 0,87) 

CALL AICUSIST NCO¥M. MYT, BLRV, MIC, LAIC) 

WRITELS UNIT, 961) Late 

FORMAT(//, 10, 1SHOROER BY AIC © ,33) 

CALL CLPLTZ(BLRY, BIC. NCOVM, 1, @HOLAV, am AIC, 81, 


COMPUTE CAT AR ORDERS 


ofce1 sFLoarinvtTs 
CALL PARTIST. NCOVM, NYT #T NORD] 
WRITE(TUNTT, 970) 
rormat(s//) 
COLL RELMATFYT MCOVM DEC MIN! MmMINz) 
Ivirtimint) ce -) ) GOTO 731 
WRITELIUNIT, 1801} 

1 FOMMAT(//, TOR, ‘OROBER BY CAT © O.°) 
COmTinuE 


aR SPECT & ANMALYSIS OF ¥ 


WORD TRemint 

TOC EMIMteMINZ) EO O} CO TO 788 

RVTPCSTimini) 

TR(MINZ. £0.00) MINZEMIMI 

STIC( HD eST (orm? 

SIGCHI eS (MINZ) 

MIiMnvE( rT) emins 

MINVO (2 PeMIN? 

00 780 1OMbe1,2 

TFCLIORO £0.21 and (mint EO MINZ}) CO TO 786 
WAETECIUNIT. 725) MinvelIORNO) SIE i toro) 
FommaT(/. 100. ‘FOR OROREA *.12,° AR MODEL. RVAR 
WRITECTUNIT,. 730) LEZLV, TOROD L212, roORD) 


¥,050,7CODEA, 
ace o1ocm. 
CUMSP BPR;T, 
1CLS.PaARZ, 

OR. CVEPw, STARS, 
Si. SELMEC.Simx, 


sone. neserece 
O),wamiet. et), 
TISO}, 


}. Lawl 20) 
th 


TONES, 168" 


7 at oc, 
7M Bu’, RO a’, 


- w ARIE, 


5.P50,P76) 


12.5.2x, 
Ro © 752.8) 


Ram’ ,//) 


$.c80,.C781 


OG(1-RHOee2)°) 


2.4/7, 10K, 
EouaTions:') 


TIONS 


,€80,C78) 


i} 


1 F109) 


BRSPIO 
AREPIO 
aasrlo 
aansPlg 
aansPlo 
aansPeloe 
ARSPIO 
aasPlo 
AASPI1O 
aasPrio 
AREPIQ 
BASP IO 
BASPIO 
aRSP10 
ARSPIO 
ARSPIO 
BREPIO 
ARSPIO 
aRBPIO 
ARSPIO 
@RSPIO 
ARBPIO 
aRSPI1C 
ARBPIO 
ARSPTO 
ARSPIO 
arario 
ARSPIO 
aRrsPlo 
aRSP10 
ARSPIO 
aRsPlO 
ARSPIO 


330 FORMATS / LON VIMCOEPFICIE NTS POR . 284,12, THORDEA AREF IO 
} 3 aRerig 
c COMPUTE BR CORPFPICIENTS PROM CORRELATIONS aRSPiC 
c 


CALL CVARPI CORR. NCOVM, MINVC(IORO) ST, ALPHTP) 
CALE CLOLTILALPHTO MInvClIGRO}, +, 4nalPn. 61,1) 
HeL1ORO 26.2) CO TO T80 

SIG( LOMO) *SIG(1ORD}e TWoP) 


c 
c COMPUTE aR BPECTAAL DENSITY 
c 
CALL AMSPLALPHTP ,SIGEIORO) .MINVE{IORD) weAkOS. St FT) 
c 
t CREATE DELTA SEQURNCE FOR AR SPRCTRAL DENSITY 
c 


CALL OMEMAV(FY NFREOS, NSEAS, CT. ST. LABART} 
CALL PREPSP(FT NFREOS. 1.3 O.S87, SPMIN, SPMAX) 
SPL ROSPMAR-SPMIN 
WRITELIUNIT, 760) SPMIN SPMAK SOLA 
760 FORMAT(//,1OK, 1EHSPECTAA, SPLMIN © FIZ. 8, 2K, 
TEMSPL MAN 8 FIZ. Bae SPL R 8 P12 8) 
cai $s KVTIPYT M1, 7 FP meo pELF. 1) 
CALL MARSPALPT ALPHTP NFREOS  MIMYCIIORD) .100,1.8°6,NOPM, ST. TERR) 
WF ROUNPREOS 
WAITELIUMET FIOF LEZE1, 1ORO},LLZI2Z,10R0) 
FTO FORMAT(//, FOR, 1tMCUMULATIVE , 284,92, 1EHORDER SPECTRA :)} 
CALL CUMSPLPT NPRO.1,CT NIT 
THO) CONTINUE 
788 CONTINUE 


€ COMPUTE SMOOTHED FPERIODOCRAM USINT THE PARZEN WINDOW. 


MIMUNC INV T/&) 22 
tPEMTAUNS GT NCOVM) WYRUNC TI NCOVM/2) 02 
100 Twi 
COLL Cv¥SPwie MO LOP TW NYTRUNC MFREOS CY St FT) 
CALL CVCRIFY DIY MeREOSs) 
WRITEC TUNITY. Tea) NYTAUNC 
tA4 PORMAT(/ 108 _'TRUNCETION POINT FOR SMOOTHED PERITODOCRAM =» *.13) 


CMLATE DELTA SEQUENCE FOR SMOOTHED PERIODOGRAM 


noae 


COLL OMEMARVIFY WFREOS XSEaS,CT.ST,LABW) 
MINEO SCHEME SELECT PROCEDURE 


TOP TMK. 1 > WO ARMA SELECT MODELING 
2 > ARMA SELECT USING BR FROM BURG ROUTINE. 
AR INVERTED TO LARGE ORDER ME TO PRODUCE 
COVARIANCE MATAIK FOR SELECT 
3 + ARMA SELECT USING TRUNCSTED INFINITE MA 
DERIVED USING CEPSTRAL CORRELATIONS. 
& - BOTH 2 HO 3 PERFORMED. 


naneAanannan 


IFETOPYMR,@0.1} CO TO 
(Cuneo 
WRITELTUNIT, 600) 
800 FOMmMaAT( tet /, 10K, 29f Ime), /. SON, "MINED SCHEME SELECT PROCEDURE: , 

o/ VOR. 291 tme) Sf) 

lef lorrme. Ge.3) GO TO B10 

ICMmMa1oMRed 

co TO 838 


COMPUTE BETAS OF TRUNCATED INFINITE MA FROM CEPSTRAL 
CORRELATIONS 


aann 


@10 CALL CPOMALFY NFREOS NCCVM,CT,ST,BETA,RV¥MA,SICINE) 
co To 949 

O35 ITRUNCOMINVE[ 1) 
TOCTTRUNC LT MEMVCE 2S) TIRUNCOMERNC( 2) 
IFLITRUNC.89 0) GO TO B68 


COMPUTE PARTIAL AUTOCORRELGTIONS USING THE BURG ALCORITHM 
TMEN PROOUCE AR COEFFICIENTS FOR ORDER ITRUNC 


anan 


CALL GOYAMS( VT MYT ETRUNC. Z,CT ST OTRO, ATC ALRYV) 
WITYTELTTUMNTIT, B38} RO 
aye FOMMAT(/. TON, "FROM DTARB, ROt’ FIE.8) 


ry 


iC CONTAINS PAC F., ALRVY CONTAINS RESIDUAL VARIANCE. 
y.oT e TIAL AUTOCORRELATIONS. 


anoan 


AVMACALAVITTRUNE)/RO 
WRITE IUNTT, @26) 
836 PORMAT(/, 10R “PARTIAL AUTOCORRELATIONS VIA BURCS ALGORITHM: ',/) 
COLE CLPLTOL AIC TTRUNE. 1, 44PACR,1.0,-1.0) 
CALL PAATAR(AIC JTRUNE AL PMTP} 
CALL MACVIALPHTP, 1 O,TTRUNC RRO} 
COLL CVCATR, RO. 1TRUNC) 
WRITECIUNIT. 837) RO 
O27 FORMMATI//, 10K, ISHCORFPFICIENTS AND INVERSE CORRELATIONS 
+7, ON, JIMSUM OF SOUARES OF COEFFICIENTS *© .F10.4) 
CALL CEPLYT2Z(ALPHTP MT TRUNC, 1. AHALPM, ANCOR], 81,1) 
RVMAPTAYMAS TWOP 
CALL ARSPLALPHTP. RYMAP I TRUNC ,NFREOS ST FT) 
COLL OMEMRVIFT WEREOS XSEAS,CT. ST, LaBAR?2) 
CALL PREPSPLPT WFREOS,1,1.0,87,SPMIN,SPMAX ) 
SPUR CSPMAR-SPMIN 
WRITELTUNTT, 760) SPMIN, SPMAX, SPL 
CALL SPOLT{FT wt. FRReg. OELF, 1) 
COLL MANSPALFT ALPHTP, NFREOS, 1TRUNC, 100,1.8-6 NOFM,ST. TERR} 
NTROCWFREOS 
walTe(1UNss, 802) 
842 FOAMAT(// tOM, ‘CUMULATIVE SPECTRA FOR BURG AR:*) 
CALL CUMSPUFT wERo. 1.01. N1) 
CALL ARMATALPHTP ITRUNC,NCOVM, META) 
WASTE LUNTT 888) TMA 
a4e FORMAT(//, 10%, ‘MOVING SYERAGE MODEL VIA SUBR sama 
eF10 ©. //, 10M, RIMAST 10 COEFFICIENTS OF INFINITE M. 
CALL CAPL TSI BETA. 10,1, nF TA, 49,9) 
O4f TCHR ETCHKet 


RvaRs:, 
tf) 


c OETERMINE MAR AM AND MA OROERS FOR SELECT 


KyeMCOVM/2- 5 

17K ST JO} KIWZO 

W2ORD 

CaLt MACY(BETA, RY¥Ma,NCOWM,R, RO) 

CALL CVCRIR MO, KCOVM) 

COLL MAMTLL ER 1 O. BETA NCOYM,MOIM. RI, K2.0, NYT WKM, 
*KPS,NOW:S,COFT, INDT M¥MA, 1ER) 

Pave TAs Redo aesinaser 

'¥Sewk  5%v Vay) 

GLb SH 6s POT, COPT, ALPHTP, SETA, NORDAR, NORDMA) 
weital mt? 862) 

82 FoRmMaTi s/s, TON, ° TO ITHTERPREY CRITICHAL ADD ANO OBLOTE VALUES. °. 
© wOTe Te" 4,108 ‘THEY AME COMPUTED USING TwICE THE SamPre’. 
ee See Seem | 
eB TELL UNIT 880) AVS 

ase Pin. (44. VOR JIMMESTOUAL VARIANCE FROM SUBSET ARMA © ,F 10 7} 

To {RORDBR EO. CO AND WORDMA #20 0} CBO TH 880 
-7"MORDAR €O ©} GS TO B85 


81 


Po ee OT Ie SN OS PE a TT eta oe ge, oe ae Sey Lage! Re ON pa et en te ge he. We 


CALE CLPL TIL AL Punter BOROAR 1 Amar Om 61 1) 
obs TFUNOROmMA £6 01 C6 TO 
CALL CLPLTIIBETS aOADms 1 one 
860 COntiave 
Aveteay 
Cape sicorr 
Cale PA te? we 
SPL ReSPmaAn-SPmie 
WAITE(1UNIT FAO) SOmim EPman BF18 
c b SPPL TIFT wt eevee 
MFAOTMFREOS 
wRiTe(iueiy, 078) 
e776 Pommat( // 1OR  FIMCuMULaATiva Bus 


Ya @1 tt 


etwas! 


T,Ovbt ap 
e.1.1 @ 6 


@ct ef #Ty 
SPman) 


+ sama SPECTOR ' 


TAars CUMSPIFT weRO ' CY RI) 

TPL ECHR £0 1 AnD 1OP Ten we 3) CO TO 036 
Cantinu® 

RETuUan 

eno 


eid 


SUBROUTINE OMEMNY IFT we 


Cacascvessseenesenavesesossce 


os ES@asp cr Brie 


Bees eseresseeerare 


OUTIME YO COMPUTE DELTA MEMORY BOGCUENCE FOR BRPEC TAAL 
Ensif’ Fumcrion FT 


ImPur: FT - SPECTRAL CONSTI 
NFREOS - MUMBER OF ZOUSLLY SPACED FREQUENCIES USES Im 
COMPUTING FT 
NSOAS - SessOnmar. 1.46. §£ B , KSERSHH2 FOR MONTHLY OATA ETC 
VaBer - 6 Eu. GIVING YYPE OF ImweUT BPECTRAL DENSITY 


AUKILE TARY: CYT, BT 
CUTPUT: PRINTER PLOT OF DELTA SPOUENCE 


SURBPROGRAMS CALLED CePLyy 


annnnnanananannnann 


eeuacos evenes 
COMMON /UNIT/ TUNIT WSCROEM 
DIMENSION PTET) ,_CTII),ST0 1) LABEL (20) 
MSNFT1. /KSERS 
WRITRITUNIT, C00} LABEL. REMWF 
KRFREOIKSH FLOATI WF REOS) 
INrteo 
KFREOLCKFREOST 
19 CONTINUE 
CTC VP CALOCIFT( IME T+5)) 
sTcryecr(y) 
oo 20 182, KFREO! 
CUCHPCALOG( HTL LetmIts) 
20 STU S)eSTET-teectish 
oo JO 181,KPFREO 
STLTISTOULISFLOBT(ED-CY( Tet) 
IF{ST(S) GY 3.0) STITI*3 © 
IV(STEVY. LT. -3.0) STTI}*#-3 Oo 
Jo COMTINUE 
CALE CLOLTO(SY MFREO, 1, OMOELT.3.0,-3.0) 
IFC INIT WE.O) GO TO 40 
MIT EKERES 
co 70 10 
40 RETURN 
100 FORMAT(1H1,9m, CEL TA MEMORY FUNCTION: ,//, 10K, 2086,//10K, 
. or t - LaG 4 1S BT FREQUENCY O° /10K, 
LOT 2 + taG 1 tS at FREQUENCY °.F7.8,//) 
ENO 
FUNCTION DOT( HX. Y.N) 


e INNER PAOOVCT OF TWO VECTORS. 


oence aoe eoeseas 
OMMEWSION RIS}. YOU) 

OOUBLE PRECISION C 

c10.000 

DO 10 Jetlm 

YO CrCoOBLE(x( 1) JeoBrelrvis)) 

porec 

RETURN 

ENO 

SUBROUTINE DTARS[DAT,N,NP,10PT. YX, YP, RO, AJ, RVAR) 
Cocerenveatseserseenestsessenseeeeessteasteser ere Daseeaneesees 


eeeerrassesenseoese 


BURGCS ALCORITHM (1091722) OR TOEPLITZ GCRAM-SCHMIDT ALGORITHM 
FOR ESTIMATING PARTIAL AUTOCORRELATIONS Asl(t),...,As(NP) AND 
RESIOUAL VARIANCES RVAR{(1),..., RVAR{ NP) POR AR OROERS 1 THRU NP. 


nanann 


escese ences seccsreses ae neosvess 
OIMENSION DATIVE, VEUI,RETD. PCY) Masts), VAR IE 
NOP EN ONP 
DO 1 bt), NPNP 
y¥(tp 40.0 
x(t)"0.0 
+ vP11)#0.0 
oo 2 f1.W6 
vOrpepartry} 
2 m¢detdeoarttr) 
ROFCOTIOAT BAT, W)/FLOATIN) 
SIGeRo 
00 100 srt. MP 
1F(10PT €6.2) Go TO 10 
AJC SD O- 2. SOOTIN, Y.NPNPD/S(DOTIN. XR, WPNP)OCOOTIV, ¥, NPP) ) 
Go 10 20 
to ToPee-2 COOTIKNG yell vViluet} -m-s? 
Bor OOvin( ser) Ki serd M-sPeDOT( Visor) vi vett.N-s) 
As( ube TOP/BOT 
20 «CONTINUE 
avsrastu) 
SiGeSiICe( 1. -assehus) 
RvaR(J}2SIC 
00 30 Kel, NPNP 
yr(m devi Kt 
VERE EVEKPeMsueR(K) 
FO MEKRDEKIR CAD se VP (KY) 
DO 40 K22, NO NP 
KE WP NP - KO? 
40 REIRRD EMER TD 
RL1}40.0 
100 CONTINUE 
RETURN 
tno 
SUBROUTINE OTCVFIX NW REREOS MCT, ST, 


R,R,RO) 
svenessoens 


SUBROUTINE TO CALCULATE THE PIRST Met SAMPLE AUTO- 
COVARIANCES RO,MI1},. .,. RIM) OF A SAMPLE KUI1)D,..., REM) 
FROM & TIME SERIES NC.) Vid THE FAST FOURTER TRANSFORM. 


INPUT 
LP eee 2 
NFREOS : AN INTEGER (> NOM, LARGEST PRIME 
PACTOR (23) SPECIFYING THE NUMBER OF EQUALLY 
SPACEO FRECUENCIES BETWEEN © AND TWOP] AT 
WHICH THE SAMPLE SPECTRAL DENSTYTY (PER) IS 


AANnANAAKRAANnAN 


aasPie 
aasPric 
719 
arsPio 
aasPio 


arnsPeio 
aRSP1O 
ARSPIO 
ARS°10 
ARSPIOQ 
aasPeio 


aRSP10 
ARSP1O 


QaRSPIO 


ARSPIO 
aRSPIO 


ansPrlo 
° 


ARSPTO 
AaRSP16 
arnsPre 
ARSPIO 
ARSPIO 
AASPIO 
arsPio 
ARsP1O 
aRSPIO 
aRsPio 

RSPIC 
a@rsPio 
ARSPIO 
ARSPIO 
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c To BF CareurcarEed 
€ 
i4 Ourryt 
c oeaty). PRR men 1.mo alii, Rim) 
4 
¢ aun eLiaay 
¢ cr er 
€ 
c SueRoUTINGS CaiceD Oey BeCY wrt 
c 
Cerecsece estes weeneeee sveeere eeccccsaee 
c 
OlmMeNS Om AIM) Cri wPeaeos) STi nFMeOS) PERI NF REOE) AIM) 
& 
Cau. OF SP ia ww os CY. 61 PER! 
Cavs SOC vl PER mF os.m Ct st. a Mo) 
C 
aetuaa 
eo 


SvuSROuUTIng oF 


Coe mene ceencccansnseceseseneveneecessees 
r4 

c SUBROUTINE TO CALCULATE HE & LG SPECTRAL DENSITYPRR 
c Ov THE MFREOS COURLLY SPACED FRROURNCIES BETWEEN © 

c @nO YwOrl OF THe SamPce wit). RIMS. Mw ¢ WPREOS, FROM 
€ ao Time Senses & 

c 

c TWP UT 

c woth) wee) 

t WPREOS | LARGEST PRIME FACTOR ¢ 23) 

c 

c Outeut 

- PEA (I) PERI nr Res) 

c 

€ acalcetane 

c c* st 

ig 

id SUBMOUTINES CALLEO vet 

€ 

Catewrssaue eeeen weeene esansee eseenetsacen 

€ 


OMMENSTON KIN) CTC MPREOS) .STEINFREOS) PERI NFREOS) 


a 


TWOP IMSS SATAN( 1 O)*>PLORTIN) 


¢ 
OO 1 1s: mFaeos 
CTL1)00 
1 $T{b)10 
oo 2 Tet lt 
2 crtdhrewty) 
c 
CALL FFT(CT. ST. NEREOS ,NFREOS NFREOS. 1) 
€ 
00 2 1+}. mFREOS 
Do PERMIT TCTIED SCTE SD eSTCTPeSTI LD) /TWOPIN 
€ 
RETURN 
NE FCODEA(X, 1FORM, NAME) 
Covvccces secaveneceunnace ereeesnesareseane 
¢ SUBROUTINE TO CONVERT REAL VARIABLE & 
c welCm 48S @ CHARACTER F-FORMAT IFORM 
© 10 © CHARACTER ALOHOMERTC ARRAY NAME WHITH 18 
¢ im & FORMAT 
c JMPUT > WECRCH : SCRATCH TAPE NUMBER 
c Rn, TFORM 
c OUTPUT : WAME( 1}, NAMEL2) 4 CHARACTERS EACH 
Cemseses aane amen seneee teesesnoan 


COMMON /UNIT/ IUNITT WSCRACH 
DIMENSTON NAME[2). TF ORM( 2} 
REWIND NSCRCOW 
WRITEINSCACH, FORM) X 
REWIND NSCREH 
READINSCAEH, 10) NOME 
1O FORMAT{ 2404} 
RETURN 
eno 
SUBROUTINE FFT(A.B,NTOT.N.NSPAN, ISN) 
MULTIVARIATE COMPLEX FOURIER TRANSFORM, COMPUTED IN PLACE 
USING MINMEO-RADIN FAST FOURTER TRANSFORM ALGORITHM. 
@y R. C SINGLETON, STANFORD RESEARCH INSTITUTE, OCT. 19868 
BRPRAYS & BHD B OGAIGINALLY HOLD THE REAL AWD IMACINARY 
COMPONENTS GF THE CATA, BND RETURN THE REAL AND 
IMAGINARY COMPONENTS OF THE RESULTING FOURTER COEFFICIENTS. 
MUCTIVARIATE DAT@ 15 INDEXED SCCORDING TO THE FORTRAN 
ARRAY ELEMENT SUCCESSOR FUNCTION, WITMOUT LIMIT 
ON THE NUMBER OF IMPLIEO MULTIPLE SUBSCRIPTS. 
THE SUBROUTINE IS CALLED ONCE FOR EACH VARIATE. 
THE CALLS FOR A MULTIVARIATE TRANSFORM MAY BE IN ANY ORDER. 
NTOY 1S TMHE TOTAL NUMBER OF COMPLEX Of VALUES. 
“15 THE DIMENSION OF THE CURRENT VARIABLE 
MSPAN/N 1S THE SPACING OF CONSECUTIVE DATA VALUES 
WHE INDEXING THE CURRENT VARIABLE. 
THE S1CN OF ISN DETERMINES THE SIGH OF THE COMPLEX 
EXPONENTIAL, AND THE MAGNITUOE OF ISN 18 NORMALLY OWE. 

NOTE WHEN TRANSFORMING OATA TO THE FREQUENCY DOMAIN(ISHEOd}, 
frT VIELOS THE REAL AND IMAGINARY PARTS OF Z2(FREO), WHERE 
TIFRECH Sumi yer TO w) OF NI VIOEKPL Tel U- ter REO}, WHERE FREQ 
REPRESENTS ONE OF THE W FREQUENCIES ON THE GRID O TO 
CIPLC Me SsMD). 

LLSO NOTE WHEN TRANSFORMING FOURIER COEPFICIENTS BatK TO 
TME TIME DOMAIN(1SMe-1), #FY VIELOS THE REAL AND IMAGINARY 
PORTS OF Ment), FOR Jt1 TOW. 

A TAL-VARTATE TRANSFORM WITH AINI1,NZ,.N3), BI N1,N2,N3} 
1% COMPUTED BY 
Call FRTIA 


VOMZONDL NTL MTL 1) 
Cali FFTiA VOowmZeny,N2. MIOeM2 1) 
caLL FPTta MPONZONS,NI.MISMIONS, 3) 

FOR A SIMGLE-VARTATE TRANSFORM, 

MTOT ¢ © 8 MSPAN + (NUMBER OF COMPLEN DATA VALUES}, ©&.G. 
CALL FRFT(A,8,N,NLN, 1) 

THE DATA MAY ALTERNATIVELY BE STORED IN A SINGLE COMPLEX 
ARRAY &, THEN THE MAGNITUCE OF ISN CHANGED TO Two TO 
GiVE THE CORMACT IMDEXING INCREMENT ANO A{2) USED 10 
POSS THE INITIAL ADDRESS FORM THE SHOUENCE OF IMACINARY 
WALUES, £.G. 

CALL PRT(A,AL2) NTOT WL NSPAN, 2) 

AAPAVS ATIMANF), CKEMANP), STIMAXP), SKIMANE}, AND NPI MAKE) 
ARE USED FON TEMPORARY STORACE. tP THE AVAILABLE STORAGE 
1S INSUPPICIOENT. THE PROGRAM IS TERMINATED BY A STOP. 
MARF MUST BE GE. THE MAXIMUM PRIME PACTOR OF Ww. 

MAEP MUST BE CT. THE NUMBER CF PRIME FACTORS OF WN. 

1 ADDITION, 1° THE SQUAR FREE PORTION K OF WB HAR TWO OR 
MORE PRIME FACTORS. THEN MAXP MUST BR GCE. K-1. 

OIMENSTONW aft). 81) 

ARRAY STORAGE If WFAC FOR A MAXIMUM OF 11 FACTORS OF BR. 

3h MM MAS MORE THAN ONE SQUARE-PREE FACTOR, THE PRODUCT OF THE 
SOUARE-PREE FACTORS muUST BE LE. 210 

DIMENSI Im weact yt) wel 208) 
C ARBAY STORAGE FOR MAXimMuUM PRIME FACTOR OF 23 
DIMENSIOM ATE 2BF.Cmt23) BTI2F). SKIT) 
@oulrvaerence (9.199) 


AR ANRANAARAAA ARP AANANARAANRARANKRARBBONHAR AN NAANAAA 


aan 


aRSPICG 
aRBPTO 


ARSPIO 
ARSPTO 
ARSPI1O 


arsPig 
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i A tt ——— rr | 


C THE FOLLOWING TwO CONSTANTS SHOULD AGEL WITH THE ARRAY DIMENSIONS. 


manre23 i 
MAEZPS2ZO8 

If im .L7. 2) AeTURN : 
Rade CoATAN( 1.0) i 


S72"RA0/8 0 
crzscosis72) 
StI¢Sim1872) 
S1z70:s0RTIO 76) 

ImCelsy 
ieiism .ce. of GO TO 10 
872 12 
Stz08-S120 
MAD+-RaAd 
1eoe- INC 
10 WTCINCENTOT 
RSSINCONS PAN 
KSPANeKS 
WHONT- INE 
JCeRS/N 
GADFeRADe FLOAT(JC}S0 & 
100 
JF6o 

C DETERMINE THE FACTORS OF Ww 
M20 
rar 
G2 70 20 

18 Memet ARSPIO 
wracimyce ARSPIO 
KOK/16 ARSP1O 

ZO FPIK-UK/96}016 .20 OO} GO TO 18 
v3 
vues 
GO TO 30 

2S pmomet 
NPFACIMJeJ ARSPIO 
KEK Joe 

30 IF{MOOIK. JJ) .£O. ©} GO TO 26 
Je? 
ssesen? 

IF(sa .LE. K) GO TO 30 
trt(k .GT. #4) GO TO 40 
atom 

WFAC( Met }d eK 

TFtm .WE. 1) Mepet 

co TO 80 

#o 1F(K-(K/4)"4 .NE. ©} CO TO SO 
monet 
MFACIM) £2 
KeK/e 

BO KTsM 
Jr? 

60 1F(MODIK.J) .NE. ©) GO TO 70 
Memet 
NFACIM} tS 
KEK 

FO Sel luerssziezer 
IFts Le. KI) GO TO 60 

80 IF IKT .£€9. ©} GD TO 100 
sew 

90 Memes 
WRACIM) NF AC( A) 
vest ARSPIO 
1F(J .WE. ©} GO TO 90 ARSP1O 

C COMPUTE FOURTER TRANSFORM ARSPIO 

YOO SOFRADF/FLOATI(KSPAN) ARSP1IO 
COr® OFSINISO}ee? arsPio 
SO'SIN(SO*SD) ansPio 
mend ARSPIO 
teres ARSPIO 
PFUENFAC{I} .WE. 2) GO TD 400 ARSPIO ‘ 

C TRANSFORM FOR FACTOR OF 2 [INCLUDING ROTATION FACTOR) AASPTO é 
KSPANEKSPAN/2 ARSP1O 
KIaKSPANG? @RSP10 . 

210 KIEKKOKSPAN ARSPIO . 
Awca(K2)} 
OmsetK2) 
ACK2I1TEAIKK)- aK 
BIKI) Bt KK) -OK 
BURKD TALK RD OAK 
O( KK CBI KK OBK 
KKOKIOKSPAN 
TP tKK .LE. WN) GO TO 210 
KR PKR- NM 
IF(RKR .LE ye) GD TO 210 
if {kk .GT. SPAN) GO TO B00 
220 Crier O-CD 
SiSo0 
230 KIKKROKSPAN 
AKEAIKK)-AlK2) 
OrrS(Kne- (KZ) 
AUKKSATKKPOA(K2Z) 
BCRKD BI KK) oO(K2) 
MIKZ) SC VeAK-S 198K 
O(KZeSreaKel sean 
KMERZORSPAN 
IFtKK .LT. NT) GO TO 230 
KITKK-NT 
cie-c1 
KKEKE-K2 
IFtKK .GT K2) GO TO 230 
ARICY-(COsCr9sH9ss) AaRSP1O 
Sie(SoOecr-coesspest ARSP10 
cie2 {axeezestec2) ARSPIO ' 
siec 1 ARSP10 
CreCiseaK ARSPIO 
KKERKOUC ARSPIO : 
re(mm okt x2) GG TO 230 ARSPIO . 
KEERTOINCOINE ARSPIO 
KKE(KH-KSPANM)/2OIC ARSPIO 
IF({ne Le Jeeucl Go TO 220 PIO : 
co To 100 spre 

C TRANSFORM FOR FACTOR OF 3 (OPTIONAL CODE) aRSP1O 

320 KtEKKOKSPAN ansPie 

K25K1OKRSP AN ARSPIO 

Akal KK} ARSPI1O 

OKs B( KK) ansPio 
Bsa Ki poatn2) ARSP1O = 

Bs B( kK pootKz) ARsP1O 
aixcpcaneag 
S(Ke)sOKoBy @aRSPI1O 
akr-O SeAsean ARSPIO 
Qkr-0, SeBueBK ARSPTIO 
a(Kip-a( Kz) pest20 ARSPIO 

Bsr (Bier) - SC KZ} PeSIZS E ARSPIO 

aimilean-Os ARSPIO : 

emi) 
alK2} 

B(K2 OK - ad ARSPTO 

KKERZORSPAN . ARSP1O 

TRLRK LT. WHE CO TO 320 ansPig 

KR ORK NN aRSPIO 
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se 


430 


aso 


aso 


$10 


$20 


40 


Ifimne .48 SPAN) GO TD 3206 
£0 Yo 760 
SPOAM FOR FPACTOR OF & 


wFACIiT) .4E. @) GO TO se0 
KEP mMens 
RSP AN KSPAN/E 
Cros. 0 
$120 
RismaenSPan 
KIER TORSO AN 


KIenTOKSP AN 
aaPoal Kn) 
enmoeai ku) -s 
asPreainilea 
Bsmealay) 


ahr ake- ae 
Baresixn)o@(Ka) 
Oumog( aK) 
Bsrpeceiat) 
Bsmop (Kt) -BIKd) 
Oi nn} ORO oBP 
PeakP- BP 
wiiss ut. ©} GO TH 460 
ARP CAaKM- BUM 

ARM aKMOB UM 

Peekmeaum 

BKMeBKM- AIM 

VF{s. to 6.0) GO TO 460 
ALK1)SAKP SCT BKPOS) 
BK) AKPeS 1 eOKPel) 
BLK2)SBIPHOC®-BIP $2 
B(KRTPCAIF TS ZIOBIP EC? 

AKI) CAKMOCI-Bamesy 

BIRD) CAKMoe SIOBKMOCS 

WRK IORSPAN 

Rime 1€ wT} GO TO 420 
Crecr-(CO*CrosOe81) 
BSre(sOeC1-CHhesr poss 

Cis2z. O- (C®eozegier2) 
Biecwest 

Ciscreca 

CIeCress-Hiee? 

S2e2 OfCresy 

CIoczeCi-SIesgi 
SIecsaeSresze€t 

Km emK NT ONE 

tri ae L KSPaN) GO TO 420 
KK OKK KSPMNOTNCG 

trian .4& JO) GO TD 410 
Tei nsPan .€0, sO) GO TO B00 
c0 YO 100 

ARPS AKMOB IM 

AKRMEAKM- BUM 

Ones OKM-AUM 

BKM BKC ASM 

(sy we ©.0) CO TO 430 
Aled) aKe 

Biers eKe 

AiK2)1a0P 
(Kz, eBPr 
AIMS) ean 

Sika) seam 

KK OKIOKSOAN 

TF (aK i 3 wY) GO To @20 
GO "9 440 


C TRANSFORM FOR FACTOR OF & [OPTIONAL CODE) 


CICTA*S2-S7Q002 
S212. 09°C 729972 
WLUKKOKSPAW 
KICKIOKSPAN 
KIOKTIOKSPAN 
KO OKIOKSPAN 
AKPeatKIeaK 
ARMA KI) AIK 


Asmat K2) 
BuheeiK2) 
Blume B(K2)-8(K 
aten) 
{enK) 


OI KK 1  BBOBKP OB UP 
BK TAKPOCTZIOAIPSCEIOAR 


Bur snmeS727¢+Bsmes72 
a( m1) sak - Oe 
aimayran 


Orns) 
Bia peOn-Ay 

ART AKPEC2*AIP OC TZO Ad 
Ons BxPsC2eBsPasCr7-98 


Cyvsanmes2- 
MIK2 Dea - BN 


BiKd) SK 
RREKOOKSPAN 

Tr LKR tt. ww) GO TO $20 
ame ON 

1F# (me Le KSPam) CO TO $20 
co To 700 


c TRANSFORM FOR ODD FACTORS 


wewract?) 

KSPONNEKSP AN 

KSPANEKSPAN/K 

1e(K .£€0 3) 6O 10 320 
Ir-{k £9. S$) GO TO $10 
trim #0. JF} GO TO 680 
sren 

Steman/FLoar(K) 

crrcos(st) 

SdeStNiS1) 

if( se .CT MANP) GO TO 998 
CK sepst @ 

SK( UF) 20.0 

Jen 
CHCUPSCKRIKPCCoSK(KI eS) 
SKC IDCCKIKPOS)-SKIRJeCt 
memet 

Cutmpecnty) 

SRiK)e- SKU) 

Jeues 

mets 1b. KI GO TO 630 
ROKK 

RIeKKOKSPNN 

aac ai Kn 
SeeBi nn) 


16 
arerig 
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ARSPIOC 
ar ° 
ARSP IO 
ORSPIO 
ARSPICO 
ARSPIO 
AASPIO 
ARSPIO 
anseig 
19 
AASPIO 
aAsPIO 
ARSPI1O 


@AsP19 
ARSP10 
ARSPIQ 
arsPiIo 
arsPlg 
aRrsPig 
ARSPIC 
ARSPIO 
ARSP1O 
aRSPIO 
ARSPIO 
ansPlo 
aRSP10 
ARSPIQ 
arseic 
BRSPIO 
@RSPro 
aRSP1O 
ARSPIO 
aRsPlo 
ARSPIO 
ARSPIO 
aRnsPisg 
RSPI1O 
sero 
ARSP1O 
ARSPIO 
ARSP10 
ARSPIO 
ARSPIO 
AaRSPIC 
aRsP16 
RRSPIC 
ARSPIO 
aAsPio 

aARSP190 

ARSPTO 

a@rnsPiag 

ARSPIO 
aRSPIO 
aRsPlo 
ARSPIO 
QRSP1O 
ARSPIC 

aRSPIO 

aRrsPlo 

ARSPI1O 

ARSPIO 

Pt9o 

aAsPiIo 

arselo 

@RsPlo 

aRsPio 

aRsPle 

aAsP10 

aRSPIO 

aRSPTO 

ARSPIO 

AaAsPig 

aanseig 
aRSPIO 
aasetg 
ARSPIO 
aRsPlOS 
ansPro 
arsPio 
aRSP10 
a@AsPIgo 
aASPTCO 
BASEL 


660 


670 


Cae 730 


800 


810 


wel 

KieKioKSP AN 
RIOKI-KREPAN 

Jaset 
aT(SPeM( KS eMIK2) 
ARCATI VP OAK 

OTL IPeOI Kt) +BIK2) 
exsativionn 


Jeet 
AYisleaint)-a(K2) 
BtCsdep( at) -BiK2) 


KRteKtORSP An 

1F(kKt .LY. K2) GO TH 880 
1K) OAK 

O(KK) 28K 

Kis KR 

K2IOMKROKSP HN 

wet 

KRieKteRSP AN 

K2I0K2-KSPAN 

Javed 

ax 
on 
asso 0 

8290.0 

Ket 

wemed 

ameaTi Keck suleak 
OurBTiKiecKiss} oan 

Kemet 

AJTAT(K) OSK( JI} eas 
BJeMT(KeSK( su) eB 

ddr ded 

IFCvs .GT. JF} SJtss- uP 
1#(k (LT. vF) GO TO 670 
REF 

AUKt)SAK-80 

BI Ki pesKeas 

AUK2Z) FAKeBY 

SUKZ) sMBK-AV 

sesel 

wis .Lv, KI GO TO 660 
KK OKK¢KS PAWN 

TF¢KeK .2@. WN) GO TO 640 
KR OKK- NN 


IF(KK .L@. KSPAN) GO TO 640 


TIPLY BY ROTATION FACTOR (EXCEPT FOR FACTORS OF 2 AND 4) 


1F(1 .8O. mM} GO TO 800 
Keesces 

c2"#1.0-cO 

St°so 

cwec2z 

S2es1 

RKTKKOKSPAN 

akc al KK) 
A(MKPC2OaK-SIeB{ KK} 
BCRK DS ZEAKeCIEB{( KK) 
RK IKK OKS PNA 

Ir tn Le. wy) Go TO 730 
akesStes2 
S2?:SteczeCie$2 
C77C+98C2-AaK 
KATKK-NTOKSPAN 


IFIKK .LE. KSPNN) GO TO 730 


CBeCt-(eoectesoests 
SrrSte(soeci-coesi) 
CI19Z.O-(C2eezesrse2} 
Stecrest 

C2zeciet2 

KK EKK-KSPHNOIC 


16(KK .LE. KSPAN) GO TO 720 


KKEKK-KSPANOJC ENG 


IFtkK .Le. we*vc) CO TO 710 


co TO 1600 


wWPLIipeKs 

MFEKT .€9. OF GO TO 890 
MeRTOKTO 

Tete (oT. KP Kem 

vet 

wr (mers 

MPL IoD ENP (SIS F AC I) 
MPI KD ENP Kot pemrFactu) 
Jeuet 

meme) 

ets .4T. Kd CGO TO 810 
KIrMP (KOT) 

KSPANCNP (2) 

wKK2ICOt 

K2IIKSPANCT 

vel 


TF" WE. wTOT) GO TO BBO 


: C PERMUTE THE RESULTS TO NORMAL ORDER---DONE IN TWO STAGES 
f c PERMUTATION FOR SQUARE FACTORS OF WN 


C PERMUTATION FOR SINGLE-VARIATE TRANSFORM (OPTIONAL COOE) 


620 


azo 


* : eco 


AKCA(KK] 

atmm peat Ka) 
Ate2) aK 
OncB( KK) 
@(enise(K2) 
Baz) BK 

KR ORR OTC 
KIOKSPANOKZ 
Ivte2 47. KS) GO TO 820 
RIOR2-NP CD) 
Jeser 

K2INP (OT POR? 


t7(K2 .GT. MP(I}) GO TO 830 


vet 

TFdIRK LT. K2) SO TS 820 
KRAKKOINCG 

RIZSKSEPANEKZ 

Fie? .17. KS) GO TO 880 
1f(me .LT. KS} GO TO 830 
stena 

ce ro aso 


C PERMUTATION FOR MULTIVARIATE TRANSFORM 


” ese 
4 sso 


KeRKese 

AKTA(KK] 

@( RRP CA KZ) 

A(R2) AK 

Bare an) 

RqpeB( Ka} 

Baz BK 

KRSKKOTINE 

R2eKIe TW 

telex 47. K) GO TO 860 
MRORKORS- JSC 

RIERIOKRS- SC 

Teta LT. wT) GO TO 880 
KICKI-MTOKSPAR 


arseloe 
aRS?P10 


ansPrlo 
aaselo 


arsP1o 
ARSPIO 
aasPeic 
ARSPIC 
aRsPso 


aRnsPio 
ARSPIO 


ARSP10 
aRSPTO 
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RR ORK RTOUE 

Ivim2 .tT. KB} GO TO 880 
B7O RIOKZ-MP LS) 

sage 

RIeHPlLyot Pon? 

ledaz GY. BPLYl) GO Te B70 


ven 
ea0 IFiak 47. £2} CO TO BBO 
agen Kkese 


K2eRBPANEKZ 
eine .tY «8)} GO TO 880 
te(ue .47. K8) GO TO 870 
JC ena 

eso IF {7eKT¢1 .CE. mM) AETURN 
KSPHMONPEKTOS) 

C PERMUTATION FOR SOUARE-PREE FACTORS DF Wb 

seme KT 
MPact Jered 

$00 MFACIJTeMPACI J oMPac{tuel} 


vagel 
JFe(y .WE. KT) GO TO BOO 
KYaKTel 


Mm eMPACIKT)-1 
IP (WN .GT. MANP} CGO TO 998 
Juso 
PE] 
Go TO 806 

POF Uses KZ 
K2aKK 
Kenet 
KK OMPFACIK) 

P04 SaeKKesd 
Irigs .CE. x2) GO TO 902 
welded 

908 KRICNFACIKT) 
KeKtot 
KxemraciK) 
Jesel 
te LE um} GO To 808 

C DETERMINE THE PERMUTATION CYCLES OF LENGTH GREATER THAN + 

Pry.) 
co Yo ete 

ero Ken 
Keone in) 
MPU KD ERR 
te lam .WE. J) GO TO 910 
KIOKK 

Ore ured 
aKenP lt) 
Ie (ee .L7. OF GO TO B14 
FikK .WE, J) GO TO #10 
MP sdecd 
t7(s .WE. we) GO TO 914 
MARE TNCOMAKE 

C REOCROER A AND B, FOLLOWING THE PERMUTATION CYCLES 

£0 TO 960 

974 sig-t 
Ar Lars) ut ©) Go TO 924 
weer 4 

926 KSPANe ST 
1e(uy .CY. Maur) KSPANeMoxe 
JIS RSP AN 
mieMP (ds) 
KReCOKe boss 
RVERKOKSPAK 
R220 

828 KAPeRIeT 
a7iKz)salKr) 
OTinzye@inr) 
KEK T= INC 
sein mE. xe) GO TO 928 

O37 KPeKK*KSPAN 
KR2IrK1-JICOEMONPIKD) 
mKeewP (Kh) 

836 MIKE) EAIKZ) 
BierpeBina) 


wm rems- FMC 

Lo oa © e019 

Fin NE. KK} GO TO 936 
KRIK? 


TF (RK we. J) GO TO $32 
KV aKK<KSPAN 
karo 
e640 K2eK72e1 
alerprAT(R2) 
Sins psBtin2a) 
Kiem THE 
te(at .WE. KK) GO TO 940 
IFtsu .NE. 0) GO TO 926 
Iirfty .w&. 1) GO TO 924 
$So syenaet 
WTON KSPMN 
Tlewr-i NCod 
leint .C@, ©} GO TO 024 
RETURN 
c ERROR FINIGN, TNSUPFICIENT ARRAY STORAGE 
998 tSm10 
PRINT 988 
svor 
990 FORMATISENOARARAY BOUNDS EXCEEOES WITHIN SUBROUTINE FFT) 
end 
FUNCTION FIIR, MP, we) 


Creesecccsee eenene ansene seccecans 
c 
c FUNCTION TO DETERMINE THE NUMERSTOR FUNCTION IN NEWTAR 
c RADTAL FREOVENCY w 
c 
c rmpuT 
c OP MET), RNP) 
c “ 
C4 
Cossecccacareeecene pers asecoenesnestareesesesevaseeceanens 
c 

OIMENSTON Ri mP) 
c 

Ceo. 

oo 1 1et,aP 

Orercoartizd 

1 CeCoOrsMI TOS IN(OLew) 

Fiec 
c 

meTuaN 

end 

FUNCTION FIPIR, WP we) 
Cee veeee aeoae eeseeene eeense esane 
€ 
c FUNCTION TO DETERMINE CERIVATIVE OF F1 AY RADIAL 
c FREQUENCY & 
¢ 
€ reeuT 
€ orp aed... RCP) 
c « 


ARSPIO 


87 


Dimension altar) 
e 
Ceo. 
bo 1 Falla 
Olerroatir) 
1 CeCeOseoealteCosioiew} 
Fiesc 


RETURN 
ano 

SUBROUTINE FRPAREOTH NCOVM wWrREOS) 
eecses see oses 


c 

c 

c SUBROUTINE TO BUTOMATICAiLY SPECIPY & VALUE FOR © isuranios} 
(4 THAYT SATISFIES THE REQUIREMENTS FOR FFT BND THAT 18 OT LO 
c 

c 

c 


as LARGE AS NeNCOVM. 


eee ROR Perr eC ErOrererrrr rrr ivrrr rir iererrrrrrr ree ere rr reer ye 
OIMENS1ONM NROL II) 
DATA NFO/60.96, 144, 102,288.384. 480,876, 
WS TeNWeNCOVM 
1P\wST LE.1562) GO TO 10 
NFREOS SO 
RETURN 
to 1st 
20) «JF (est Le wrOli)) Co TO 30 
lejyet 
co TO 20 
30 NWFREOS TNFOIT) 
RETURN 
No 
SUBROUTINE FORIERIPF U.N. A. mA) 
Cececunecenarvassreresevessteneees esse sess eesetusanressacesnseses 
SUBROUTINE TO COMPUTE THE FOURTER TRANSFORM 
PHIIV) OF & DENSITY DEFINED OW (0,1) FORM Ve0.1, 7) 
Input 


960, 11827 


c 

4 

c 

c r.ULN VECTORS OF LENGTH MW COMTAINING FiUudiV 

c MA : MANIMUM VALUE OF v¥ FOR wHICh Putty) 18 COmeuTEeD 
ct OUTPUT A : COMPLER- VALUED VECTOR CONTAIMING THE PHI'S 

c SUBROUTINES CALLED : NONE 
¢ 
c 
c 
c 
c 
c 


NOTE: THIS ROUTINE 18 AN INVERSE FOURIER TRANSFORM DESIGNED 
TO PROOUCE TH AUTOCORAELATION FUNCTION OF A TIME 
SERIES CIVENM THE SPECTRAL OENSITY FUNCTION 


Peer er eee P Peer eee rere err ererr rrr rr rerrrr errr irri rrr rrrrrrrrr rrr ry 
OlmMeNsTOn FIN) UN) 
COMPLE A(MA),2 
TwWOP]:S SATANM{?.) 
OwerLOaTin) 
oo 20 Im:1, ma 
Fimatm-y 
AL1M) © CMPLEIO, 0.) 
O©0 10 3:1. 
TACMPLRLO., TwORPLEFIMeU{ TI?) 
10 BEIM) EAC IM) OF LDP OCEKP(Z) 
AlLIiM)peat IM) /sFLOaTIN) 
20 CONTINUE 
RETURN 
END 
FUNCTION FOFNC(X, IDON} 


Cesena wvaseeeeasnaseenasenreserasasauessesenttasensenzane 


c ROUTINE TO COMPUTE THE VARIOUS DENSITY-QCUANYTILE FUNCTIONS 

c IwePuT: 

c Mo + VALUE AT WHICH THE FUNCTION 15 TO BE COMPUTED 
c TDOH - IMOICATOR FOR THE CESIREO FUNCTION, 

c MUST @E IN TME EXCLUSIVE RANGE 1-11 

ce weosese seseeeee Bereenrecereeseeressessessace 


COMMON / CNSTNT / P1,TWOP1,RSOTPI, RPT 
DATO EPS / 1. 8-14 / 
GO TO(t.2,.3,4,8,6,7,8,9,10,119), 100N 
C COMPUTE THE WORMAL 
’ CONTINUE 
Ft kw LY ©o5) GO TO 101 
IF( Mw .GT. .999) GO TO 101 
CALL MOWRIS(H,T.1ER) 
FOFNCTEXP(-O S#TeT) # RSOTP) 
RETURN 
101 FOFNC + EPS 
RETURN 
C COMPUTE THE EXPONENTIAL 
2 CONTINUE ° 
FOFNCS1.-x 
RETURN 
C COMPUTE THE LOGISTIC 
3 CONTINUE 
FOFNC © mett. = RD 
RETURN 
C COMPUTE THE DOUBLE EXPONENTIAL 
4 COMTINUE 
Formc + 8 - aUs(x-.6) 
RETURN 
C COMPUTE THE UNIFORM RECIPROCAL 
s CONTINUE 
Forme » (1. - K)ee2 
RETURN 
C COMPUTE THE CAUCHY 
6 COMTI MUE 
FOFNC © SIN(PTeX)eo2 © RPI 
RETURN 
C COMPUTE THE EXTREME VALUE 
7 CONTINUE 
F(m .8@0. 4.) GO TH 102 
rorNc « = V.pemeocit. + x} 
RETURN 
102 CONTINUE 
FOrFNC «© EPS 
RETURN 
€C COMPUTE TRE LOG NORMAL 
s CONTINUE 
tr(e .80. t.) GO to 103 
irim .820. ©.) GO YO 103 
CMLL MOMRMIS(K,T. TEAR) 
foruc s T © {.8ste1} © RSOTPI 
ForNC * ExPl-FOFNC) 
ReTuRn 
103 FOFNC © EPS 
RETURN 
C COMPUTE THE PARETO 
’ CONTINUE 
Forme © (1.-N)oots. sBETAPI/SETAP 
QeTuRN 
C COMPUTE THE WEIBULL 
10 CONTINUE 
tein .€0. 1.) 6O TH 106 
Wormc « (5. -R)et-MALOGl) -Nbjecit. -SETAWI/SETAW 
AETUAN 
104 CONTINUE 
Forme © EPS 
RETUAN 


aasPie 
1e 
le 
te 
1e 
19 


AaARnsSP10 
ars?lo 
P19 
ARSPIO 
BASPIO 
ARSPIO 
arnsPio 
aRSPto 
ARSPIO 
ARSPIO 
@RSP10 
arsPlo 
arsPto 
AASPIO 
sP10 
sPlo 
ARSPIO 


aRsPio 
ansPig 
arsPig 


ARsPro 
arsPio 


88 


- a eR SE EES EE ee ee ee ee eee 


. € ComeurTe The MALP LOGISTIC ARSPILO 
4 CONT NUS aAsPio 
9ORNCe) -Nee? aasPlo 
¥ te RETURN aRsPie » 
ano ARSPIO 
Gg SLUOROUTING FIBAPLUL.V. RLF wm) aasPic P 
rat Cecceece nese eonses seenecegreaseeseseseos AaBPIO 89 E 
“a € aRerte * 
ie c SUBROUTINE TO PERFORM LINEAR INTERPOLATION OW ¥ YO OBTAIN F AT a 16 
ae 4 THE Mok VALUES. aRrePrioe . 
¢ aRnsPio . 
i ¢ DePuy; Us VECTOR GP VALUES AT wHICH ¥ 18 £VALUATED ARSPIO ; 
\ c Y - PUNCTION VALUES TO INTERPOLATE arselo 
F c Ro- VALU AT WHICH INTERPOLATED PFUNCTION TO BE ARSPTO 
c var eo aRsPlo 
c N+ OIMENSION OF VECTORS U AND Vv ARSPIC 
¢c mM + DIMENSION OF VECTORS x AND F ARSPTO ) 
c aRSP1O 
s c WOTK. att SOSCIBSA VECTORS MUST BE ORDERED aRSPIO : 
he rd aRsPlg 5 
i t OUTPUT F + INTERPOLATED FUNCTION VALUES aRSPiC 
. c ARSP1O0 
. € seers seeeaan esceas avesseusuneroenn euvassscce aRsPioc a 
‘ee ¢ ARBPIC . 
* DIMENSION UlwD, VIN) .RIM) FIM) aasPso aA 
. 1F(w. 20.) CO TO 100 ARSP1O ki 


gaat 
00 Go 1e1,M 
yO IFiNt)-vVtEII}}20, 80,80 
20) «65Ftt1 WE.S} GO TO SO 
FUP pov vbotegadevitpredc chy -ulsddsCutcad-ulrdd 


co To 60 
he DO FCP MEDD VPOEVETTID-VEED HPD EMET I UCT Et DISEUCTTE-UbL a ddd : 
, GO To 60 
. ao wlidbevittl 
- co To 60 


Bo pisses 
TFERE.LY.W) GO TO 10 
pie 
co To 30 
60 CONTINUE 
100 RETURN 
END 
SUBROUTINE I1CODEAIK, TFORM, NAME) 
Cece enee ase tS ae SSGe eee SESE Hates eSFESESSTESESEESEESS 
c SUBROUTINE TO CONVERT INTEGER VARIABLE K 
c wHiCm mS 8 CHARACTER 1-FORMAT IF ORM 
c JC & CHARACTER ALPHAMERIC ARRAY NAME WHICH 18 
c Iu B-FORMAT. 
4 
c 
c 
c 


abut  . WSCRCON : SCRATCH TAPE NUMBER 
K, 1PORM 
ourPut WAME( 1) .WaMEt2) > « CH CYERS EACH 
PTTETTSTTCITORI TIT TTOT LTTE TTT Tee Por 
COMMON /UNTT/ SUNTT.NSCREN 
DIMENSION MAME12), 1FORM( 2) 
REWIND NSCACH 
WRITEUNSCROH, FFORM)K 
AEWINO SCREW 
REDOINSCACH, 10) NAME 
410 PORMAT( 204) 
RETURN 
eno 
SUBADUTINE KSDIO0,U.N.0M,UM,0P UP) 


Ceasseveenaesseostennteaecaseansnesaves 


eaesecseaseure 
c SUCROUTINE TQ COMPUTE KOMOLCOROV-SMIRNOFF STATISTIC FOR 

© Tre OFVIATIONS Diu)-v. UPPER AND LOWER BOUNDS ARE GIVEN. 
c tweut o,uU,m 

c OurTeuT 
c 
c 
c 
c 


OF .uP : man [¢) OEVIATION, OP. wHIChH IS at UsUe 
Om,um : man (+) OFVIATION. DM, wHICH 1S at UrUM 
SUBROUTINGS CA 

Bu eeccereevaces 
OsmMenStow Olu) ,uin) 
SON « SORTIFLOATINI J 


om © fOlt)d) - wed) ‘. 
+ up « utd) 

om + DP 

uM 5 UP ‘ 
: pO 16 ' N 


I 2. 
ore (Ofnd) + ult) D 
1F(OY .L@. OF) GoTo 1 
up ec vit) 
op + oT ’ 


DM + OF 
10 CONTINUE 
oP s SONeDP 
om 5 SONSOM 
RETURN 
eno 
SUBROUTINE MACV(BETA,SIG.NO,R,RO} 


Coanccnn mene ees eres tetas nha taecesssseseesessescatenesee 


. c . 
¢ SUBROUTINE TO CALCULATE THE AUTOCOVARIAWCES RO,RI Tt), : 
. c -RUNO] FOR A MOVING AVERAGE PROCESS OF ORDER WO 
c with PARAMETERS BETALI],...,@ETA( NO). AND SIG (RES VAR) 
¢ - 
c 1eeuT 
. c WO. BETACI),.... BTA NO) SIC 
e 
¢ oureur : : 
c mo mity,...,MtNo) 
7 : 
Pry € SUSMOUTINES CALLED wOone 
¢ 
Ceeacecncvaeseeassetegaasesasee eenes ecetee 
ce c 
, DIMENSION BETAINO) RIND) 
c 
i cet . 
Pe po t r3.m0 
. 1 crComevattpenerait} 
> morcesic a 
oe % 
2 po 2 tet, Ne % 
Tatts 
1 ter. @0.N0} Go TO 2 
momlene-t 
Oo 3 set, nomi x 
3 CoCeSRtal sp sBETAL vet) 
- 2 webdbecesic 
x c 
ie aetuan 3 
- tao 


8 
m1 
eeaons 


ROUTIRE MAMRSLIRY RYO BETA, MR, MOTM, KI, KZ, KYOP, 
.P , WAME COP, IND. Rvan, TER) 


esos rserertsscccerectsnacoservetseasttases 


e 
€ 
€ SUGSROUTINE TO PEAFORM SUBSET MINED SCHEME ESTIMATION 
c CCCedv( Tiemtupecriy 


input 
MA RVOLAYCIN},..., RVR) : 18ST meet AUTOCOVARIANCES 
or y¥ 
mOTe : 18 MA JS WECATIVE, THE COVARIANCES CBOTAINED 
ARE WOT PRINTED 90 
MDIM : DIMENSION OF MATRIX & 1m CALLIRG aasere 
PnOcRam aaseio 
eetaiiy..... BETACMA) : COEFFICIENTS OF MRITN arnePic 
ORDER MOVING AvERA REPRESENTATION OF ¥ ARSP1O 
RVAR: STANDAADIIED RESIDUAL VARIANCE FOR MA MOOEL aasPio 
arsPe1o 


htm, Ot 


a 


mt KZ 


4 


MUCRER OF LACS OF ¥ AND E RESPECTIVELY TO Bf 
te. UDED IN THE PULL MODEL. 

Ts (K1. LT. MR-MA), THE MOST IMPORTANT LAGS OF ¥ 
ARE CHOSEN FOR INCLUSION BY THE MAGNITUDE OF 
TMEIAM BUTOCOVARIANCE. THE SAME PROCESS 18 
CARRIED OUT FOR LACE OF & IF( KZ. LT.MR-MA), 
1P1K1.20,0) THE RESULT 18 & PURE MOVING 

AVERACE STEPWISE REGRESSION. 1FI1K2.20.0), THE 
AESULT 15 A PURE AUTORESRESSIVE STEPWISE 
GRESSION. 

K1,.@0 -JI), OM (K2.E0.-32) OR BOTH, THE Lacs 
SPECIFIED BY USING THE VECTOR NAME AS INPUT 
(LAGS OF VY BEING POSITIVE AND IN THE FIRST 

Jt POSITIONS. THE LACS OF € BEING NEGATIVE 

ANO IN THE NENT J2 POSITIONS) ARE POACED INTO THE MODEL. 


KTOP © INDICATOR (1P KTOP + ©, THE BEST MODEL 
1S CMOSEN AUTOMATICALLY, IF KTOP © J, THE 
MOST IMPORTANT LACS ARE FORTED 1wTO THE mODEL) 
t SAMPLE SIZE 


ourrut 
Ke NUMBER OF PREDICTORS IN CHOSEN MODEL 
WAME(1).. -MAME(K10K2) ©. VECTOR OF LAGS IN FULL 
MODEL 


NOTE : IN VECTOR NAME, YLACS ARE POSITIVE, E LACS ARE 
NEGATIVE 


cCOF : COEFFICIENTS OF CHOSEN LACS 
IND : LAGS CHOSEN [ACCOROING TO WAME) 
TER {t MEANS ALL DIAGONALS ARE it 1.E-8 FOR 
SOME CYCLE. © 3S NORMAL RETURN) 
SUGROUTINES CALLED : SEIRES, MAK 


ANNAN ANANNANANAAANAANANANNAAANAANNANANNANANNANAAAANLF ANNAN AANA 


DIMENSION RV¥(1) ,BETAI(MA) A(NOTM.NDIM) waMEe(t),COF{1), 
1IWDI1),REV{300), REEL 300) ,.wK( 300) .NwKi 300) 

DIMENSION RYE( 300) 

DATA NOUT/6/ 


FIND SUTOCOVARI ANCES 


ann 


MOPTsMR 
109710 
sF ime GY O} GO TO 110 
1OPTet 
MRe MR 

140 CONTINUE 
ie (erm GY. 1) GO TO 103 
WRITECNOUT, 104) 

1OR FORMAT SON, mR LT.) TN Cumusi’) 
GO To e9 

103) CONTINUE 


OO 108 Itt,mR 
REVOLT) savaRmegetai!) 
RVELTI20.0 
MEE(TIFOL0 
108 CONTINUE 
MEYOIRVAR 
REEOCRVAR 
IrFttorr.#0.11 GO TO 109 
WRITE HOUT, 106) 


106 FORMATI NHI, 9 2K, IMLAG, 10K GHRYY{V) 10x, OMREV(V), AasPig a 
110K GHAVE(V) 10K, OMREE[V) 100, ARSPIQ 5 
TEMPUMLVI/S1TOX, BOT IM-)) aRnsPro a 

1140 ARSPIO . 


WRITECNOUT, 107) 13. RV¥O, REVO, REVO, REO, REYO ARSPIO 
107 PFORMAT(10K,15,2K.$6F16.7,2K)) 
PYMIREVO 
OO 108 121,mR 
PUMIPY He (REY{( 1) es 2} /REVO 
108 IFCPVAN.LE.0.9999) WRITE(NOUT. 107) T,my(1), REVEL), MVEtis 
*ROOll) Pen ARSPIO 
ARSP1O 
aRsP1o + 


DETERMINE VECTOR NAME 


nanan 


108 KK FTABS (KI) + TABSI KZ) 
IF IKK £O.0) RETURN 
00 1 terlan 
ba t vet KK 
1 alt. gpeo 
mteme 


mt 


ann 


re{xt.ct.o) Go TO 10 


we 48.0 


ane 


kis LABS(K1) 
co t0 26 
10 tFUKS €O mt) GO TO 18 


O.LT MULL Tm 


ann 


OO tt Trt ms 
wet r)eass(ay(t)) 
Tt WwKOD Peo 
OO 12 Let nt 
CALL MARI WK), WKMaAK, MIND) 
WAK( MIND) +1 
120 WKIMIND)&O. 
tte 
BDO 13 Term 
tTr(mwx(?) €9.0) GO TO 13 
WAMELKI-Ftetjel 
Ttetres 
12 CONTINUE 
ce 10 70 


ana 


arnsPle 


ann nan 


nan 


aan nannn ann 


ann 


nan 


anaann 


ana 


20 


32 


36 
36 


ao 


a2 
ai 


ay 
46 


to 


$5 


6? 


$9 


eo 


63 
@r 


1° 
oe 


OG 16 Det a 
MAME TP emt- 304 


m2: 
fe(x2. 67.0) SO TS 30 
a2.L8.0 


R201 aBs{Kz) 
co To 40 
lfix2z.260 m1) GH TS 36 


Oat m2.b7 mt 


OO 31 ber.mt 
whi i)saBsS(Revit}) 
wwkt1)2o 

OO 32 te1,K2 

CALL MAKI WK, M1, WKMAR, MIND) 
MWK EMIND Det 

watmino)eo. 

Diet 

OO 32 1et.ms 

1¥(MwKil) 0 0} GO To 33 
MAME LK TORZ-Ttelje-d 
Vvejyiet 

CONTINUE 

Go To 60 


a2 80 m1 


00 36 141 KZ 
NOME( KMS OTs K2- 9 


FORM MATRIX A 


IPL (mt we. 0}. aN. (K2.Ne.0}) CO TO BO 
1Pinz? €0.0) GO TO 45 


PURE MOVING AVERAGE 


OO at 164,K2 
Teename(t) 

oo 4a) uetyt 
Tre1aeS(1t-Namets)} 
17132. £0.06} GO TO 42 
A(t, spsRee(i2) 

GO TO 41 

all, v)sREeo 

ats, 1)sali sa} 

co Yo 60 


PURE AUTORECRESSION 


O00 46 Tet.Kt 


OO 466 velit 
P2s)aeSti)-Name( yd) 
JF(12 £0 ©} GO TO 47 
ach. slemvet2) 
GO TO 46 
ail, ,uviiavo 
als, tpeati, dd) 
co to 60 
MINEO 
DO G4 151,KK 
Dvemame(t) 
DO 84 sel, 
JI OnaMe( J} 
H2epyest 
Vaeraes(ii) 
bartaas( st) 
T3rtaes(13-14) 
If{12 GtY.0) GO TO 86 
aev: 
1f(t3.€0 0) GO TO $B 
atl, slemev(i3} 
WFC T19 91) GT.O) ALT. JPERVEC TS) 
co Tro $4 
all. sl eR€vo 
co TO 64 
tr(}1 GT ©) GO TO $8 
Ree 
1F11T3 €O 0) GO Ta 87 
ati, sie meetia) 
co To $4 
atl, Js) sREEO 
co To $4 
avy 
1f(13 €0 0) GO TO $9 
AtT, vbeMv(t3) 
co TO 84@ 
ALT. s)eAvo 


als. b) alia) 


MeRKeS 
BEM Ww) ERVO 

OO 61 Tt Ke 
Litwamett) 

1F¢i’ Gt ©) GE TO 63 
120-33 

Atl. mpemevis2) 

Ga to ’ 

att a) tie} 

atm bisati, wm) 


CALE SELREG 


INPUT TWICE SAMPLE SII FOR BREUMENT 6&6 TO PERMIT MORE 
VARIABLES TO BE SIGHIFICANT. 


C2evez 

CALL SELREG( A, NDIM KK, NAME KYOP.L2,KP,COF,1ND. IER} 
17inp 20.0) GO To ee 

PO Fo 18d Ke 

IPC emOolr) GTt.O)} CoF(tts-corti} 

conmrinue 

cCoaTtinuve 

oan MOPT 


ReTUAN 
end 
SUSROUTINE MAR( EK, &. XMAx, IND) 


ARSPIO 
50 
10 
ARSPIO 
ARSPIO 
ARSPIO 
ARSPSO 
ARSPIO 
ans 


91 


a 
easecce 


SUBROUTINE TO FIND THE MAXIMUM VALUE {NMAN) AND THE 
TNOEX GF THE Maximum VALUE (N09) OF A VECTOR » OF 
LeNCtTH wm. 


INPUT + 
es re 20 


OuTPuT 
uMAK, IND 


SUBROUTINES CALLED - NONE 


nanannnnanananaaan 


DIMENSION KIN) 


n 


ue@amer{?) 
twoet 
teiw £O.1) RETURN 


oo 1 382, 
SPURMAMLLT. MEDD) INDE? 
+ eMaweK( 16D) 


RETURN 
ENO 
SUBROUTINE MAXSPAISPEC,ALPH WEREOS NP MANIT DEL. NOFM,P IER) 


Ceeeenacsenseesasnenonasesanteuasesanaerceuasseceeseevecs 


SUBROUTINE TO OCBTAIN THE NOPM PERIOCOS. P{1},....P( NORM) 
OF THE MARIMA IN {0..26) OF THE AUTOREGRESSIVE SPECTRA 
SPECIND. .., SPEC(NFREOS} CORRESPOMOING TO THE AUTO- 
REGRESSIVE CORPFICIENTS ALPH(1),... ALPH( NP). 
Input 

WP ALPM(1), 2. ALP wie) 

WFREOS -: UMBER OF EQUALLY SPACEO FREQUENCIES 


WEEN © AND TWOP) AT WHICH SPEC 1S EVALUATED 
sPpect(t}..... SPECINFREOS) 

MONTY » MAXIMUM NUMBER OF ITERATIONS Sh NEWTAR 
1f MANIT IS NEGATIVE, NO PRINTING 18 DONE 


DEL : CONVERGENCE CRITERION IN NEWTAR 
NOTE : THE MS VAR POR THE AR PROCESS 18 INPUTTED IN 
SPECI MFREOS) 
OuTeuT 
NOFM : NUMBER OF MAKIMA FOUND (NDFM.LE.NP} 
Pit)... ., Pt NOrM) : PERIODS OF Maxima 
YER(1),. 2. TER(NOFM) : RETURN COOE FOR EACH MAXIMA 


(SEE NEwraR wARITEUP) 


SUBPROGRAMS CALLED . MEWTAR MACY FI _LFIP 


AnanannnnnnnananaAnnanaAnannann 


OIMENSTON SPEC(NFREOS) .ALPOH( NP), ST( 280), 01280), 
1A( 250) .P INP), TERI MP) 
OatTa MOUT/E/ 


c 
c 
TwOoPls8 *STAN(S.0) 
ONWIFLOAT(NFREOS ) 
HOF MAO 
NMIAMFREQS/4& 
10PT1O 
MOPTsmMaAxtT 
Fe (manit.¢t,0) Go TO 110 
MAMI T) -MAKIY 
10P Te 
110 CONTINUE 
c 
c FIND STARTING VALUES 
c 
oo t $52.1 
IPCUSPECES- 1). LT. SPEC(I)) AND. {SPEC(1) GT. sPeci(ie+t)}) CO TO 2 
co Tro 1 
2 TF{NOrM. EO NP} GO TO 98 
NOFM:NOFMe 1 
ST(NOPM)E(FLOAT(I-1)¢TWOPS)/ON 
DUNOFM) sSPEC{ I} 
1 CONTINUE . 
c 
is 
1FU(NOFM.EO 0} GO TO 100 
c 
¢c 
9s sF(rorPT €0.13 GO TO 118 
WRITE (NOUT.§$) 
BS FORMATI//, 10OR ‘ANALYSES OF AR SPECTRAL MAKIMA :°) 
WRITE(NOUT. 6) 
O FOMMATIIGR IMI, 1ON SMSTART, TOR, SHFINAL SH. SMNUMIT, 
VON THSPEC( I), 2K. IMIER/10O¥, OSL IM- 7) 
118 CONTINUE 
c 
c FIND AUTOCOVBRI ANCES 
c 
CALL MACVIALPH, 1. MPR. RO} 
c 
c FIND MAXIMA 
c 


00 10 t71,NORM 
CALL NEWTARISTIED RLF mMANMIT OF: POLI} wumMIT, reRtty) 
wetteats) we.s) Go to 
PrrPtrt 
csR0 
00 30 sti, wD 
30 CeC+2 MRI JISCOSI FLOAT sie 
OCT) SPEC (MFREOS I) /( TwoPlec} 
STUT)#TwOPIsS TINT 
PUL) eTwOoPtsPi 1) 
18 1FT10RT #0 1) GO To 10 
WRITECMOUTL IG] ELSTOLD PETE MUMET OCU, TERUT) 
16 FORMAT( ION, 18,20 ,.793.9.30, 713. 7,110, an FIs. 7,18) 
10) «CONTINUE 
00 $0 t'1,NOFM 
Catt MAK(D.NOFM, OmMax MMI6] 
STCL PoP i emtNy 
$0 Ofmminlso. 
oo $1 5r1,NOPM 
St Prides tie) 


c 
c 
100 «CONTINUE 
MARKT TeMOPT 
c 
c RETURN FOR NO MAXIMA 
c 
RETURN 
aD 


SUBROUTINE MOWRIS (P.¥. TER} 


ee 


ansPilo 
aasPio 


92 


aRSPI1O 
ARSPIO 
ARSPIO 


arsPrio 
ARSPIO 
aRSP1O 
ansPlo 
AAsPla 
ARsP19o 
aRsPlo 
aRSP10 
ARSPIO 
BRSPIO 
ARSP TO 
arsPig 
ARSPTO 
ARSPIO 
ARSPIO 
aRSPI1O 
ARSPIO 
ARSPIO 
aRSPIO 
aRSPTO 
BRSPICO 
ARSPI1O 
ARSPIO 
aansPic 
ARSPIO 
arsPio 
arsPelo 
ARSPIO 
ARSPI1O 
ARSPIO 
ARSPIO 
aansPic 
ARSPI1O 
ansPrle 
ARSPIO 
ansPio 
ARSPIO 
ARSPIO 
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MICROCOPY RESOLUTION TEST CHART 
NATIONAL BUREAU OF STANDARDS-1963-A 


Ceceevessccccrsasersen es eeeeenssnessenssseessererosenseeeresese 


€ CUANTILE PUNCTION POR WiO.1) 
oa. Coe weeneasennceace @ceecccaseseacesssssseesese 
+ Common /UMIT/ TUNTT, MSCREM 
x | vee at Y,O087), D026) C123) WIL M2, MS, me, a2 
oe vivave CVO. COUN Carers Bese, Catrs,ee 
DATA O/. F972EEEI7TES1ADEO, 1200G7HEISHAIIOS, 
-) .O41607819934072100. 00272666706437162, 
« oocee 3473023 000088655218898, 
. ©00020381812 000004327271658, 
. 00000093808 1413, 000000206734721, 
. - 00000004515 8689. 000000010416680, 
. . 088000002371 BO!t, 0000000008 28, 
. ©99000000!12 .0090000000078138, 
. - 0000000000087 0000000000018 81, 
, . ©00000000000374, cococcoeeooj‘ss, 
. 00000000000002!, c90000000000005, 
.- ° 000000000000001/ 
. BOTA C4 HIPS MSOTA ISHS, - 0162726072818 E7EE4, 
x ° 000433866472 0002146388 70074, 
« 9000024257510 , 600003072168 1080, 
° + ©00000012466062. 000000062406808, 
. ° +. C000000OOE4O172B5.- 000000001423208, 
. ©000000000743 ©00000000033585, 
. + 00900000000 + ©090000000000810, 
. ©00000000000083, 000000000000020, 
. > €00000000000002/ 
pata o/. 79709020083. -.623107004309065, 
. + 00637 6097508.-. 0008 226681. 
. > ©00010961022307, 000028108847025, 
. 0000 10862336088. 0000027546412330, 
° 0000000324684 + 9000000206306337, 
: ‘ + s090000083 + 900000017684010, 
. + 0000000039012 000000000186932, 
' ° ©0000000072779373. 000000000132817, 
. ©00000000031834, 000000000001870, 
. + ,00000000000203 ©COKCCO0OOOCOCDSS. 
: . = ©00000000000220.- o00000000000010, 
. ©00000000000013, 900000000000008, 
: . ,000000000000001/ 
OATA MI MZ, HS MERZ -1 61304623733, 
. 2.8 O1rz231er7e $6945763132983, 
e 2. 2E7HIS THO IETS, 5 414213 IITII/ 
OaTaA Mime / 4,810 / 
Re t oPoP 
tea 2 © 
SlLGmMa +s StSH Cs. KS 
tel woY.(k.GY -1 ano x i7,.3.3) GO TO 48 
aOsinx) 
eT 8) GO To 278 
zezs 32-1. 
22 
oy 
1 
et 
eo, 
a] 
s alieer) 
's © 6 © B11PhP OL 82) © He 
* u@ © we 2 -HI 
«ome 
oe oY 
© c@ze- 4 
YFtLez Le ~) GO TO 1$ 
CO TO (20,36).4 
20 vy 4 J # 86 * SIGMA 
GO TO #0 
29 B © SOMTL-acac(1 -292)) 
Lege Gt -99741 69 186 36 
woe MTee oN? 
1pm 1 26 
tos 2 
noe 46 
co TO 10 
3o we mee 8B + HE 
19h 6 ot 
woe 2e 
toe 2 
co To 10 
35 vos 86 8 KG * SIGMA 
ao ve Atev 
RETURG 
as vos SIGMAGKING 
fem « 129 
WRITELIUNTT. told 
YOU FORMATIS////////1R,. TSMBCCOTERROR CONDITION IN MONRIS#t0ee 
metuen 
ewo 
SUBROUTINE MEN( XN. RMIN. TWD) 
Co seme re oesresnenenerssovreeneennaseneseseesensseeseser 
c 
¢ SUBROUTINE TO FIND TRE MINTMUM VALUE [NMIN) OND THE 
c INDEX OF THE MINIMUM VALUE (IND) OF A VECTOR x OF 
c LENGTM M 
c 
c berpur 
c Womtty. 2 OM) 
c 
c outeut 
c uMIMN, IND 
c 
t SUBROUTINES CALLEO NONE 
c 
Ceonee venees senese eeees eeares sense 
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t eP 
t 
a] c 
c 5 
fot c 
be. c 
be c 
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c 
c 
€ 
t 
: c 


OIMENSIONM KIN) 


uminmaxty) 
Iwoe1 
tw €0.1) RETURW 


oo t te7.m 
TO Cami GT eC ehh Ueaed 
meinen (ino) 


RETURN 
exc 


SUBROUTINE MESPSL(COF, IND STE. KP _NFREOS CT. ST, SPEC) 


eeoetenenseeesasesnaeseesseeeeeseneee 


imPuT 


RP : NUMBER OF VARIABLES I” MODEL 
raoti), -dMOL Ke) 
{¥ LAGS ame POSTTi VE © LAGS NEGATIVE) 
corti. COP URRY COEFFICIENTS OF CHOSEN 
vant aces 

StC . RESTOUAL VARIANCE OF CHOSEN MODEL 
wrneos NUMBER OF FREQUENCIES OE TWEEN 


© amd TwOorl AT WHITH SPECTRA 18 TO BR CaiCcuiaTted 


ROUTINE TO CALCULATE SPECTRA FORM SELECT MIXED SCHEMES. 


TMOICeS OF CHOSEN VARIABLES 


4th 


aasPio 
AASP10 
P10 
ARSPIO 
aasvioc 


aRSPIO 
aRSPILO 
arsPlo 
aRsPlo 
ARSPIO 
@asPig 
aRSPio 
RSPIO 
RSPIO 
aRSPIO 


aRSPIO 
ARSPIO 
ARSPIO 
aRnsP1o 
ARSPIO 
ARSP1O 
arsPio 
sro 
sPi9o 
aRSPIO 
BRSPIO 

asP1o 
ARSPIO 
arsPero 


GRsPig 
argsPlgo 
aRSPI10 


aRsPio 


93 


a - C TL 


ec Cureur : 
e SPECI)... . SPECI NFREOS) a 
c a 
c SuBAOUTINES CALLED : FFT a 
c a 
c sees eeen eses eseee sere eee 
€ r 94 
OIMENSION COP(KH). TWOLKP) CTI mPAEOs) STI NFREOS). a 
18PEC(uraeos) a 
€ a 
c 
vac C/le.catants o}} 
1PIK® GY 0} GB TO 20 
DO tO ter, wFRROS 
1o | «sPpeclijerac 
co To 98 
c 
¢ 
20 06 30 1513, HFRROS 
cTtibeo 
STIT)e0 
3O «=SPRECII) FO 
c 
CTUtpar. 
svcsdes. 
00 BO 121,KP 
Vietwolr} 
WI11.L7.01 GO TO 40 arsPio 
svTerievrecorii) arsPie 
co TO 80 


ao) odle-dt 
Crt bles) Corts) 
$0 CONTINUE 


c 
c 
CALL FRFTICT, SPEC NFREOS .WFREOS WFREOS, 1) 
c 
DO 6O Ltt, NFRECS 
SPECITPOCTCIPCCTE SD oSPECIT eosPEc(1) 
to crtriso. 
c 
¢c 
CALL FFT(ST. CT. NFREOS MFREOS WFREOS, 1) 
c 
DO 70 §51,MFREOS 
70 SPECIIDOLPACHSPRCCIDISISTISPOST(RIeCTIS Dect Es)) aAnePiO 
c ARSP1O 
c aASP10 
99 CONTINUE ARSPIO 
RETURN ARSPIO 
END aRnsPie 
SUBROUTINE MEWTARIST RNP MAKST,OEL, 2.NUMIT, LER) ansPio 
Cesevesnsarerevevesev assesses ranesosesesessessisesunoe AASPIO 
c ARSPIO 
c SUBROUTINE TO OFTERMINE THE RADIAL FREQUERCY 2 
c OF THE AUTOREGCRESS IVE SPECTRA CORRESPONDING TO 
c RIT}... , RUMP) CARL.) OBTAINED BY CONVOLUTION OF 
c acveuly), . ALPHI(NP)) FOR STARTING RADIAL FAEOQUENCY 
c VALUE ST. 5 
c 
ce IWPUT 
c WP ST RIT)... RIP) 
.c¢ mMaxit MAXIMUM NUMBER OF ITERATIONS 18 
c MEWTONS ALCORITHM 
c O8L : CONVERGENCE CRITERION 
c 
c ourrpur : ARSPTO 
c 2 ARSPIO 
c NUMIT ; NUMBER OF ITERATIONS WEEDED (16 arsPio 
Ld NONCONVECEMCE NUMITEMAXITO1} ARSP10 
c TER ' NORMAL RETURN arnsPio 
c 2 ZERO FREQUENCY ENCOUNTERED ARSPIO 
c 3 ZERO DENOMINATOR ENCOUNTERED aARSPIO 
c a MONCONVERGENCE INDICATOR aRSP19O 
c arsPloe 
c SUBPROGRAMS CALLED : F1,FIP ARSPIO 
c ARSPIO 
CHeesesseeaeesessaansesresesssaseeseseseteseesetanaveseneese aRSPIO 
c AaRSPIQ 
OTMENSIOWN RI NP) aRBPIC 
DATA NOUT/6/ 
ce 
c INITIALIZE 
c 
1EAet 
vet 
wOesT q 
© q 
c CEGCIN ITERATION 
c 
1 Iftwo.tt.1.8-18) GO TO @ 
Carita, mP, wo) 
CPF IP (RNP, wo) 
iF (aest(cr).17.1.€-28) Go to & ARSPIO 
whewO-16/GP) aRsPlo 
CASS (wn-wO) /aes (wo) 
wtic.4.7.081) Go to 2 
1717 €O.mMaxIT} GO TG 3 
ves tes ARSPIOQ 
wOrwn ARSPIO 
co Ta 1 QRSPI19 
c 
€ CONVERGENCE RETURN 
c arsPia 
2 wuUMETalt ARSPIQ 
town aRSPIO 
AETURN 
c 
¢ WONCONVERGENCE RETURN : aRsP1Q 
c ARSPIO 
3 MUMITIEMANITOS ARSPIO 
1ERra AaRSPIO 
2awn 
RETURN 
i3 
€ ERROR ARETURK 
4 
a WRITE(NOUT, 10} . 
10 PORMATIION, “ZERO FREQUENCY IN NEWTAR’} 
TERS? 
RETURH 
c 
8 WRITE( NOUT i] 
3S FORMATLIOR, ‘3ERD CENOMINATOR IN NRWTAR’) 
AETURR 
i4 
Eno 
SUGROUTINE F 


i 4 
t 
c SUBROUTINE TO CETEAMING AVUTORECARESSIVE COoRFFICIENTS 
c ALOMALID, ... ALPMAENP) SIVEM PAATIAL AUTOCORRELATIONS 


PARTI Ve, .. PARTI MPD 


impure: 
mP PARTI ND, ..., PARTI MP} 


GutPuT : 
S.Pmatid,... ALP MaAleP) 


anePio 

SuGROUTINES CALLAO : NOME aAsPio 
aasric 

cesar as eenseccsee sare soneaseeseresanesseeesesaensencensseen aasPio 


Nnaannnaannanan 


OIMENS ION PART(NP) aL Pnai ne) aAsPio 
ARSPio 
wPet aAasPIOC 


nnn 


ALPMat Ti sPaRT(t} 
IFinP BO. t3160 TH 98 


we GT 


FLk ORDER IMP tw RECURSION 
00 20 INPe2, uP 


VENO ALPmMAl OD, ... aLPMALING-1) 


ann annanna 


ALOHAL INE} © 

ClrALPMAL INP) 
Le tMPs/2 
oo 
INP MeN 
TEMPO SLOMALSIOCHALPMAL ING MS) 
ALPRALINPMIDEALPMALINP MID OCOAL PHA LU) 
avpmaltsls TEMP 

COMTIMUE 


CONTINUE 
RETUAN 

ano 
SUBROUTINE 


econ es sess . 
SUBROUTINE TO FIND BIASED RESIOUAL VARIANCES RVARII1),..., 
RVARIM), CIVENM PARTIAL AUTOCORMELATIONS PART(1},.,.,PART(M) 
and VARIANCE RO. 


bwePuTt 
Mm, PARTIT},..., PARTIM) RO 


OuTPurT 
AWVAR(T),..., RVAR(M) 


SVUGROUTINES CALLED : NONE 


nanananannannaann 


DIMENSION PART(M} .AVARIM) 


ta) 


RVBAL 1} eMOF(F.-PART{ 1) 202) 
1FEM 80.1) RETURN 

oO y Ts2,m 

MVARUITD ERVAR(T-t)50 1. -PART( 1) e682} 


RETURN 
eno 

SUBAOUTINE PARZTIAVAR, M,N. CAT, WORD) 

oer eas aneesseane ARSPIO 

ARSPIO 

SUBROUTINE TO DTEAMINE THE CROER OF AN AUTOREGRESS IVE ARSPIO 

PROCESS BY PARZENS CAT CRITERIA arsero 

ARSP10 

: aRsPic 

mM, MVARIT),..., AVAR{M) ; STANDAROIZEO RES VAR ARSPIC 

FOR ORDERS 1 THRU mM. ARSPIO 

MN : SAMPLE SIZE ARSPIO 

ARSPIO 

oureur : ARSPIO 

MORO : CETERMINED OROER ARSPIO 

CaT(iy,...,carim) arseto 

aRsPro 

SUGROUTINGES CALLED : MIN ARSPIO 

ARSPIO 


ANAANNAANAANKRAANnKAANN 


OIMANSTON RVAR(M) ,CAT(M) 


n 


OneFLOAT(N) 
dO 1 Fe1,m 
Cro. 
pO 2 vest 
CeCol et. + (PLOATIII/ON) S/AVAR( J)? 
Csc/om 
CAT(SPSC- C09. -CPLOATE TI SON) D/RVARiT)) 
CALL MINICAT,M,CATME ,NORO) 
FC CATME.GY.-1.) WORDEO 


RETURN 
eao 
SUBROUTINE PEGRRAIR,AO.M, CORR) 
aes eees 


SUBROUTINE TO TRANSFORM AUTOCOVARIANCES RO.AL1),...,R1M) 
TO SUTOCORMELATIONS CORRI1),...,CORRIM} AND PLOT THE 
VECTOR (1. , COMMIND,..., COMMiM-1)} VIA CLPLTS 


TmPUT : anrsPro 
M.RO,AL TD... RIM) ARSPIO 

ARSPILO 

euteur : aAAsPlo 
Comat), ..., COMRIM) aRSP106 

aRsP10 

SUOROUTINGS CALLED . CLPLT1 MAK, min, @RSPIO 
aRrsPero 

coceseeeneseeeseeseeeeseseseeeseneseateseae 


AAnnAnaAnnannannann 


Ormaenstom Rim), CORRIM) ARSPIO 
OaTa wOuT/6/ ARSPIO 
ARSPIO 

waite(wour,2)} ao 

PORmMAT(IOR, PHRIO) «© ,816.7) 

sete ey.) To 10 

Comnitsemi tp sro 

warreimout,&) connii) arnsPic 

PORMATIION, “COMR(1) © *, 818.7) ARSP10 
Tuan ARSPIO 

CONTINUE ARSPIO 

OO t tetle ARSPI1O 

COMALS emi} sae ansvieg 

CALL CLOLTOI CORAM, 1, 4NCORA,1.0,-1.0) ARSPIO 


Le ket a ee OL Te ee ee ae ee eS gee ee Pe 


®ETuUan 
anup 
SUBROUTINE PLIMEL MR. WE ML START. BND. CMARL , PAaRRaY ,.YMIn, VMAK, VINE) 
DIMENSION STARTINL) ,ONO( mL) Canc (mr), 
wmemi os we ansPrie 
pwc 8 1. / PLOATINGMI} anerie 
“ 1e« wR oe 8 aage 
é ; 96 | 

00 10 the 1, mL ansrie 
if ¢ STATE ILE - @NOtILI } 20.30.40 aansria 

20 BIRTCTLD © amautiSrar( iid, vere) 
@nupise) ® AMIN ANDI Te) .vman) 
cero BO 

ao STAFCIL) © AaMIMIIBSTATIIL) vman) 
eworin)s «© AmMABIIENDE IL) .¥MiM) 

bo BS © FEWOLILD - STRTCILED © PHC © vine 
S 9 CHTATITRY - veer devinge t & 


two ¢ ENTUS} 

bo $0 3 8 2, mem 

s Ss ¢ OS 

ie ( IMTtS) .€0. IMD } GOTO 60 

ted © TMTIS) 

VUNRP TINO, 1) © CHARLIE) 

eo ComTinvue 
core 160 

Jo im ¢ S¥MTETLD 8 
imo oe (STATUTE? 
00 70 t © 2,NCMI 
PARRAYIWAPT-TND,1) © CWARLI TL) 

10 COWTINUE 

10 CONTINUE 
RETURN 
eno 
SUBROUTINE PPLOTIN V.NY. NCHAR,KZ,Z,N2,12, 2OMAR, 10RO, VMN, vMR, 

a MLIN® STAT O,CMARL LX LY¥,LC LE) 

Ceececeoaneesaeseeses aevtese eens sesesesasseeesas 

ROUTINE TO DISPLAY PRINTER QUANTILE-BOK PLOT 

repur 


VMIN .OR. SYRT(10) .Ge. veant GoTo 10 
vmMtw) © VINE © 1.6 


HR + VECTOR CONTAINING THE VALUES J/iNvet) WHERE 
9t1,2,3, ..,.4¥ 
Y + VECTOR OF SIZE NY YO BE PLOTTED 
MCHAR - CHARACTER FOR VY IW PLOT 
Lm byoe ¥ TORS OF SITE 2 CONTAINING THE LABELS 
FOR HX ANO VY RESPECTIVELY 
LC - VECTOR OF SIZE FO CONTAINING THE CAPTION FOR 
THE PLOT 
2 - OPTIONAL @cTor OF SIZE M2 TO Be PLOTTCO 
RZ + ABSCISSA FOR 2 
NLINE STAT, END - VALUES FOR SUB. PLINWE WHEN 2 18 & LINE 
VMN, VMK - MIN AND MAX VALUES FOR ORDINATE OF PLOT 
190K - £OvAL ZERO TF BOX PLOTS ARE NOT waNTED 
12 - €Ovart 1 1F VECTOR 2 18 TO BE PLOTTEO 
JCHAR - PLOTTING CHARACTER TO BE USED FOR THE 
ZI VECTOR. MUST BSE CIFPFERENT FROM 0. 
Li - OP TIOWAL LABEL TO FOLLOW CAPTION 
Joro « 1 1f Y 18 GRADERED FROM MIN TO max 
a) IF v 18 WOT OROERED 
* 2 FOR WORTTONTAL ZERO-LINE (FOR 10 PLOT) 
aencee eeors 
COMMON /UNIT/S IUNIT,NSCREM 
DIMENSION PARMAVIS! 61), CHARI(10) _XINV) Vimy 


nNANnAANNAARANAANANANAAANAN 


e LEOM2O) RIEU D. TIMBILLTC ZO). MRRCNII, START] 
DATA BLANK, DASH us .CcaPl : ARSPIO 
DATA CHBR1,CMAR/'O','°O*,'°2 ARSPTO 
OaTa €°$ / 1.€-10 / aasPrio 
vine « & aasPero 
mIwc # 6 aRrsPlo 
wRoe 1 ¢ 10 © SINC ansPio 
wt 8 1 6 10 © KIWE ansPic 
wePpioe WROO 8 aasPio 
WARS © WR Tt aRnsPlo 
weet © we; ARSPIO 
Fc os FLOATINGMI} ARSPIO 
WUNTT © FUNTT arsPro 
CwAR{( 1) © KCHAR sP10 
Te { CHAR(1) .€0. BLAaWK ) CHAR(1) © Chany ARSPIO 
ree1ORO-1) 2,3,.4 aRSPio 
2 CALL MINTY, MY, VMIN, TER) aasPro 
CALL MAK( VY. WY, Vmax, Tea) sP10 
witz .€0. © } Ga Yo 1 aRSPIo 
COLE MIME ZT, WZ, IMIN, TER) arsPrio 
CALL MAKIZ, NZ, 2MAK, PERI aasPic 
TEC ZMIN LT. VMIM) VME st 2M sPro 
TEC IMAN. GT. VMAX) VMAX? ZMAX arsero 
coro 1 aasPrio 

3 vein oe vit) 


vman os YONY) 
wrw(rz £oO. ©) Goro 1 
zeim oe Bit} 
wax» FtNw2) 
IPCIMIM LUT, VMK) VMIN © MIN 
TFC 2Max .CT. VMAN) VMAN © IMAX aRSPlo 
coro t ARSPIOQ 

4 vMIN © veRN 
veaR « ven 

1 comTinue 
AONGE «© VMAK - WMIN 
tf ¢ RANGE ,LT. EPS ) GOTO 88 
VINC ¢ FLOATEWAMS) / RANGE 
90 10 te 2. MRM 
PARRAY(],1)2C API 
ParmRav(1.NC} © CAP! 
oo Jt2, NOME 
PARAAY(T, J) SBLANK 

s CONTINUE 

10 CONTINUE 
bo 20 set. Ne 
Parmav(i.s) * DASH 
PARRAYINA, J} eDASH 

20 CONTINUE 

16 Jeo 
s*0.0 
De 28 354,NC, KING 
seuet 
xJvisbe 8 
s\se ARsPio 
PARAAY(MR, 1) 2PLUS ARSPIO 

28 Continue aRsPro 
be 30 1 fe 1,08, VINE arsPre 
PARRAV(MRPT-1,1) © PLUS aRrsPig 

30 ConTinug aasPto 
Je ( MEINE .GE. 1) CALA PLINEI(MR HO WL INE STRY, END, CHAAL, 

a PARRAY, VMIN, VMAK. VINC) 

to CenTINUE 
oo Reiwy 
i] txpere ¢ 9.8 
TMOVELVIKD-VMIMpevINE © 9.8 
tre? we. 1360 76 70 
IMORE(CT(KD- VIM PeVINe © 1 
PARRAY(WRPT-I1NOV, INOX)2CN 
PARRATIWRPT-1NOZ, INBR) SICH 
IPL I1NOR. BO. SHOV)PARRAY(MRPI-INOY, NDR) CHARI 10) aRnsele 


so te 46 

CONTINUE 

1P i Cuamitd .60. NCwaR ) COTO 44 
Bo 42 109.10 

se ipaanar{ -tmOv,1NOK)}. €0.CHaRt1)1)G0 TO 63 
Continue 

PaRRav( MAP L- INDY, IwONPecmaRt i} 

co TO 48 

1711.80 sojpdet-s 
PRORAVINRPT-IMOV, EWON} SCHARI1O1) 
Continue 

iF 


es vime eg 
.Impn2) #0. Cmhanti} } GoTo 68 

PARRAVINAPT-INOZ, FHDEZ) © OMAR 

GOTO 47 

PARAAY( NAP Y-1NOR, TWOKZ) © CHARI10) 

continue 

wWhITEIWONIT. 1OO1LC 11 ansPrio 
Ove «© EPLOATIJSING) / FLOATIWAMI1}} © RANGE a 1° 
YHALLE © vmAr ¢ OVE 

OO 3200 Let MR, SINC aR3P10 
THUMP ARMAYIWRPIT- 1,8) WE. PLUB) PARRAYINAPH-3, 9) s2ORO ARSPIO 
CONTINUE aRSPIO 
oo 6 i aRsPtO 

: LUS). OA. (PaRRAv(1,1).€0.2E€0}) GO TO SO 

WRITE NUNIT. FOR) (PARRAYI1, JI, J0t, wo) aaseio 
Go TO 68 ARSPIO 
VRALUE ® YEALUE > DY arsPio 
WRITETNUNIT, THODYRMLUE, (PARAAY(T, J}. vet No) ARSPTO 
CONTINGE 

WRITEINUMIT, LIOIVMIN, [PARRAYINA, J} ,491, NC) 

WRITEI NUNIT VES EKG 

WASTECMUNIT, T2OILA LLY 

co TO $98 

CONTINUE 

WRITEDNUNTT, 199FKC 

WRITRINUNTT, 200) 

FORMATi' A VALUES’) 

WRITECMUNST, JOSDEME TI Det, wy) 

FoanatT( tx, 10F10.8) 

WRITECMUNIY, 210) 

vormati(*’ v values’) 

WRITECNUMTT, 2O0S}(4VE1),T89 Nv) 

FORMATI’ ERROR IN PRPLOT POR THE PLOT OF *, 2086) 

RETURN 

PORMATI OMI, 15K, 2000,/,16K,2006,/) 

FORMAT(I3N, 6141) 

PORMAT(IN, F10.3,2K,61A1) 

FOUMAT( N24. 399F3 1.303) 

FORMAT( IER, “ABSCISSA 15 °,284,4%,° , ORDINATE 18 °, 208) 

eno 

SUBROUTINE PRPLOTC(W.K.V,2.M, 1CAPT, NAMW,NAMK,NAMY .NAMZ, MM, SMIN, 
*Sman} 

CocsaetseevesesecapsesseaneresessstseearaseaesesesSersersesevensesses 


SUSAOUTINE TO PLOT LOCAL OUANTILES OF PERIODOGRAM FROM ansPio 
SUOROUTINE SLOCAL. OCUANTILES ARE PRINTED, LOC(OUANTILES} @rnsPio 
are PLOTTED ARSPI1O 
ARSPIO 

INPUT Mw, KY, 2 - WwW TS OROERED ON INPUT AND KRiTisxiwits), ETC. ARSPIO 
ICAPT - LITERAL CONSTANT FOR CAPTION OF PLOT aRsP10 
NAMw. NAMK .NAMY, NAMI - © CHAR, LABELS FOR wx, ¥,2 ARSP10 

mm - wIDTH OF PLOT (¢#101) ARSPIO 
SMIN.SMAX - MIN OF (X,¥,2) AND MAX OF (K,¥.2) BRSPIO 
ARSP10 

SUBROUTINES CALLED: FYERP, MAK, MIN arsPio 


savoeeuesasecersen weesnesseese eeasv 


c 
c 
¢ 
c 
c 
i 
c 
c 
c 
c 
c 
c 
c 
e 
c 


COMMON /UNIT/ IUNTT.NSCREN 

OLMENSION wiN) KEM}. YONI. ZINN), 1CaPt (20) 

OUMENSTON w114G), R1146)} VEG), 21066) actor) 

OATH BLANK, DOT, STAR. SL. PLUS. SM/IM , 16. ,1M0, OMT THO, THM 
1OPTK#O 

te (m.GT.23) CO TO 20 

WATTELTUNIT, 10) N 

FORMATIION, ‘SAMPLE SIZB OF °.12,°' 18 TOO SMALL IN PPLOTG,‘) 
co 70 999 

COMTI WUE 

AMINOALOGISMIN) 

AMANTALOG(SMAR) 


CREATE wi VECTOR OF E£OUALLY SPACED w AND INTERPOLATE TO OBTAIN arsPio 
co SPONDING K,Y¥, NO 2 VALUES sP10 


Oecrlwimi-mwi1))/65.0 a@aasPle 
00 20 121.406 aRSP10 
WHT bewC uber Gast s-tssoNec ARSPIO 
COLe FTERP (Ww, RK wt, RIL ARSPIO 
CALL PT tw, Vimr,vi, ARSPIO 
CALL FYEAP iw, 2,WwI, 21,8 ARSPIO 
ARSP10 

IMTTEBLTZIE AL AASPIO 


Ons (mm 19/2 

DO 40 sit, mM 

avisispor 

WAITE(LUMIT, 48) 1CaPt 

PormaTi/s/.410,20a6,/} 

WRITE(TUNIT, SO} NAMW, WAMK MAMY ,WAMZ, (AL(J),J21,MM} 
PORMATI/S/. 16K, 86,60, 84,60, A4, 68, AG/ION, 4O19M-), 2%, 101A8) 
oe 60 vt1. me 

aclu) Be ann 

But iysse 

Avtmapase 


Rvs. 2elamax-AMia) 
TOERV LT. 4. +20) TOP THet 
ansPie 
or e190 
ete 
06 100 vs1,66 aasPic 
1P{ 100TH £0.15 6o TO 706 arnselge 
Chez efFALOGIAtis) )-amimrsay - 1. aaseio 
C22. e(MLOGIVtls))-samiadsay , 
C32. LAL OST Its)? -amMINDAY 1 
co Te ao 
C120. 
C2eo. 
C3ro. 
Kweomeicres. 992.8 
mv eOMmeICze1. ped 6 
K2soMe(C3e4. 942.68 
ac(an}eoor 
ALIRT)©STaR 
SLiKT) 007 
SPURS. RO RY) ALIRND Se 
IFC mR. BO KEP SLIDES 
IPURY. 20. M2} ALIRVDsBM 


ao ge ge ag ee gg ee fae ee ae 


WRITE LL UMIT. 
9@ FoRnmatTi tos 
arianbe 
OUCKV) OL ame 
OLCa2) 8 am 
100) «6—CoNtinue 
CO 140 101, om 
180 «aclibeoor 
@.tiberPius 
Ou lem) ee.us 
WATTEC FumiT 120) Paci). bet, mop 
1200 «FORmAT( ION ao(in-) 2m 10tat) 
PMAR CAV OCAMEN 
VEL LUMI. 1301 amin oman 
120 FORMMATIG7R #10 6.608 710 4) 
998 CONTINUE 
RETURN 
two 
SUBROUTINE PREPSPO (SPEC WN REOS. LOPT OLY. WK. SPMIN, 


Cocenccerscccaressscesesesenessvassocenessessrasesanes 


OL mils. WHLGDL Vdd ad BEC GT CALCUL, tet, my 
© 6.32. FS &,tK,. 79 &, 18, -O.2n 18a} 


c 
4 OuTiIne TO OiviIde SPEC Br Civ. Tame 186. AND CLIP 
¢ TO -6.6 BLSO MAX AND MIN OF RESULT 1S FOUND. 
if 
€ TmPUT 
c weReos srecti). > Cimpngos) Orv 
c sory 1 MEANS AEOS ARE © TO 2°71 
c © MEANS #08 aff © TO PI 
i4 
€ OuTPur 
€ wetid. WK (Mt) WHERE MI TH IMPREOS/2Z5¢1 OA NFREOS 
¢ Semin Seman 
i4 
Ceoauenravesesenssacenseveesene eeaanacesennsasece eee 
c 
OMMENSION SPECI MFREOS) KI NFREOS) 
c 
c UML TIALIZE 
c 
SPmimes. 
SPmans-6 
Bora. OCIidIV) 
CmMineDIVeERP(-6 0) 
em Dlveagrete oO} 
w NFREOS/2)¢1 
1F(10P7 BO Of} MICNFREOS 
c 
c Tam€ 1.OC and CLIP 
id 
00 JO fet wt 
IF1SPECTID GT man) GO TO 10 
we t(srecti) Lt Emind CGO TE 20 
waitdeaioc (t)b)-00 
IPiwKIT) &T SPMIN) SPMIMewKi I) 
TF iwett) CT. SPman} SPOmMARewK (1) 
co ro 30 
10 CONTINUE 
we (tsb. 
SPmManes 
co T0 30 
20 CONTINUE 
wei rir-6 
SPMInNe-6 
30 CONTINUE 
c 
c 
RETURN 
eno 
FUNCTION OF TND(O.N. OUANT) 
Coeccservevevascsoncuvessessenateceassanesneasteevenasenenses 
c eureose TO FING THE PERCENTILE VALUE OF O 
c ay OUANT. 
ec HP UT: 
¢ QO - VECTOR OF SIZE Ww 
c Ww + NUMBER OF VALUES IN ¥ 
c quan? VUANTILE VALUE 
Ceeeaese e eeeee seeste sereenee 


COMMON /PARMI/ KMIN 
DIMENSION O(N) 
RP « FLOATI Met) © QUANT 
tot TMTER) 
Pt wR 
aF{ 1) we. 0) Soto + 
QFIMO © (4 -P} © RMIM © Be Ott} 
RETURN 

’ tect our. wd) GOTO 2 
Orrmo =: Olm) 
RETURN 

2 OFiImD 2 {t.-P) © OF1) © MP & OLEet) 
RETURN 
eno 
SUBROUTINE OLOCALIB,UL& 


SUCROUTING TO COMPUTE AND DISPLAY LOCAL QUANTILES 
lwPur - OBTA SET FOR WHICH LOCAL QUANTILES ARE OFSIRED 

+ ABSCISSA FOR OROINATE &Z (Rh 18 A FUNCTION OF U) 

10W OF VECTORS XK AND UW 

oF LOCAL CATA TCMES TAKEN FROM DATA BATCH K 
WOME - MEADING FOR PLOT 
TPLTV - 1 TF VERSETEC PLOTS DESIATD, O BIHEAWISE 


ourrurT PRINTER OND VERSETEC PLOTS OF LOCAL QUANTILES €28,680, 
+ GPUID. ©,78, {8 THE P-TH PERCENTILE 
vaLUE MW The DATA SATCM KUT), RETOID, ©... METOReN), 
Te@ ABSCISSA VALUR UGII) 18 YAKEW TO BE VITeK/2) 

1? & 18 OOD AND THE AVERAGE OF UlLIen/2-1) and 

ulbeoms2) Ff K tS Even. 


SUOPROCRAMS CALLED: QUICK, OF IND, ICOOBRA,. MIN. MAX, PPLOTE 


NaAanAAnAnANnAAnNnAANnANANAN 


seesee evener aces 
COMMON /UMIT/S TUMIT MSCREM 
OCIMENSION RIND, UIN), wamMEt 20) 
OrMEnstow vireo). c28t 0}. 6801580) C765 
Cimenston .aet( 20), LABEL (20,2) 
DATA £AG1/°LOCR’, kL OU’. ANTI’, LR P*,  BAITG’, 'OOGR’,*‘am: °, 
orm et, Uae & 
OATH RAMU, WAMGZE WAMESO NAMGTE/' Ui’,* G2B8’,' GBO’,* GIB’/ 


c 
ce 18 K 900 OA EVEN? 
e 
loooemon(n, 2) 
KOZeK/2 
RSTeM-Kes 
c 
c ComPuUTE LOCAL OUANTILES 
ce 


BO 20 Jet RST 
mger-s 
po 10 set,e 


ansPrio 
arsPrig 


ARSPIO 
BaSPIO 
arsPro 


ARSPIO 
ARSPIQ 
arnserio 
ARSPIO 
arsPio 
ARSPIO 
ansPio 
arsPlo 
ARBPIO 
ARSPIO 
aRSPTO 


anePrio 
OaRSP16 
ARSPIe 
ARSPIO 


98 


BMINTONI 1) 
BO GO 341.8 
Ezteli} © CFI Noten. mop1,exuii)) 
oo Continue 
WRITE NUNIT, 92) EXU,ENTO 
02 FORMAT(//TI10.° WU *,38.8 
° 4T10.°300U)", 3N.81F 
RITE(HUNTT, 900} 
CALL DESTAT(Y.W, WAME. LAG. LOUART,.0.0,528, 080,576) 
wuLe? 
Gisdecicaseinase-t) 
TCASE(ZsLICasetnase) 
COMPUTE RAW SPACINGS (WKICLITTLE ©} 
CALE CTOPOL ON, U, NOP) WK, BPCFAC AVLE) 


e.9ud, 
wmyssst 


POR FPERICOOGCAAM, CHECK FOR SINUSDIOAL BENAVIGUR 


nan on 


IP(torprT @0 Of GO TO 98 
CALL MAXN( WK! Oo, OmMAKY, LNOO} 
OPCONI INDO /TIOSRSO 
WRITRINUNTT. 903} Omany, OF 
O08 84 FrINOO.N 
OPPONISI/TIOCHEO 
VOreIPLOATITI-O. BI/FLOATIN) 
WRITEI NUNIT. 304) 1,U0P,0P 
CONTINUE 
CONTINUE 
C4 COMPUTE AND PLOT WEICHYED SPACINGS FOR CASE ICASE 
avicro 1: oO 
0O 100 1 * 1.N0 
wae2t1) © ForNCciULiol}, INUL} 
AVACFO + AVEGCFO © ALOG(wWR2I1)) 
100 «CONTINUE 
AVLGPO © aviGro / FLOATING) 
CALL WSPACE{ WHS .D NOP) WK). WK2.U, Sw, aviwn) 
SwWiG «© ALOGI Sw} 
WRITE( NUNIT, 902) AVLE, AVLWK AVLGOFO, Sw, SWiG 
c PLOY CUMULATIVE WEIGHTED SPACINGS WITH D* AND D- 
CALL KSOID.U. NOP! OM, 
CALL FCOORA( DP, SHIP? 
CALL FCODEA(UP anir? 
CALL PCODEALOM, BHIFT. 
CALL PCODEAlUM. BHI F7. 
Lasrotaiercaselt) 
Lasro(s jsicase(2) 
CALL PPLOTIU,D.NOPT,BLK.0..0.,1.0,81K,2,0.,1..9,0.,1. ASTER, 
a WAMU .NAMCWS.LABIO. (ASS? 
WRITE I MUNTT, 200) 
200) «FOMMmATI//////) 
800 FORMAT(///, TSO, FULLY NOM-PARAMETRIC ANALYSIS‘/TEO, 31(tme)) 
901 FORMAT(//T38%, ‘SUMMARY OF AR PARAMETRIC SELECT ANALYSIS '/ 


. TIS, 401 9M-1///} 
902 FORMATI(///,720,°A¥, LOG SPACINGS’ TAO, 
° ‘AV LOG Ww SPACINGS’ .T63,°AV. LOG HYP. FO-, 
* TAS, ‘SIGMA ZERO’, T103.'°L0C SIGMA 2ERO'// 
° .T20,8(8e,618.0)) 
903) FORMAT(/,10X. "MAXIMUM JUMP FOR OULNTI FUNCTION «= , 
** CORRESPONDING TO Of(U) © °,F12.8,//,14",°2 vu’, Fi 


oF, VOR, BIVGIM- 2) 
goa PORmMATI NOX, 16,310.78. 3,3K.F18.6) 
RETURN 
eno 
SUBROUTINE QUICKIN, T) 
CaeseeesesesattseseeessssesessessesaseseneantsGuacessteegaeessetesaceeseaes 
QUICK SORT THIS ALGORITH 18 ALSO REFERRED TO AS A PARTITIONED 
EXCHANGE SORT EXPECTED RUNTIME 1S PROPORTIONAL TO NOLOGZIN) 
ALTHOUGH THE WORST CASE 18 PROPORTIONAL TO Wee2. 
REFERENCE: CONALD €. KNUTH- THE ART OF COMPUTER PROGRAMMING VOL 3. 
1mPuT 
4%, > VECTOR TD BE SORTED OF LENCTH W 
oureur : 
x SORTED VECTOR 
SUSROUTINES CALLED : NONE 
eereesesesteesate erase etessseeeeseesseseneseteasesessestesseusasresnaens 
Rear. TIN), 
INTEGER IP ,CV(16),14V096) 07, TOP 
tvirjyas 
tv(ijen 
Pst 
1o IF(IP.LT.1) GO TO 76 
AG LOCH IVC EOS Ev CLP eS 0T.83 CO TS 20 
Go TO 26% 
20 1PelP-4 
co To 10 
28 LPsiv(IP)-1 
1uPeivire) 
yertipue) 
30 IFC CLUP-LP) LT.2) GO TO 45 
LPeLPot 
IF{T(Le) LE.¥) GO TO 30 
TErUPLeT (LP) 
38 TRC etuP-ch).L7.2) GO TO 40 
TuPelyu 
IF{T(IUP) GE.v¥) GO To 38 
TULP PETE TUR) 
GO To 30 
a0 1UPstTuP-) 
as TEP sey 
IPCC TUP-LV( TP) ET. CiVvGtO)-10P)) Go To $B 
co To 60 
SE LUC LPO PCL IP) 
Iv(sPotperue-s 
LVEIP Der UP ed 
co 70 70 
GO LVCIPO1)TtUP et 
TOL IPerpeRvire) 
TVET pe ruP-1 
70 IPerTeeos 


annnnannan 


76 RETURN 


SusaduT Merman 


Ceessecs seene 


SVSAOUTINE TO PRIWY SIGMIPFICANY RELATIVE MINIMA (AND THEIR 
twOtces) oF Caritt. . avin) SIGNIFICANT RELATIVE 
MINIMUM TS CEFTTHED TO OE A RELATIVE MINIMUM WHICH 18 AT 
ROAST CHC (USUALLY 1/SAMPLE SIZE) LESS THAN THE PREVIOUS 
valug. 
1eP UT 

m.CarTit),. carina), 
outeur 


mint - 19C8n OF OVERALL MINIMUM OF CAT 

min? IMNOEK [OTMAR Than MiNi) OF Smarr est? 
STCHIF ICANT RELATIVE MINIMUM {IF THERE ARE 
any) GR INDEX OF 2HO SMALLEST Cat VALUE. 


SUCROUTINES CALLED : MIN 


NnOaANnANnANnANnANnANNANnAAnA 


eeeaanecesasruesusese 


BasPic 
1@ 


ie 
aasrio 
arsPrio 


ARSPIO 
OarsPig 
aRsPlg 
19 
ansPilc 
arsPig 
aasPIO 
RsPio 
ASPIO 
ARSPTO 
ARSPIO 
10 
RSPITO 
ARsPia 


ars 
ARSPIO 
aansPio 


aasPrig 


BASPTO 
° 
arsPio 
aRrsPelg 


arsPio 
arsPie 
ansPte 
ARSPIO 
ansPio 
aRsPio 


ARSPTO 


100 


a a 


am Tai al 


nnn 


€ 


c 


c 


AANANANnAARANAANKANKnAAAn 


aMinmeomtt) 
oo 80 1t1.8 
e th » OF twOlom. wort. enuli)) 
oo CONTINUE 
WRITET NUNIT. 02) BRU. BELO 
82 Foamatis/T10.° wo + ,3B 
. 7710. °20tUIs. 30.8 
WAITEL NUNIT. 800} 
CALL DERTAT(Y ©. MAME. (8, LOUART,0,0.K26,"80,278) 
Kaseetmure? 
Weasecsieric 
Icasedzpevic 
COMPUTE Raw SPACINGS water tTTLe Oo} 
CALL CTOPFOION,U,. NOP) WK) SPCR AC AVL EK} 


i), 
mr.ssst 


POR PERIDOGERAM, CHEEK FOR SINUSOIDAL BENAVIOUA 


IP(1OPT €0.0) GO To 98 
COLL MAR[ WK) NO. OMSEU, TNDO} 
OPrOn( INDO) /TIOSISO 
WALTEINUNTT,. 803) Omany, OF 
OO 84 18tNDO,N 
OPsONTTD/TIO®ORSO 
vormr(Proatii)-O. 6 I/F LOATIN) 
WRITEI NUNIT, 9063 1,U0% 07 
96) «CoNTtinu® 
8 CONTINUE 
COMPUTE AND PLOT WEIGHTED SPACINGS FOR CASE ICASE 
AVLGFO © oO 
oO 100 1 + 1.N0 
WRITT) © FPORFNC(ULIe1). IMUuL) 
AvVLCPO © AVLGFO © ALOC(wK2{1)) 
100 «CONTINUE 
AVLGFO + AVLGFO / FLCRTING) 
CALL WSPACE (WHS 0. NOP! WK! , WK2,U, Sw, Avion) 
Swit + alocisw) 
walTEINUMIT, 902} AVLR AViWK AVLGFO. Sw. Swi 
PLOT CUMULATIVE WEIGHTED SPACINGS WITH De AND D- 
CALL KSOID,U,NOP1, DM, UM. OP. UP) 
Catt FCODEA(OP .BHiF? 4) 
Cart PCOoDEaluP Mtr? 4) 
Catt FCodea(om tr? a) 
CALL PCODEA(UM, SHIFT. 4) 
Lasiolse)escase(s) 
ry 1019 )etcaset2) 
CALL PPLOTIV,O,NOPT.BLK,O.,0.,1.0,81K.27.0..9.,1,0..1., ASTER, 
ry WAMU, WAMCWS.LABIO, LABS) 
WRITEIMUNTT, 200) 
200 Fommatis/////) 
900 FOoRmaATI///, TSO, °FULLY Now 
@0t FORMATI//T3%. "SUMMARY OF 
° TIS, 40 ImM-)///) 
002 vrormaTt///,t20,‘av. (OG ACINGS’,T40, 
. L0G w. SPACINGS’, T63,°AV. LOG HYP. FO-, 
° “SIGMA ZERO’, T1039, ‘LOG SIGMA ZERO‘ // 
. .T20,8(8",618.81) 
903 FORMATI/. 100, "MAXIMUM JUMP FOR OUANTILE FUNCTION © °,#12.6, 
** CORRESPONDING TO O1U) © °.F12.8,//, 98m, °4 ,16K,°O°, 
o/, Om, BC tm-)) 
906 FORMATE ION, 18,9K.F78.3,2K,F98. 8) 
RETURN 
eno 
SUBROUTINE OUITCKIN, TI 
Bee ROTHHSCESFSHSPHESHSESSSSSSTHSHSTHERETSTSTSESSSTHESReseeTEeeseeseeesesesseuasses 
QUICK SORT THIS ALGORITH 18 ALSO REFERRED YO aS A PARTITIONED 
ExCHANGE SORT. CTEO RUNTIME 318 PROPORTIONAL TO NeLOG2ZIN} 
ALTMOUGH THE WORST CASE 5S PROPORTIONAL TO wHee2. 
REFERENCE: DONALO E. KNUTH- THE BRT OF COMPUTER PROGRAMMING VOL 3. 
Iweut 


METRIC ANALYSIS’ /THO,31( tHe) ) 
METRIC SELECT ANALYSIS ‘'/ 


X,N : VECTOR TO BE SOATED OF LENCTK Ww 
ourpuT 
& : SORTED VECTOR 
SUBROUTINES CALLED : NONE 
ecees eevce 
Rear T(M),Y 
IMTECER IPL UCVENGS.TVE16) LP. tue 
Lvdupes 
Ivtt)en 
THe 
to rF(se.4T +) CO TO 78 
1S TFCOTVCIPD-EVETP PD. LT.9) GO To 20 
co to 28 
20 1063P-4 
Go Te 10 
25 LPrLVC IP): 8 
1uPesvE IP) 
vevTCrurP} 
3o TFC (TUP-CPP 0.2) CO TO 48 
LPsLpes 
TFCTOCLP).LE.Y) GO TG 30 
THtuPleTiLe? 
38 tri(lvP-t]).Lt.2) CO TO 40 
TUPtSUPet 
tec TisuP) GE.v¥} Go To 38 
TEAPDETELUPD 
co To 30 
ao 1UPeTUP- 8 
as Ti lUPley 
TFECTUP-LVC IPI AF CI¥EIHP-tUP)D} GO TO BB 
co to 60 
SS UVC TPO) ouvETe) 
Iv(t tieiv 4 
AVCTP DS LP ot 
co To 70 
6O LV¥(1POT EE TUR Os 
TuCrPorpetv( ips 
PvGtP)esuP-9 
70 [PetPer 
co ro 18 
78 aeTuan 
eno 
SUGROUTINE RELMNICAT,.. DEC, MIN), MINZ) 


aecses eee ences 


Sesecsveouseceneseseese 


SUBROUTINE YO PRINT SIGMIF ICANT RELATIVE MINIMA (ANDO THEIR 
Imorces) OF Cart(t),. cat(m) A SIGNIFICANT RELATIVE 
Minimum tf 6 twao 790 ATIVE MINimum wHICcH 
Least OC [USUALLY t/SAMPLE SIZE) LOSS Thaw THE PREVIOUS 
va.ug. 


rmeeuT : 
m.CarTt(y),...,catim), oe 


ourput : 


mtmt : INDeEr oF © ALL Minimum OF Cat 

mind? : INoas (OTN THAN Mimi) OF SMALLEST 
SIGMIF ICANT RELATIVE Minimum (17 THERE ARE 
ANY} GOR INOEE OF 2ND SMALLEST CAT VALUE. 


SUBROUTINGS CALLED : MIM 


aagPig 
a ie 
aasPio 
aasrio 
aasPio 
aRarig 
aRSPLO 
aaePre 
anePio 
anePre 
arsPrie 
arerPre 
aRsPrle 
aansPio 


aRrsPio 
aasPig 


arsPle 
SPIO 
sP19 
aansPlio 
ARSPIO 
aasPico 
ARSPTO 
aRSPI1O 
aRSPIO 
aasPio 
aRSPIO 
a@ansPio 
ARSPTO 
ARSPIO 
@RSEPIO 
aRSP1g 
ARSPIO 
@aRSPIO 
aAsPio 
aasPtgo 
@RnsPio 
ansPio 
aRsPig 
aAASPI1O 
@RSP10 


100 


4 


€ 
€ 
€ 


1° 
2° 


nan 


an 


AnnannonannanaAnnanannnannnaAaAnaanananaan 


100 


nanan 


OimMeNSION CAT( HM), RELL 20) Minti 20) 
O4Ta NOUT/E/ 


FIMO PIAST Two mimIMA : 


COLE MIMI CAT. M, Cmtier, mim) 
CaTiMINi}se. 6 

CALL MINI CAT, M, Cmim2, MIN2) 
CATIMING A CMITNY 


FIND SIGMWIFTCANY AMLATIVE MINiMa 


t#tm.48.2) GO TO BO 

£40 

Meta: 5 

BO 70 1127, mmMt 

ecstatil) 

ioltce Ge Cari l-1}-08C) OA (CC GE. cati1e+))} CO to 20 
17(1 @0 mint) CO TO 20 

heuer 

MWrtkrece 

Minimisded 

walteimour. 10) &.1.6€ 

FORMATI ION OmMIMIM( 12,466) © (13,28, PNVALUE ©, F718 7) 
CONTINUE 


OlmG mlM2 


tet. €0.0) GO Ta 80 
CMLL Mim( ROL. L, CMINZ, MINZ) 
MINZeMIMINIMI M2) 


CONTINUE 

WRITE (NOUT. 60) MINT, CMIND mMIN2, Cmluz 
FORmMAT( / ION, FHMEMI © (19,20. t2MCBT(erNy) © . 618 F/t1OK, 
VTi. © TTL BLP ANCAT(MIN2) © 718. 7) 


ReTURN 
eno 
SUBROUTINE 


SUGROUTING TO PERFORM SUBSET RAECRESS ION 


TWPUT 
@ : COVARIANCE MATRIX TO GE ANALVYIEO 
NOIM : DIMENSION OF & IW C i Rocram 
x > NUMBER OF INDEPENDENT VARIABLES IN FULL 
MOOEL 
NAME > NAME(T) «VARIABLE NUMBER OF VARIABLE OW 
1tm O}aCQNAL OF a 
KTOP ; IMDICATOR 10 MEANS SHLREG AUTOMATICALLY 
CHOOSES BEST MODEL, KTOP « * FORCES BEST M 
VARSABLES INTO MODEL) 
KTOP » -1 MEANS TERMINATE AFTER FIRST DELETION 
1f RES VAR .LT. 0.100 OR TERMINATE LF RES VAR 
LT.9.001 
1P KYOP Ce.99, CHISO VALUES FOR Ps. 90 ARE USED AND 
KTQP-100 15 USEO TO DETERMINE OTHER OPTIONS. 
IF KTOP GE.199, THEN CHISO VALUES FOR Pt. 80 ARE USED, 
aS 15 THE VALUE OF KTOP-200 
L : SAMPLE SIZE (IF & 1S NEGATIVE, HO PRINTING 
1S OONE) 


ourrur 
KP: MUMBER OF PREDICTOR VARIABLES tu CHOSEN 
MODEL 
IND : VECTOR TELLING WHICH VARIABLES (ACCORDING 
TO MAME) HAVE Gh CHOSEN 
cor : VECTOR OF COEFFICIENTS FOR CHOSEN MODEL 
JER . ERROR IMDICATOR (1 MEANS ALL DIAGOWAL 
ELEMENTS OF A.LT.1.8-8 AT A CYCLE. O 38 
NORMAL RETURN) 


SUBROUTINES CALLED : NONE 


DIMENSION AC MOTM,WOIM) WAME(K) ,COP(K), TMOIK), 
1CHISO(62), CBlLI00}, 1NB11 100) 
DATA wouT/6/ 


1OPTPr) 

LLOPT OL 

tF(L.6t.0) GO To so 
werk 

1OPTP1O 

CONTINUE 

RTOP1eKTOP 
1PiKTOP.GE.98) GO TO 2 
CMLL CH1IP(1, CHISO) 

co To 6 

1F( KYO, ce.198) GO TG 3 
Cat CHIP(2, CHISsO} 
XTOPeKTOP- 100 

co To 4 

CALL CHtP( 3, CHISO) 
KTOPCKTOP-200 

CONTINUE 


Oo 10 1e1,K 
eCBlbieco. 
rwOrlt)»o 
rEReo 
te(10"1TP.£0.¢) GO TO 81 

WALTE( HOUT, 69) K RTOP 

PORMAT(//, 10K, "REGRESSION ESTIMATION STAGEWISE SUMMARY’ /10x, 
1° MUMBER OF VARIABLES IN FULL MODEL «© °,13/10R, 

1°MAXIMUM NUMBER OF COEFFICIENTS © *,33/) 

wairtetwour, 1) 

PonmaT( ion, 78(1m-)) 

cowtinue 
Wenel 
PORMATIIOR. 


€10.3,2")) 


F268 638 
Riso 
FOr)! &-6 
REeo 
aPeo 


Stace 1 AUALVZE MATRIX & \ 


atoreo. 
wmibrzeese 


Re See ls 


arsPig 
arnerio 
ansete 


aasPio 
arselto 
aaselo 
e190 
ansPio 


aRsPlg 


aasPrio 
ARSP1O 
arnseio 
aRSPIO 


aRsPlg 
ARSPIO 
ARSPIO 
@RSP19 
aRsPio 
@aRSP10 
arsPio 
aaseto 
ARSPIO 
aRsPio 
ARSPIO 
ARSPIO 
aRsSP?10 
ARSPtoO 
ARSPIO 
arnsPrto 
ARSP10 
ARSPIO 
ARSPIO 
a P10 
arsPig 
arsPig 
ARSPIO 
AnsPlo 
aaseio 
ar 9 
ansPie 
aRrsPic 
arspig 
ansPie 
a @ 
° 
ansPio 
aansPIo 
arsPto 


ARSPIO 


ansPie 
ansPig 
aRsPie 
arsPio 
AREPIO 


aasPrig 
BRSEPIQ 
ARSPIC 


101 


vMarso. 
amine® 
fee 


RTOL COUNTS NUMBER OF DIACONALBS.6T. 


aTOLSO ansric 
OO 110 IKTet Kh a ie 
Veal int. wpeatm, rat} aasrie 
ansPeie 
CRECK BP VARIABLE IKT 18 1% MODEL anselo 
ansrio 
Flv 47.0.) GO TO 108 aasese 
COLIKT) OO, ansPio 
PwOTLIKT) 20 ARSPIO 
co TO 106 ansPio 
COLIKTI ALERT. Mw} aRSPIO 
DwOrlTaRTied aRSPIO 
KRPeKPet aasPio 
ARSP10 
CHECK DIAGONAL ALIKT,IKT) : ARSPIO 
ARSPTO 
106 IFC MLIKT,1KT).LY. TOL) GO TO 110 ARSPIO 
RTOLSKTOLS! ARsPlo 
vev/a(iKt, 1K) arselo 
iPtv 47.0.) GO 76 60 aasPic 
RIOTERTOTOY AAsPIO 
co TO 80 ARSPIO 
AABPIO 
NO VMIM, WMI: ARSPIO 
ARSPIO 
TFC TABS V)-aBSl¥MIN)).GT.O.} GO TO 110 aasPio 
vMINGY ARSPIO 
MMIMEIKT ansPro 
co TO 110 ARSPIO 
ARsPi9 
NO VMAN.NMAN : aRnseio 
aRrsPio 
1F((V-veaK).L€.0.) GO TO 110 ARSP1O 
vee v . ARSP1O 
weaned KT ARSP1O 
110) «CONTE NUE aRSPIO 
ARSFIO 
ARE ALL DIAGONALS. LT TOL > ARSPIO 
ARSPIO 
IFUKTOL.€0 ©) CO TO 999 aRsPio 
ARSPIO 
CALCULATE AND PRINT CAITERION : ARSPIO 
ARSP1O 
PMD KP ed aRSP10 
RYDER KP aRsPio 
1F (KT! €0.0} CO TO 120 ARSPI1O 
FVeCHISQUKTS) ARSPI1O 
co To 128 
wed €-10 
TIF(MBS(VMIN) OPMT)/IF2EA( Nh) ) 
TIE(RTOTOPHII/IF teat Mw)? 


1FIMMIN. EO ©) CO TO 130 ARSPI1O 
KOTEMAME (MIM) ARSP1O 
co 70 140 ARSPIO 
xO100 ARSPIO 
IF (wman. €0.0) GO TO 180 ARSP1O 
KOZ: NAME [MAX ) ‘ ARSPI1O 
60 To 160 ARSPIO 
e020 aRrsPioc 
IF (1.20.0) GO TO 170 ARSPI1O 
KODeMAMEeI KE} ARSPIO 
co TO 180 ARSPIO 
170 «KO3"0 @RSP1O 
ARSPIO 

@kAIKE : a@nsPio 
ARSPIO 

180 PHIOFLO2EKE SPIO 
PMTTSOMT OC? ARSPIO 
PUTZEPNIO-2 ARSPIO 
PMIIH2 eFLOATIKP)/FLOAT(L) ARSPI1O 
TICALOCI AIM, MI IOPMIS sP10 
Tar(vmManePMItt/(2. 9am, m)) aRrsPro 
TErivetwePMbai/si2.emin, wd) ARSPIO 
ARSPIO 

spto 


. 
PRIMT COEFFICIENTS ANG VARIAB 


1F(10"T° £90.60) GO TO 52 ARSPIO 
WRITE(MOUT, 190) ARSPIO 

190 FORMAT( 10K, 12NVAR IN MODEL, 10", SHCOEFF) ARSP1O 
$2 CONTINUE ARSPI1O 
1F(KP £0.0) GO TO 1846 ansPioc 
ARSPTO 

ARSPIO 

P10 


TFC 1OPTP £O.0} GO TO 101 

waiTe(woury, 192) 12,ce(t) 

FORMAT( 10" ,112, 62,710.86) 

CONTINUE 

1F(10PTP. £0 0) Go Te 8S 

wRITE(NOUT, 108) 

PORMAT(ION, OX, SMCRIT 200,6K. MHCRIT CEL, 7K, 7H 

1T7R, PMFPE 200,70. 7HFPE OBL) 

WRITE{ HOUT, 1661 12,71, 73,740,768 

Forwmat( ion, 61F12.6,2%)) 

WRITE( HOUT, 1875 

PORMAT(IN, OR, IK, THVAR BOO. SX, THVAR DEL, OH, SHVAR LAST, IK, ansPro 
19HNGO. PREOS, SH, PHRES VAR Tueo. CYC) ARSPIO 
WRITE{NOUT, 188) KOZ, KDI. KOS, KP ACM) KE aRSPTO 
ronmat( ex, 4112,2%,F710.8,3%,110} ARSPIO 
WRITE{NOUT,1) sPto 
CONTINUE ARSPIO 

ARSPIO 
ARSPIQ 
STAGE 2 : CHOSE PIVOT OR TERMINATE 


200 «IFLTKTOP GT.0). aNd [KP LT. KTOP) AND. iwman.€0.0)) GO TO 919 
CPCI ATOR CT. 0) and. [KP -RTOP) And. i wmax we .O}} CO TS 300 
IPCIRTOP.6CY 0) and. (KP. 8@.KTOP}) GO TO 380 BRSPTO 
aRSPIO 
automaric : aRsPro 
ARSPIO 
Cmece FORM OBLETE : ARSPIO 
ARSPIO 
ARSPIO 
CHECK RYOPS-1 CONDITIONS - f aARSPIO 
AASPIO 
iF(atoe we.-1) Go TO 228 ansrig 


Anannanan 


TFL LRP LT. KE). BND Fain. w).LT.0.9)) CO TO 280 aREOTO 
1Ftaine,@).b7.0@.0091) GO TO 2238 


CuRCR OELETE COMBSTIONsS FOR KIOP WE. -t 
226 IPCC MOR (ves). Ce. F2eAI he WI sR) OR. (KI. BO. wmMInN)) GO TO 230 


vas 


aan ann 


22B KV 
60 %9 240 


Cweck a00 


nnanan 


230 «fF (RTOY LE Fie(Aiw. MI sPHI)) GO TO 280 


ves 


ann 


300 Ktenman 
340) «CONTINUS 


stace 3 . PIVOT 


annannn 


PlvOtet. sMURT, KI) 
DO 380 URTEI, A 

BSO ALK IKRTDEMEKT, VET) OPI VOT 
20 370 IKTe1,h 
TeTIKT.£0.K1) GO TO 370 
TEMP CALIRY KE) 

DO 360 JKTs1,4 

180 ALITY, JRTPSALIRY, SKT}-TEMPOMI KT, JKT} 
ALIKRT Mt) e-TEmmPer vor 

370 «CONTINUE 
ALKE, MUP OPT VOT 


c 

c METURN TO STAGE 1 

c 
wCeKces i 
1P1KE.LT.100) GO TO 100 

e 

c TEAMINATIONS 

© 
WRITE( MOUT, O13) 

912 FORMATI/,. SOX, "TERMINATION OF SELREG SINCE KC.GE.100°) ; 

6O To 280 

c 

c 


918 WAITE(NOUT, 920) 

O20 FOMMAT( /. ION, “TRAMINATION SILWCE ALL REMAINING VARIABLES HAVE © PAR 
TT2AL VARIANCE’) 
co TO 280 


c ‘ 
c 
999 WAITEINOUT, 1000) 
1000 PORMAT(/, JON, “TERMINATION SINCE ALL DIACONALS 147. TOL "3 * 
igaer 
200 KTOPERTOR! 
Leccopr i 
agtuan ° 
c 5 
c 
eno 
SUBROUTINE SLMNIRP, IND. COF, ALPHA, BETA, WORDAR, WORDMA) 
Cescee vesesese ARSP1O 
c ARSP1Q sf 
c SUBROUTINE TO TRANSFORM MINED SELECT SCHEME PARAMETERS ARSPIO 
c RP, INO.COF TO FULL MIRED SCHEME PARAMETERS HORDAR.NORDMA, ARSPIO 
c avona, BETA ' ARSPIO 
c ARSPIO 
c Tweur ARSPI1O 
c xP: Un OF LAGS IN SELECT SCHEME arsPlo 
c imo(t), .1WOIKP) ; LAGS (AR Lacs «+, ma Lacs -} ARSPITO 
c Coriid,....COF(KP) : COEFFS CORRESPONDING TO Lacs ARSP1O * 
c ARSP10 ' 
c OuTpur : ARSPIO 
ec NORDAR,NOROMA . AR ANO MA PART ORDERS ARSPIO 
¢ ALPHA(T),..., ALPRA(WORDAR) : AR CORFFS ARSPIO ‘ 
€ eeratih, .CzTat(woromal : ma cogres arseto 
. ARSPIO 
Cecece ecee ease eves ence as0e eee aRsP1o 
c ARSP1O * 
OIMENSION ENDIKP),COFI KP) aLPmal)), SETA(1) ARSPIO . 
c ARSPIO “ 
c FIND WORDAR.NORDMA - ARSPIO 
e BRSP1O . 
WORDARO ARSPIO . 
NORDMALO ARSPIO 
IF (KP.@20 ©} GO TO 99 ARSPIQ 
DO 20 105, KeP ARSP1O 
trerwots) ARSPI1O 
1F(1T.GT.0) GO TO 10 ARSP1O é 
IPC-11.6T.NORDMA) NOROMAr-1I1 ARSP1O 
So 10 20 aRsPilo 
10 «CONTINUE ARSPIO . 
IPC11. GT. MORDAR) WORDARETI ARSPIO ' 
20 «CONTINUE ansPioc 
c ARSPIO 7 
e ZERO CUT ALORS AND BETA : ARSPIO . 
¢ ARSPIO : . 
TP(MOROBA.£O.0} GO TO 40 ARSPIO 
00 30 1 RORDAR sP1O0 
acemat(i)9o.0 ARSPIO 
30 «CONTINUE ARSPTO 
a0) «CONTINUE ARSPIO 
1FtwOROMA.£O.0) GO TO 60 ARSPIO 
00 Bo 121, NORDMA ARSPIO i 
GETall}so,0 aRSPLO 
so continug ARSPIO 
$0 continue ARSPIO « 
c ARSPTO F 
c PUT COEFFS INTO ALPHA, RTA . ARSPIO . 
e ARSPIO 
00 80 101,00 ARSPIO 
tresmein) aRsPloe 
wtit CY.e) Ge Te 76 aRsPto 
dse-dt 
airslccerts) 
18 
70 «cenTinus 
OLPmaiti)scorti) ARSP1O ‘ 
60) «cONTInuE ARSPIO 
c aRSPTO 
e Pinte - anseie 
€ aRsPic 
99 CONTINUE ARSPIO 
aeTune ARSPIO 


ane ansPio 


SUSUCUTTWE SVevi seer WREST EY BY NOT 


Cocessaesesenevasesssesessesseenncosseccoenssennesesenes 


SUBROUTINE TO CALCULATE THe AUTOCoOVaRTaANe MO.Mith, .., Rhee) 
PROM THE CORRESPONDING SPECTRAL DENSITY ©) 

EVALUATES AT THE NFREGS BOUALLY BPACED FREQUENCIES 

GRIWEEN O ANO TWOP]. 
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inpur 
NPREOS CINTEGER WHOSE LAACEST PRIME FACTOR I8 
© 23). m. seeciit, SPECI NFREOS) 


outeur 
wmo,arty). ~ Rim) 


AUNILE TARY 
cry. .s7T 


aasPig 
e 
aRSPIO 


SUBABUTINES CALLE vRFY 


Peewee eeesvaceceeresassersascasnenesnneceneasteeeseessse 


nnannnannnnnnnnnnanan 


DIMENSION SPECINFREOS) CTINFREOS} ST( wWFAROS) RIM) 
c AASPIO 
PACe(S CATAN[1 OF )/FLOATINFREOS) AASPICO 
c @RSP10 
OO 1 tat wWFEREOS 
cT¢rpespec(i) 
STC LI«O. 
CALL FFTICT ST, WFREOCS NFREOS WFEREOS,-1) 


ROUFACECTI1) 
00 2 1let.m 
2 RiLperacecr(iel) 


RETURN 
euo 
SUBROUTINE SPPLTISPEC NFREOS FFREO OELF.10PT) 


SUBROUTINE TO PAINTER PLOT THE SPECTRAL OCUANTITIES 
SPECI YT}... .  SPECINFREOCS) WHERE THE FREQUENCIES ARE 
GIVEN BY PFREQOCIJ-1)°DRELF, Jel, NPREOS. 


LePuT 
MFREOS PFREQ OELF 
srpecti)...., SPECINFREOS) 
10PT : 1 MEAMS LOC SPEC PLOTTED, SPEC PRINTED 


SUGAROUTINES CALLEO : MAX, MIN 
ARSPIO 


weeetoceneneasaseseseeseeseresesseenssesee saeeseeecsesen ARSPI1O 


naaanananannannan 


OIMENSTON SPECI NFREOS) .ALI61) 
OATA WOUT, BLANK, COT, K/6.9N (IK. uHe/ 


scare 


nan 


Ie(SPECc( TE. LY 4.@-8) SPECI 1) =sPect2) 
1f(10PT wWE.1) GO TO § 
smaxeavocispec(s)) 
SMINTSMAX 
OG t 192, NFREOS 
Cracoc(srec(s)) 
PEC. CT. SMaAK) SMAREC 

’ TF(C. LT. SMIN) SMINEC 
SRANGESCSMAN-SMIW 
co vo 7 

S COth MAX(SPEC NFREOS ,SMAX, IND) 
CALL MIN(SPEC NFREOS .SMIN, IND) 
SRANGCE*SMAXK-SMIN 

7 CONTINUE 


c INCREMENT FOR 2 PACES 


wiwres 
TF (MPREOS CT.120) WINTs(WFREOS/120)01 


IMITEIBLTZ2€ BNO PRINT MEADING : aRSPIO 
ARSPTOQ 

WRITE (NOUT, 10) 190 

10 PORMAT(/, 1ORK, 1%, SHFREQUENCY 6X, GHPERIOD, BX, ARSPIO 
VGHSPEC/ ION, BAC IM-), aK, OFC IM.) ARSPIO 

DO 20 145,61 ARSPIO 

20 BLED) FOL ANK 


nan 


PLOT 


ann 


IFUSMANGE .GT.1 €-28) GO TO 8 
wALTE(NOUT,. 9) 
9 FORMAT(/, 10%, RANGE.LT.1.€-28 IN SPPLT’) 
co To 99 
@ CONTINUE 
TwOP1s& SATANI1.0) 
DO JO 1:1, NFREOS .NINT 
PRON (PFREOSFLOATII-1)°ORLF)/TWOPT 
PER 
IF(FRO.GT.1.£-10) PERSS./PRO 
SSesPec(i) 
IF(TOPT €0.1} SS ratoc(ss} 
Ce2. (1 (S8-SMIN) /SRANGE)- 8} 
Jt30 21004. 501.8 
aLlsyen 
WRITE(MOUT, 28) FRO, PER, SPEC(1), aL 
25 FORMAT(ION, 710.65, 20,F10.5,2N,F10.8.6xK, 61A0) 
ALC J) eOLanK 
30 «6<CONTINUE 


¢ 
c 
99 CONTINUE 
METUAN 
eno 
SUBROUTINE WSPACE (WHS CWHS NOPT, KS. FOOO,U,SI1GO,AVLWK) ARSPTO 
Coccccccavecseseves Pee easosanenvesscarscesensesessaes ARSPIO 


c SUBROUTINE TO COMPUTE OLU). CUMULATIVE D'S, AND SIGMAO 
€ POR TME MODEL OC{U) *MUcsSICMAsoo(U) 
c ImpuT : ARSPIO 
¢ FO.NQ : VECTOR OF LENGCTM NO CONTAINING FOU} ARSPIO 
c FOOO : HYPOTHESIZED CENSITY OUANTILE FROM FOFNC ansPic 
c Uo: VECTOR OF LENGTH NO CONTAINING U VALUES AARSPI1O 
c ourpur : aRSPIO 
c WHS =: VECTOR OF LENGTH NO CONTAINING DIU) aRsPio 
c CwRS : VECTOR OF LENGTH NO CONTAINING THE a te 
c CUMULATIVE O'S ARSPIO 
c S1¢o0 COMPUTED VALUE OF SIGMAO «© CwxS( NO} aRSPIO 
c aviws AVERAGE LOG OF LITTLE BiU) ARSP1O 
4 SUOROUTINES CA aasPio 
Cececee aeve a ecescacnsesecazece ° ARSPTO 
DIMENSION FOOCO(NOP1). NS(NOP1) _UINOP1) ,wEHS(WOPT) .CwnS(NOPT) aRsP1O 
ROTMOPT-1 ARSPIO 
CwRsiiteo. aRsPro 


20 


00 10 bs 1,80 

wEGL I) © Fe@OtR? © BEIT} 
Cwesites) © Cwneii) ¢ westi) 
CONTINUE 

er 

Pw es FLOATING) 

Bt 6 1 sCwnsimort) 

02 + oF © ON 

oO 20 1¢1,R@ 
@as(ilewestij © 92 
CwrsitdeCwneti) ¢ Bt 

& © 8 © MeOGiwasii)! 
CONTINUE 

Careiworiiet 

aviws + $/PH 

8160 © +. /02 

aeTuan 

eno 


be Oe Cait Ce ee Ee eR See On ee Sree eee ee 


1 
aaserig 


ae ae 
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